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Scope of JET PROPULSION : 


This Journal is a publication of the Ameri- 
can Rocket Society devoted to the advance. 
ment of the field of jet propulsion through 
the dissemination of original papers disclosing 
new knowledge or new developments. As 
used herein, the term “jet propiilsion” 
embraces all engines that develop thrust by 
rearward discharge of a jet through a nozzle 
or duct; and thus it includes air-con. uming 
engines and underwater systems as yell as 
rockets. JET PROPULSION is open to con- 
tributions dealing not only with pro, ulsion 
but with other aspects of jet-propellec flight, 
such as flight mechanics, guidance, tel« :neter- 
ing, and research instrumentation. In- 
creasing emphasis will be given to the sc entific 
problems of extraterrestrial flight. 


Information for Authors tien 


Manuscripts must be as brief as the pro; =a 
presentation of the ideas will alloy xX. Nati 
lusion of disp ble material and « oncise- 
ness of expression will influence the | ditors’ 
of a manuscript. In te ms of Labc 
standard-size double-spaced typed p»ges, a 
typical maximum length is 22 pages of text 
(including equations), 1 page of ref«rences, 
1 page of abstract, and 12 illustrations, 
Fewer illustrations permit more text, and vice 
versa. Greater length will be acc.ptable 
in exceptional cases. 

hort manuscripts, not more than one 
quarter of the maximum length stated (or full 
articles, A qualify for publicat:on as 
Technical Notes or Technical Con ments. 
They may be devoted to new develo; ments 
requiring prompt disclosure or to comments 
on previously published papers. Such manu- 
scripts are usually published within two 
months of the date of receipt. ; 

Sponsored manuscripts are published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
publication according to the above-stated 
requirements as to subject, scope or length, @ J@CESs 
but which nevertheless deserves widespread ause: 
distribution among jet propulsion engineers, 9°” 
may be printed as an extra part of the Journal J pheno 
or as a special supplement, if the author or §- 
his sponsor will reimburse the Society for @ ae! 
actual publication costs. Estimates are rent! 
available on request. Acknowledgment of §* 
such financial sponsorship appears as 8 
footnote on the first page of the article. 
Publication is prompt since such papers are nolect 
not in the ordinary backlog. Ear! 

Manuscripts must be double spaced on one a 
side of paper only with wide margins to allow §mdiat; 
for instructions to printer. Include a 100 to 
200 word abstract. State the authors’ 10 Cor 
positions and affiliations in a footnote on the and ex 
first page. Equations and symbols may be 
handwritten or typewritten; clarity for the high t 
printer is essential. Greek letters and unusual 
symbols should be identified in the margin. If Jecent 
handwritten, distinguish between capital and f gase 
lower case letters, and indicate subscripts and 5“ 
superscripts. References are to be grouped at 000 
the end of the manuscript and are to be given 
as follows: for journal articles: authors first, 
then title, journal, volume, year, page numbers; 
for books: authors first, then title, publisher, 
city, edition, and page or chapter numbers. 
Line drawings must be clear and sharp to make 
clear engravings. Use black ink on white 
paper or tracing cloth. Lettering should be 
large enough to be legible after reduction. 
Photographs should be glossy prints, not 
matte or semi-matte. Each illustration must 
have a legend; legends should be listed in 
order on a separate sheet. f 

Manuscripts must be accompanied by 
written assurance as to security clearance 1D De 
the event the subject matter lies in a classifi tel 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility rests A as 
with the author. 

Submit manuscripts in duplicate (original 
plus first carbon, with two sets of illustra- 
tions) to the Editor, Martin Summerfield, 
Professor of Aeronautical Prince- 
ton University, Princeton, N. J. Preprints 
of papers presented at ARS national meetings 
are automatically considered for publication. 
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After receiving a B.S. degree in Mathematics and Physics 
fom the District of Columbia Teachers’ College, the 


author was employed by the Aerodynamics Section of the 
National Bureau of Standards. In 1947, he joined the 
Aerodynamics Department of the Cornell Aeronautical 
laboratory where he carried out research in propulsion 
and high temperature gasdynamics. He received his 
Ph.D. in Physics from the University of Buffalo in 1955, 
and is currently head of the Aerophysics Laboratory, 
\erophysics Department, Space Technology Laboratories, 
Inc. 


Introduction 


ERODYNAMIC heating problems associated with 

hypersonic re-entry of ballistic missiles and satellites 
into the atmosphere have been greatly complicated by the 
uecessity for considering new mechanisms for heat transfer 
caused by radiation and atomic surface recombination 
phenomena (16).! These processes, previously unencountered 
in aerodynamic phenomena, arise because a high speed 
mentry vehicle is surrounded by a sheath of highly luminous, 
ligh temperature air composed of a complicated mixture of 
nolecules, atoms, ions and electrons. 

Early attempts (22,27) to assess the magnitude of the 
ndiative contribution to heating indicated values equivalent 
0 convective effects and stimulated extensive theoretical 
ind experimental studies of the emissive characteristics of the 
ligh temperature air components. This article describes 
reent progress in the determination of the radiative properties 
ii gases (primarily air) in the temperature range 3000 to 
(000 K. The review and bibliography are intended to 
erve only as a guide to current investigations and conse- 
wently are not complete. Electric arcs, for example, are 
mown to produce temperatures in, and well above, this 
nnge. A comprehensive survey of recent advances in the 

ermination of radiative gas properties must necessarily 
clude mention of are studies (10,11) as well as investigations 
{flames (48) and explosions (25,48). 


Determination of Radiative Properties of Gases 


A gas emits radiation as a result of rotational, vibrational 
id electronic transitions from excited energy levels to 
wer energy levels. The emitted radiant energy corre- 
ponding to these transitions is distributed over a wide wave 
igth region. The total radiant intensity emitted from a 
iven quantity of gas is obtained by summing the radiant 
ntensities corresponding to each of the individual energy 
een The simplest approach to the determination of 
e radiative intensity of gases is based upon the determina- 
ion of overall emissivities or absorptivities for gaseous 
ixtures as a function of pressure, temperature and path 
mgth (33). This approach involves only limited appli- 
htion of the fundamental relations determining the intensity 
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of absorbed or emitted radiant energy,” although it has been 
developed to the point where reasonable predictions can be 


made for practical applications. 


Such calculations require knowledge of the following 


factors (17): 

1 The number of molecules present in each of the various 
energy levels in thermal equilibrium. 

2 The transition probabilities for each of the possible 
energy transitions. 

3 The frequencies for the energy transitions. 

4 Thespectral lineshape. The spectral line shape at high 
temperatures is primarily influenced by gas motion (Doppler 
broadening) and the presence of neighboring particles (pres- 
sure broadening). 

In the temperature and density range of interest for 
missile and satellite re-entry, the most important air com- 
ponents according to concentration are Nz, Oo, NOo, NO, N, 
O, O-, O2+, NO+, O* and electrons (14,24). Contri- 
butions to the observed radiant flux from these components 
can be conveniently classified as transitions yielding con- 
tinuous spectra and transitions yielding discrete, line, spectra. 
Contributions to continuum radiation can further be classi- 
fied as free-bound, wherein the electron-ion collision results 
in attachment, and free-free, wherein the electron-ion or 
electron-atom collision results in emission of energy without 
union between electron and particle. 


Absorption Coefficients for Air Components | 


Calculations of the absorption coefficient of air in the 
temperature range 2000 to 18,000 K were first carried out by 
Magee and Hirschfelder (25) in connection with studies of 
radiation effects of atomic explosions. They considered 
photoelectric absorption by N, O and O-, and absorption 
by NO:. The absorption spectra of the Nz and NO band 
systems were not taken into account. When a state of 
thermodynamic equilibrium exists, the emission coefficient 
is related to the absorption coefficient through Kirchhoff’s 
law, and either the emission coefficient or the absorption 
coefficient can be employed to describe the radiation. 

The possible significance of NO as an important contributor 
to radiation in this region was pointed out by Teller (22) 
and Meyerott (27). In calculations of absorption coefficients, 
one of the important parameters, which is a measure of the 
intensity of a spectral line, is the oscillator strength, or f 
number, which relates the actual energy absorbed by a 
spectral line to the energy absorbed by a classical electron 
oscillator. 

The major uncertainty in calculations of the NO contri- 
bution is related to the integrated oscillator strength or f 
number to be assigned to the NO molecule (27). Use of f 
values of 0.1 indicated that the NO molecule would be one of 
the most important constituents for absorption. However, 


2 In principle, rigorous solutions to the problem of emission of 
radiant energy from an assembly of gaseous emitters which 
are in thermal equilibrium can be evolved, based on quantum- 
mechanical considerations. In practice this is difficult because of 
uncertainty with regard to many of the fundamental atomic an 
molecular parameters; see for example (33). 
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emissivity measurements by Keck made at Avco Research 
2 Laboratory, using a shock tube to produce the radiating gas, 
demonstrated that NO was not an important contributor 
(21,22,38). This result was substantiated by the absorption 
measurements made by Weber and Penner in cold NO (45). 
_ This work, combined with measurements by Marmo (26) 
- and Hurowitz (28), indicated that the f values could be as 


small as 0.002. 
Meyerott (27) carried out the first extensive calculations 
of the absorption coefficients for air components. In this 
- work, which has now been modified to reflect the effect of 
recent experimental studies (28), 0.2 was taken for the 
= f number ‘for all molecular electronic transitions. 
He considered the N- first and second positive systems and the 
N.+ first negative system. Contributors to the continuous 
absorption which appeared significant were photoelectric 
absorption by O-, O, N and N; and free-free absorption by 
electrons in the fields of both positive ions and O. 


Experimental Determination of Absorption 
Coefficients 

Experimental determinations of absorption coefficients 
using room temperature gases are of limited value for studies 
of high temperature radiative properties, because many 
species make contributions to the absorption coefficient by 
absorption from excited states. Other techniques must be 
employed to evaluate contributions of excited species. 

One of the most promising of these techniques centers 
around the use of the shock tube,* and this technique has 
been used extensively by Keck, Kivel and Wentink (21,38,46) 
of the Aveo Research Laboratory. A confined region of gas 
is shock-heated to the temperature, density and thermo- 
dynamic state of interest to hypersonic flight problems by 
reflecting a strong shock wave from the closed end of a 
shock tube. In this program, measurements of absolute 
intensities were carried out as well as measurements of time- 
resolved emission spectra. 

Conclusions from these studies at 8000 K and sea level 
density were that in the visible and near-infrared regions, 
little structure can be observed. Above 5500 A the radiation 
was attributed to the N2 first positive system. Fig. 1 shows 
the measured wave length dependence of radiation from air 
at 8000 K. Significant contributors to radiation were 
determined to be NO (vy band), Ne (second positive), Ne 
(first positive), CN (violet), O~- (free-bound) and O + e 
(free-free). Positive spectroscopic identification was obtained 
for the Ne and CN band systems. One of the more signifi- 
sant results of these studies was the determination that in 
flight at a Mach number 20, at an altitude of approximately 
100,000 ft, the radiative heat transfer from a 1-ft radius 
sphere would be only approximately 10 per cent of the 
aerodynamic stagnation point heat transfer. 

Although the work of the Avco group has clarified many of 
the important questions about the absorption properties of 
air in the 8000 K temperature range (23), Meyerott (28) 
points out that additional work is needed in other temperature 
and density ranges and with different gas mixtures to isolate 
particular contributors. The relative importance of different 
species changes with temperature and density. For example, 
at 8000 K and 10~-* normal density, the relative importance 
of the N.* first negative band increases, O~ becomes unimpor- 
tant, and the photoelectric absorption from the excited 
states of N and O becomes a contributing factor. At 12,000 K 
and normal density, in addition to the Ne first positive and 
second positive band systems, the N.* first negative system, 
and O~ photoelectric absorption, the photoelectric absorption 
from the excited states of N and O, as well as free-free absorp- 
tion by electrons in the field of positive ions, contribute to the 
total absorption. At 12,000 K and 10-* normal density, 


3A detailed discussion of shock tube techniques i is presented 
later in this paper. 
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Figure 1 Monochrometer-photocell measurements of the distri- 


bution and absolute intensity of the radiation from 1 cm of air at 
approximately 8000 K and p/p. = 0.85, compared to the esti- 
mated radiation 


the molecules are dissociated, and the only remaining con- 
tributors are free-free absorption by electrons in the fields of 
positive ions and photoelectric absorption from the excited 
states of N and O. The contributions of many of these 
species are as yet poorly evaluated, and the possibility exists 
that a constituent other than those presently considered 
might make a significant contribution to radiation. 

Shock tube studies of the radiative properties of high 
temperature air have also been carried out by Wurster, 
Glick and Treanor at Cornell Aeronautical Laboratory 
(49,50) in the temperature range 4000 to 7000 K. From 
absorption studies, the gross absorption coefficient for air 
was obtained. Typical features of the spectrum observed 
were: NO @ bands at 6000 K and normal density; strong 
continuum absorption below 2600 A at 4000 K and four 
times normal density; CN bands at 3900 A and sequences 
of the second positive system of N» in emission at 7000 K, 
and high vibrational transitions of the O2 Schumann-Runge 
system. 

Such studies of the radiative properties of high temperature 
air will be of increasing significance for hypersonic photo- 
reconnaissance and satellite re-entry applications (13,28). 
They may also provide improved recognition techniques for 
countermeasure applications. 


Recent Theoretical Studies of Radiation 
From High Temperature Gases 


One of the outgrowths of the recent interest in determina- 
tion of the emissivity of high temperature air has been the 
renewed interest in the development of theoretical approaches 
which will permit more accurate prediction of the radiation 
from gaseous mixtures at high temperatures. Jarmain, 
Fraser, Nicholls and Turner (19), at the University of Western 
Ontario, have made extensive studies of molecular band 
intensities. In particular, they have carried out a ictal 
investigation of the Nz first positive system and have com- 
puted overlap integrals for the NO 8 and y bands. Kivel, 
Mayer and Bethe (22) have extended the latter results. 

Breene* (2,3,4,5) at the Aerophysics Laboratories of the 
General Electric Company is undertaking an extensive 
program to calculate diatomic vibrational wave functions, 
molecular electronic wave functions, atomic wave functions 
and wave functions for free electrons in the fields of various 
atoms. The primary objective is the development of 4 
theoretical approach which will permit the accurate pre 
diction of radiation from high temperature gas mixtures. 


*R. G. Breene Jr. (private communication). 
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Radiation Phenomena in Rocket Thrust _ 
Chambers 


Radiation phenomena similar to those encountered in 
ballistic missile applications are also of importance for the 
design of cooling systems for rocket exhaust systems. As 
emphasis is placed on achieving greater specific impulse and 
thrust values in chemical rocket engines, combustion chamber 
temperatures tend to increase to values approaching 3000 K. 
As combustion chamber temperatures increase, greater 
percentages of heat will be transferred to the nozzle wall by 
radiative processes. The problem, therefore, exists of 
determining how large a fraction of the total heat transfer is 
produced by radiant energy. At the present time, only the 
transitions corresponding to the infrared vibration-rotation 
bands make important contributions to the total radiant 
heat transfer at the temperatures of interest in connection 
with combustion chamber studies (1,33). 

As a typical example of this problem, studies were carried 
out at the Jet Propulsion Laboratory for the RFNA-hydrazine 
propellant system (29). Although water vapor was the 
only component of the combustion products capable of 
emitting energy, experimental measurement of the emissivity 
indicated that, in terms of heat transfer to the combustion 
chamber, radiation could contribute 10 to 30 per cent of the 
total heat transfer. 

In many cases, the problem of radiant heat transfer in 
rocket exhaust systems is more serious than the corresponding 
ballistic missile problem in that the combustion chamber 
pressures tend to be extremely high, between 300 and 600 psi. 
Since the density is correspondingly large, the radiant energy 
in the rocket combustion chamber can contribute a much 
greater fraction of the total overall heat transfer. 

The problem of the accurate determination of radiant 
energy contributions in rocket combustion systems will, 
in the future, be an extremely complex one, since new fuels, 
both liquid and solid, are being considered which will operate 
at increasingly higher combustion chamber temperatures. 
In many cases very little information is available concerning 
the chemical composition of these fuels at high temperatures, 
and, consequently, the problem of determining the radiative 
contributions will become increasingly difficult. 


Use of Shock Tube for Study of High 
Temperature Radiation Phenomena 


The shock tube is now recognized as one of the more 
promising of the experimental tools for studies of high 
temperature radiation phenomena. Kantrowitz was one of 
the first to recognize its importance for the creation of high 
temperatures in a controlled manner. At Cornell University, 
studies of radiation produced by high temperature argon 
were initiated, and high speed spectrographic techniques were 
leveloped (36). 

The combination of the development of the shock tube and 
the refinement of high speed recording techniques (18,34,37, 
43) has enabled investigations of radiation phenomena in 
regions which were previously inaccessible. Recent develop- 
ments in shock tube driver heating techniques involving the 
use of plasma generators for inert driver gas heating promise 
ichievement of enthalpy and temperature conditions even 
higher than those currently obtainable by hydrogen and 
combustion driver methods (40,41,42).6 The development of 
this technique will permit the controlled investigation of 
radiation phenomena not accessible at the present time by any 
but microsecond duration electric heating techniques (20,52). 

Many current investigations of chemical and physical 
phenomena are based on studies of the emission or absorption 
tharacteristics of high temperature gases (8,15,30,44). In 
thock tube studies of kinetic phenomena in air, Camac, 


1. Lees (private communication). 
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Camm, Keck and Petty (6,7) have determined the concentra- 
tion of O. in the ground vibrational state by measuring the 
gas opacity to Schumann-Runge radiation. Such studies 
yield information on both vibration and chemical relaxation 
ratesinair. Roth and Gloersen (39) have studied the visible 
continuum emission from xenon in the temperature range 
6000 to 11,000 K and have determined the activation energy 
for this emission corresponding to the energy of the metastable 
state of the xenon atom (39). Windsor, Davidson and Taylor 
(47) have studied the energy exchange between translation 
and vibration and rotation in CO through study of the infrared 
emission. 

The fundamental knowledge concerning internal structure 
of gas molecules, collision phenomena, reaction rates and 
relaxation processes formerly supplied primarily by studies 
of the radiation from electrical discharges and combustion 
processes,® is now being greatly extended through the use of 
shock tube facilities for the generation of high gas tempera- 
tures under controlled density and temperature conditions. 
Information derived from these studies is making a signifi- 
cant contribution to the solution of the problem of radiation 
from high temperature gases in thermal equilibrium. For 
an excellent review of the use of the shock tube for the in- 
vestigation of physical and chemical phenomena the reader 
is referred to the paper by Penner, Harshbarger and Vali (34). 
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Ignition of Combustible Mixtures by Hot Gases 


Bureau of Mines, U. S. Department of the Interior, Pittsburgh, Pa. 


_ Fuel-air mixtures are ignited by a jet of hot gas which is 
heated in a ceramic furnace. The igniting gas is con- 
tinuously injected into the cold explosive mixture. Lu- 
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minous reactions can be seen to occur in the center of the 
hot jet if its initial temperature approaches a critical 
temperature. On further raising the temperature of the 
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jet, the luminosity increases and a laminar luminous pen- 
cil extends up to 30 cm above the furnace exit. A flame 
can then be seen to originate from the top of this column 
and propagate into the cold mixture outside. Two types 
of experiments are performed: The ignition of cold fuel- 
air mixtures by hot neutral gases, such as nitrogen or car- 
bon dioxide; and the ignition of cold fuel by hot air. In 
the second case only a diffusion flame can be ignited which, 
under favorable conditions, floats on top of the hot air jet. 
These ‘Shot gas ignition temperatures” bear little relation 
to the “‘spontaneous ignition temperatures.”’ 


Introduction 


HE ignition of combustible mixtures by hot burned 
gases from flames is of importance in many appliances, 
such as ramjet burners or rocket combustion cnambers. In 
the ignition of methane-nitric oxide mixtures by pilot flames, it 
was found (1) that the size and temperature of the pilot 
fame determined ignition. For small pilot flames (volume 
flow of the order of a few cc/sec), the parameter responsible 
for ignition was the energy flow of the flame (cal/sec), whereas 
for large flames, ignition became independent of flame size 
and was determined by a minimum temperature of the flame 
(the requirement being expressed in cal/cc). Intermediate 
pilot sizes had to be characterized by both energy and tem- 
perature statements for ignition of the surrounding mixture. 
It was first proposed to extend this investigation by burning 
pilot flames in hydrocarbon-air mixtures. However, both 
extremely small and larger but highly dilute, i.e., cooler, pilot 
fames invariably led to ignition, and minimum conditions for 
ignition could not be investigated in an apparatus allowing 
only relatively small flow rates. It was necessary therefore to 
heat gases in a furnace, and then lead them into explosive 
mixtures. It is difficult to heat small gas flows in a furnace 
because of inevitable heat losses at the furnace exit. The 
present investigation deals, therefore, with the minimum 
condition of fuel-air ignition by relatively large flows of hot 
gas. The influence of velocity of the hot gas on ignition could 
not be investigated for reasons explained later. This in- 
vestigation was limited to ignition by a steady stream of hot 
gas. Ignition by sudden bursts of hot gas is outside the scope 
of this paper. ' 


Experimental Arrangement 


The source of ignition was a steady flow of hot nitrogen, 
carbon dioxide, air, etc. The gas entered the furnace from 
below. After being heated, it flowed into the explosive mix- 
ture or pure fuel, which in turn flowed slowly through the outer 
vessel (Fig. 1) to avoid dilution by the jet of hot gas. This 
jet of hot gas was surprisingly coherent; for instance, a sheet 
of paper held 10 cm above the jet was punctured initially with 
avery small hole followed by ignition of the paper. The 
furnace consisted of an inner ceramic tube through which the 
gas flowed. This tube was surrounded by a spiral of platinum 
wire embedded in ceramic paste. A larger tube enclosed the 
inner tube and was sintered to the inner tube with the ceramic 
paste. Additional ceramic tubes around both these tubes 
served as insulation (Fig. 2). Gases could be heated to a 
maximum of 1400 C before the furnace was damaged due to 
melting of the ceramic and platinum wire in localized spots 
within the furnace. 

The temperature of the gas jet was measured directly above 
the furnace outlet with a platinum-platinum rhodium ther- 
mocouple before and after each ignition experiment. The 
thermocouple wire was 0.005 in. thick; a small bead formed 
the junction (Fig. 2). Furnaces 1 and 3, which had an 
inner tube of 4 mm internal diameter, were suitable for heat- 
ing flows of around 30 to 50 cc/sec. Flows above 100 cc/ 
see showed temperature profiles with layers close to the wall 


* Numbers in parentheses indicate References at end of paper. 
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Fig. 2 The construction of the furnace 


Platinum wire 


hotter than those at the center. Flows below 20 cc/sec 
were unsuitable; the layers close to the wall were much 
cooler than those at the center, due to the cooling of the last 
5 mm of the inner tube which had no heating coil, since it 
protruded through the water jacket. Furnace 2 was suitable 
for slightly larger flows, since the inner tube was larger (7 mm 
diam). 

To keep results comparable, flows of 35 cc/sec were always 
used (unless stated otherwise) for furnaces 1 and 3, and a 
flow of 60 cc/sec for furnace 2. Helium and hydrogen jets 
were difficult to work with, because heat dissipated very fast 
at the outlet of the furnace, and this could be overcome only 
by larger flows. In the initial stage of an experiment, the 
igniting gas flowed through the furnace, and steady furnace 
conditions had to be attained. At this stage, an inert gas, 
such as nitrogen or carbon dioxide, was used for the outer 
gas stream instead of the explosive mixture to be ignited. 
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a ss observed. The inert atmosphere was introduced again, 


_ The temperature at the center of the hot jet was measured 
aml above the inner tube. The thermocouple was then 
removed, and the explosive mixture was substituted for the 
Ignition or nonignition 


and the temperature checked with the thermocouple. The 
_ outer explosive mixture to be ignited by the hot jet flowed at a 
steady rate, usually 100 to 300 cc/sec. No difference in 
minimum ignition temperatures could be observed between 
these flow limits. The outer explosive mixture was kept as 
close as possible to room temperature. However, as the cool- 
ing water in the furnace jacket usually was at a slightly higher 
temperature, the gas may have been closer to 40 C than to 
room temperature. 

Ignition temperatures with a hot nichrome wire were meas- 
ured for comparison with the hot gas ignition temperatures. 
A 40/1000-in. wire 9 cm long was suspended horizontally in a 
tube 9 cm in diameter and 50 cm in length. The explosive 
mixture flowed past the wire at a rate of about 200 cc/sec, 
corresponding to a flow velocity around 3 cm/sec. Our pur- 
pose was not to make a detailed study of ignition by hot wires, 
but to obtain ignition temperatures for a variety of fuels 
under comparable conditions. The temperature of the wire 
was measured with a pyrometer and reported directly as 
measured, i.e., the black body temperature of the wire at 
6500 AE. No correction for emissivity was applied, because 
the wire acquired an oxide coating following some preliminary 
experiments, making the emissivity of the surface doubtful. 
The wire was heated by slowly increasing the electrical cur- 
rent, following the temperature rise of the wire with the py- 
rometer. The ignition temperature was defined as the tem- 
perature immediately prior to ignition, detected by the 
sound. All gases used were taken from cylinders (Matheson 
Co.). Carbon monoxide was of commercial grade. 


The gas ignition temperatures were reproducible within 
+10 C, for example when checking the value for ethane fol- 
lowing measurements with propane or a similar hydrocarbon. 
The ignition temperatures, however, were not the same for 
different furnaces and depended on the history of each furnace. 
Because of the heat losses of the hot gas jet in the exit of the 
furnace, the temperature across the jet was not a constant, 
and the resulting profile depended on the conductivity of the 
sintered material on the furnace top. After prolonged use, 
the furnace usually developed cracks in this sintered material, 
leading to smaller heat losses and an improved temperature 
profile, which in turn led to a lowering in ignition tempera- 
tures. After aging, all furnaces built during this investigation 
gave ignition temperatures agreeing to within +20C. Anew 
furnace, however, initially gave values up to 150 C too high. 
Intentional aging would have been beneficial, but it was 
difficult to predict at which temperatures aging or total 
destruction of the furnace would have occurred. 

Accordingly results are reported with the following limita- 
tions: Furnaces 1 and 2 were aged before the measurements 
recorded in Tables 2, 3, and 4 were made, and ignition temper- 
atures were reproducible throughout. Furnace 3 was new 
when used for the recorded data, but aging occurred during 
the experiments (see remarks in Tables 2 and 4). 

The temperatures recorded were measured directly with the 
thermocouple and potentiometer. It is believed that heat 
losses of the thermocouple due to conductivity are very small 
as a great length of wire was immersed in the hot gas jet. 
No radiation corrections were applied, since they are un- 
certain and may falsify the results. Asa guide, the corrections 
in Table 1 are believed to be approximately valid for an air 
or nitrogen jet of 35 cc/sec, from a 4 mm diameter tube. 

Ignition temperatures with the hot wire were reproducible 
to within +10 C or better. Difficulties were encountered 
only with very rich ethylene mixtures, rich acetylene mixtures 


Accuracy of Results 


a b c d 
Fig. 3 (a) Flame developing with a 35 cc/sec hot air jet in ethane, 
(b) Luminosity within hot argon jet in stoichiometric ethane-air 
mixture. (c) Luminosity within hot argon jet in a very rich 
ethane-air mixture. (d) Luminosity within hot argon jet in a 
very lean ethane-air mixture 


Corrections for thermocouple, deg C 
« 


Table 1 


Temperature 


as measured Correction 


700 


1100 
1300 


and all hydrogen mixtures except those near the lean limit. 
In these cases, catalytic surface reactions contributed to the 
wire heating. No results are reported for these mixtures. 


The process of hot gas ignition will be described first. Two 
types of ignition flames must be distinguished: 1 Flames 
which develop when a jet of hot air flows into cold fuel, and 2 
flames which develop when hot air or hot inert gas flows into a 
cold explosive mixture. 

1 A jet of hot air at slightly below ignition temperature 
introduced into a cold fuel, such as ethane, leaves little visible 
trace of reaction. When the air temperature is raised, 4 
small flame suddenly appears floating about 5 cm above the 
furnace tube. This flame has a bluish base topped by a 
bright green section (Fig. 3a). When the temperature of the 
jet is reduced, this flame disappears. When the temperature 
of the air jet is increased above the minimum ignition tempera- 
ture, this flame may become self-propagating and may flash 
back onto the furnace tube forming a diffusion flame with 
abundant carbon formation. Reduction of the air tempera- 
ture no longer extinguishes this flame. For the 7 mm jet, 
this flame has a greater tendency to flash back and sit on the 
furnace tube, whereas smaller jets (4 and 2 mm diam) give 
rise to floating flames which give the impression of stationary 
self-ignition flames. 
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2 An inert gas, such as nitrogen, flowing into a cold ex- 
plosive mixture gives rise to a luminous column where slow 
reactions take place at temperatures below the minimum 
ignition temperature (Fig. 3b). Increasing the temperature 
of the jet makes the column longer and brighter. Finally, the 
column branches out at some point up to 30 cm above the 
furnace tube, and a flame propagates throughout the vessel. 
Areport and the blowing off of the lid are further evidence of 
ignition. Ignition can also be approached by making the 
outside mixture too rich and slowly reducing the fuel content. 
The luminous column is quite bright, with a purple color for 
the first 10 cm of the jet, followed by a bright green region. 
As before, ignition occurs from the top of the column. Figs. 
3e and 3d show luminous columns with a jet temperature 
high enough to ignite stoichiometric mixtures, the actual mix- 
ture being either just too lean or too rich for ignition. 

Table 2 shows ignition temperatures for hot air jets flowing 
into pure cold fuel. 1n all instances, small flames are formed 
which float between 2 and 5 cm above the furnace exit. In 
some instances, this flame may subsequently flash back to 
form a full diffusion flame. 


Table 2 Ignition temperaturesin deg C for a hot air jet 
into cold fuel 
Furnace 1 Furnace 2 Furnace 3 
(aged) (aged) (new) 
sd Jet. diam Jet diam Jet diam 
4 mm 7mm 4mm 
Fuel 35 ec/sec air 60 cc/sec air 35 cc/sec air 

methane 1180 1135 1180? 
ethane 920 905 990? 
propane 980 960 1070 
n-butane 960 ae 1075 
ethylene 840 840 890° 
propylene 1040 1020 1085 
iso-butylene 1030 1110 
carbon mon- 

oxide 745 760 765 
hydrogen 655 670 665° 
acetylene 700 7007 

* Furnace aged 

> 920 C after aging 

© 845 C after aging 


Furnaces | and 2 were used at high temperatures for pre- 
liminary experiments, and therefore were aged, giving lower 
values than furnace 3, which was used directly for final meas- 
urements. The hydrogen and acetylene experiments in fur- 
nace 3 were run following an experiment (referred to in 
Table 4) in which the jet temperature was raised to 1300 C, 
with subsequent aging. The methane experiments were last 
in each set. 

In one experiment, 35 cc/sec oxygen was substituted for air 
in furnace 1. This reduced the ignition temperature of pro- 
pane from 980 to 960 C. With hydrogen and carbon mon- 
oxide, it is possible to reverse the ignition experiment and 
cause ignition in air with a jet of hot fuel. Table 3 shows the 
results. 


Table3 Ignition temperatures in deg C for air into fuel 


and fuel into air with furnace 1 
Hot jet Cold medium Ignition temp. 
air 655 
air 820 
air CO 745 
CO air 860 
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Fig. 4 Change of hot gas ignition temperature with mixture 
strength 


Table 4 gives ignition temperatures of fuel-air mixtures 
ignited by a jet of hot nitrogen. The temperatures recorded 
correspond to the mixture compositions which gave the lowest 
ignition temperatures. 


Table 4 Ignition temperatures in deg C for a jet of hot 
nitrogen into fuel-air mixtures 
Furnace 2 Furnace 3 

Fuel in fuel-air (aged) (new) 

mixture No jet, 60 cc/sec Nz jet, 35 ce/sec 
methane >>1200 Close to 1425 (aged) 
ethane 1015 1110 

(1050 after aging) 

propane 1110 1200 
n-butane 1095 1220 
ethylene 985 110 
propylene 1150 1305 
iso-butylene 1155 1295 
carbon monoxide 795 905 
hydrogen <755 not tested 
acetylene nottested not tested 


Methane-air mixtures could not be ignited. At 1425 C 
(furnace 3) a strong luminous column formed in the hot jet, 
indicating that the ignition temperature was near, but further 
raising of the temperature made the furnace break. Furnace 
2 broke above 1200 C when trying to ignite methane-air. 
Hydrogen-air mixtures gave a violent explosion when ignited 
with a jet of nitrogen at 755 C, so the accurate ignition tem- 
perature could not be measured. Acetylene was not tested 
for the same reason. 

An attempt was made to measure not only the lowest 
temperature at which a 35 cc/see hot nitrogen jet ignites the 
most favorable fuel-air mixture, but also the ignition tempera- 
tures over the whole range of mixture compositions. Fig. 4 
gives ignition temperature (furnace 3) vs. mixture composition 
curves for butane, ethylene and carbon monoxide. The 
curves are quite flat, i.e., ignition temperature does not 
change much with mixture composition, making it difficult to 
locate the minimum accurately. Hydrocarbons seems to 
ignite slightly better when the fuel is rich, with the exception 
of ethylene for which the minimum is on the lean side. With 
propylene the minimum corresponds to a stoichiometric mix- 
ture, whereas with butylene the minimum is again on the 
fuel-rich side. Carbon monoxide ignites much more readily 
when the fuel is rich. 

The ignition of fuel-air mixtures by hot air, instead of nitro- 
gen, also has been investigated. (See Fig. 5.) For pure fuel 
the values are identical to those in Table 2 (furnace 3). The 


addition of air even in rather large amounts affects the igni- 
tion temperature very little. 


The resulting flame obtained is 
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Fig. 5 Ignition of fuel-air mixtures by hot air (35 cc/sec). 
- Arrow marks stoichiometric point 


simply a suspended diffusion flame. As the mixture composi- 
tion approaches the rich limit, the ignition temperature be- 
gins to rise and continues steadily upward toward the lean 
limit. Within the explosive range, a flame moves through 
the whole vessel after ignition, instead of remaining suspended. 
Table 5 summarizes the ignition of stoichiometric fuel-air 
mixtures by a 35 ce/sec jet of hot air. 


Table 5 Ignition of a stoichiometric fuel-air mixture by 
hot air (35 cc/sec) furnace 3 (new) 


Fuel Ignition temp., deg C 
methane not tested 
propane 1160 
n-butane 1170 
iso-butylene 1215 
carbon monoxide 870 


The influence of the nature of the hot gas on ignition has 
been further investigated by substituting carbon dioxide, 
argon or helium for nitrogen or air as the igniting gas. Table 
6 gives details for ethylene-, ethane- and carbon monoxide-air 
mixtures ignited by hot nitrogen, argon, helium or carbon 
dioxide jets. 


Table 6 Ignition temperature of most favorable mixture 
by inert hot jets, deg C (35 cc/sec) furnace 3 


Carbon 

Fuel sisal Nitrogen Argon Helium dioxide 
ethane «1290 
ethylene 1110 1165 1265 1075 
carbon monoxide 900 925 1000 870 


The foregoing values are the result of separate experiments, 
and the values for nitrogen are not identical, although close, 
to those reported in Table 4. The flow of helium was ad- 
justed to 100 ec/sec for the ethane and ethylene experiments, 
and to 60 cc/sec for the carbon monoxide test, because the end 
effects of the furnace influenced the helium temperature too 
much with a flow of 35 ec/sec. The furnace tended to over- 
heat for low helium flows. Carbon dioxide seemed to be as 
effective as nitrogen for ignition. Argon required slightly 
higher temperatures, but helium did not ignite well at all. 
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The oxygen index of the explosive mixtures has little in- 
fluence on the ignition temperature. An increase of the 
oxygen index from 0.21 to 0.25 in an ethylene-air mixture 
ignited by hot nitrogen reduced the ignition temperature by 
less than 20 C. Higher oxygen indexes could not be con- 
sidered, because the explosion became too violent. 

The results of the wire ignition experiments are summarized 
in Fig. 6. The ignition temperatures are not very sensitive 
to mixture composition. With the exception of carbon mon- 
oxide, all ignition temperatures are lower for lean flmes, 
However, the slope is more pronounced for methane, propyl- 
ene, butane and butylene than for the other gases. Hydrogen 
and acetylene could only be measured satisfactorily on the 
very lean side for reasons already explained. The oxygen 
index was altered for an ethane- ‘air’ mixture. The ignition 
temperature is affected very little by this change (see Fig. 7), 


Little information is available in the literature on the igni- 
tion of explosive mixtures by hot gases. Wright and Becker 
(2) appear to have been the only investigators to attempt this 
problem. They heated air or nitrogen on quite a large scale 
with a heat exchanger and circulated fuel-air mixtures around 
the hot jet. They were able to heat the inner jet to just 
above 1000 C. This was only sufficient to ignite an outer 
acetylene-air mixture, which required an inner jet temperature 
of 752 C. Propane could be ignited only when the outside 
mixture also was heated. Thus hardly any data are avail- 
able for comparing Wright and Becker’s work and the work 
reported here. Their value for acetylene is slightly higher 
than our value for hot air into pure acetylene, but thev do 
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Fig. 6 Ignition of fuel-air mixtures by a 0.04-in. nichrome wire 
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Fig. 7 Ignition of C.H,-‘‘air’’? mixture by hot nichrome wire. 
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Table 7 
Minimum ignition energy 
as ignition Spontaneous Stoichiometric Optimum 
temp.” ignition temp. all = mixture mixture 
Fuel deg C with air with air with air 
; deg C 
methane = 537 03 0.28 X 1078 joules 
ropane 980 493 0.40 
| ethylene 840 490 0.11¢ 
propylene 1040 458 7 0.41° 
iso-butylene 1030 4300 
carbon monoxide 745 609 
| hydrogen 655 572 0.02 0.02 © 
| acetylene 700 305 0.03¢ 
| “Interpolated from Figs. 185, 187 in (4) 
’From Table 2 | 
From Calcote in (5) (Flange electrode) 


not say whether this value is for hot nitrogen or hot air. 
Moreover, their mass flows were very large compared to ours. 
However, Wright and Becker were the first to note that a 
luminous reaction precedes the spreading of the flame. 
They call this reaction an “‘initial flame.” We doubt whether, 
at least in our case, this is a flame, as it does not have such 
flame characteristics as self-propagation. Accordingly, we 
prefer to call it a luminous column. 

It is tempting to compare ignition temperatures as meas- 
ured here with other flame characteristics, such as spontane- 
ous ignition temperatures (as measured in furnaces), minimum 
spark ignition energies, etc. The values found in the litera- 
ture on spontaneous ignition temperatures vary widely; 
the values recommended by Scott et al (3) are recorded in the 
second column of Table 7. They are compared to the lowest 
values obtained for ignition by hot gas, i.e., when hot air is 
injected into pure fuel. It is clear that little or no correlation 
exists. For hydrocarbons the hot gas ignition temperature is 
very roughly double the spontaneous ignition temperatures, 
when expressed in deg C, whereas that for hydrogen or carbon 
monoxide is only slightly higher. Even within a series of 
hydrocarbons, there is no correlation. Table 4 shows that 
this also holds for ignition with hot nitrogen. 

Table 7 also gives minimum ignition energies as reported in 
(4). If there is a correlation between hot gas ignition and 
minimum spark ignition energy, one would expect the rela- 
tion to be with the energy for the most favorable mixture. 
However, paraffin hydrocarbons show very little change. A 
certain correlation does appear to exist for the values taken at 
stoichiometric ratio, especially as hydrogen and acetylene 
have very low minimum ignition energies. If carbon mon- 
oxide also had low ignition energies this would be an important 
proof. However, this is very unlikely, and no measurements 
could be found in the literature. Any anticipated correla- 
tion breaks down, if we consider that an increase in oxygen 
index has very little influence on the hot gas ignition tempera- 
ture, whereas it affects the minimum ignition energy greatly 
(see (4)). 

It is because of this lack of correlation that the hot wire 
experiments were performed. A comparison between Fig. 6 
and either Table 2 or Table 4 shows that the correlation is 
nearly perfect. This holds also for the values obtained with a 
different oxygen index, as in both instances little change is in- 
dicated. This close correlation suggests that an important 
characteristic of flame ignition has been found. 

Correlation with lean limit temperatures leads to some 
agreement as hydrogen, carbon monoxide and acetylene have a 
low lean limit temperature compared with paraffin hydro- 
carbons. (For values see aaa (6).) However, methane 
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does not fit inte the order, since its limit temperature is much 
the same as that of other hydrocarbons. It is also known that 
the oxygen index has little effect on the lean limit. 

The ignition process itself will now be considered in more 
detail. It is important to note that when hot nitrogen is in- 
jected into a fuel-air mixture, a luminous reaction occurs in the 
core of the hot jet. The hot jet gradually loses heat to the 
surrounding gas, but at the same time mass transfer takes 
place, i.e., nitrogen flows radially from the jet and is replaced 
by an explosive mixture. This explosive mixture will react if 
the residual temperature of the jetishigh enough. Ifthe heat 
produced along the jet is less than that lost by conduction, the 
reaction dies out, and there is no ignition. If the initial jet 
temperature is increased, the luminous column becomes longer, 
and a propagating flame originates at the center of the jet. 
Theoretical calculations for a two-dimensional model (hot gas 
in contact with cool explosive mixture) have predicted this 
behavior (Marble (7), Cheng (8), Spalding (9)). Wright and 
Becker were able to observe it in their high-speed ignition 
arrangement. It has now been shown to be applicable in a 
wide variety of experiments. 

Therefore, ignition is dependent on the reaction rate at a 
given temperature level in the core of the jet. The ignition 
temperatures in Table 4 measure the relative reaction rates of 
the fuel-air mixtures. However, the reaction rate is not the 
only factor influencing ignition. If it were, a jet of hot air into 
cold hydrogen would have the same ignition temperature as 
hot hydrogen into cold air. In fact, the temperatures differ 
by 165 C. Mass transfer must be considered. More air 
must diffuse into hydrogen to form a reaction mixture than 
hydrogen into air, for reasons of stoichiometry. More in- 
formation as to the importance of diffusivity and the specific 
heat of the jet is available from the experiments with hot 
argon, helium and carbon dioxide. Carbon dioxide is not 
essentially better than nitrogen, and the heat content or 
energy of the jet does not seem to be of great importance. 
The ignition temperature of argon is only slightly higher than 
that of nitrogen. However, the high diffusivity of helium is 
clearly responsible for the large increase in ignition tempera- 
ture as compared toargon. The higher jet velocity which had 
to be used with the helium jet is not responsible for the in- 
crease in ignition temperature, since the influence of jet 
velocity on ignition temperature is small for laminar jets. 
This subject will be considered in a separate paper. 

That the oxygen index is of little importance is also under- 
standable, since temperature controls the reaction rate to a 
greater degree than does concentration, and dilution with 
nitrogen decreases the ignition temperature only slightly. 

It has been known for some time that gas ignition by hot 
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wires occurs in the wake of the wire, at least as long as the 
wire does not give rise to catalytic effects. This ignition also 
will be determined by whether the reaction rate is greater 
than the heat loss by conductivity. The close similarity 
between hot wire ignition and hot gas ignition becomes under- 
standable. More difficult is the question of why hot wire 
ignition is so insensitive to mixture concentration, whereas 
hot gas ignition shows a minimum. A definitive answer can- 
not be given at the moment. It can only be pointed out that, 
whereas hot gas ignition depends on a reaction rate in a given 
volume, hot wire ignition can expand the useful volume of 
reacting gases along the wire axis. 

Although diffusivity-as well as reaction rate determines igni- 
tion for hot gas and hot wire ignition, it is rather surprising 
that so little similarity exists between those modes of ignition 
and spontaneous gas ignition. Hot gas and hot wire ignition 
require higher temperatures than does spontaneous ignition. 
As different fuels have very different activation energies, a 
rise in temperatures entirely reverses the reaction rates, i.e., 
carbon monoxide seems to have a higher reaction rate than 
propylene at 800 C, whereas at lower temperatures the reverse 
is true. 

More information is needed to establish the influence of jet 
diameter and jet velocity. In this investigation, the diameter 
was varied from 2 to 7 mm, and ignition temperature showed 
little change. The results with the 2 mm diameter were not 
reported here, as they were less accurate due to the difficulty 
of measuring temperatures in such a limited space. Much 
larger diameters must be used in the future to investigate this 
effect. The velocity of the hot jet could not be altered sig- 
nificantly in our experiments because of the change of the 


Investigation of the diam. 
eter and velocity effect in the future is hoped for. 


temperature profile of the hot jet. 


1 wish to thank Mr. A. E. 
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the ignition experiments. 
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A digital computer program for engine transient study 
has been developed to obtain an evaluation of pressure 
surges in the propellant system during startup and shut- 
down. This technique is intended to supplement and 
eventually replace analog computer start and shutdown 
transient analysis which does not provide realistic simula- 
tion of propellant lines and combustion dead times. The 
mathematical model used for propellant line representa- 
tion was based on the water-hammer theory as origin- 
ally developed by Allievi, i.e., hydraulic oscillations 
appear as pressure surges caused by the conversion of 
kinetic to potential energy, where the kinetic energy 
changes are a result of system impedance changes taking 
place during engine start or shutdown. Line equations 
and component transfer functions have been programmed 
in a general manner, allowing routine parameter changes 
to be executed without the need of major reprogramming. 
The analysis determined the values of instantaneous pro- 
pellant flow and transient pressures at points of interest 
in the system for successive time intervals. 
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Nomenclature 
A = area, ft? 
b = pipe wall thickness, ft 
c = velocity of sound, ft/sec 
c* = characteristic velocity, ft/sec 
Co = specific heat at constant volume, Btu/lb/F 
Cp = specific heat at constant pressure, Btu/Ib/F 
= constant 
D = pipe diameter, ft ae 
E = pipe modulus of elasticity, lb/ft? af @ 
e = water-hammer slope, Ib-sec/ft* 
e’ = modified water-hammer slope, sec/ft? 


= function denoting a pressure wave 

= function denoting a pressure wave 

= acceleration of gravity, ft/sec? 

= film coefficient of heat transfer, Btu/sec F ft? 

heat of vaporization, Btu/Ib 

= moment of inertia, ft-lb-sec? 

= constants or bulk modulus, (min/rev)?, lb/ft? 
thermal conductivity, Btu/sec F ft 


length, ft 


molecular weight, lb/mol 
= interaction index 
= Nusselt number 
= pressure, lb/ft? 
= Prandtl number 
= volume flow, ft?/sec 
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= heat flow, Btu/sec 
gas constant, ft/R ; 
resistance coefficient, sec?/ft?-] 
Reynolds number 

gear ratio 

temperature, R 

turbine power transmission losses 
time, sec 

volume, ft* 

= fluid velocity, ft/sec 
= weight, lb 

flow rate, lb/sec 
= LOX quality 

= density, lb/ft® 
= difference 

= efficiency 

= pump or turbine speed, rpm 
= torque, ft/lb 

= kinematic viscosity, ft?/sec 
= valve or injector admittance, lb-ft?/sec? 
= combustion time lag, sec 

= pump constant, sec®-ft/lb 

= angular velocity, radians/sec a > 
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Subscripts 
= burned propellants 

= thrust chamber a 
= pump discharge 
= flow 
= fuel or filled condition 
= gas 

heating 

= inlet 

= injector 
= line 

= liquid 

= metal 

= efficiency 

= oxidizer 

= pump 

= pump suction 

= turbine 
= time or thrust chamber throat 
= turbopump 

= valve 
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Introduction 


METHOD was developed to analyze transients in liquid 

rocket engines by means of a high-speed digital compu- 
ter. The basic mathematical model was developed and its 
physical counterpart is shown schematically. An example 
illustrates the use of the method on a pressure fed system, and 
the results are compared with test data. 


Object of Transient Analysis 


The reproducibility of desired starting and shutdown tran- 
sients is one of the most important factors in designing a rocket 
engine with a specified performance. Almost all malfunctions 
and performance shortcomings result from an unexpected de- 
velopment during the starting transient. Conversely though, 
not all unusual transients will result in component failures. 
Therefore, an understanding of the mechanism of transient 
phenomena is required to predict system tolerances and their 
effect on the prescribed performance. The analysis relates 
pressure and flow oscillations in the propellant lines to tran- 
sient phenomena during starting and shutdown. 


Mathematical Model 
The physical system analyzed is shown in Fig. 1. The fol- 
lowing is the breakdown required for mathematical represen- 
tation of the main components: 
The mathematical model for the propellant lines was based 
on water-hammer theory as originally developed by Allievi, 
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¢ Fig. 1 Basic rocket engine schematic used for digital computer 
representation 


= 
| 
ie., hydraulic oscillations appear as static pressure surges 
caused by the conversion of kinetic to potential energy. 
From (1)* the conjugate water-hammer equations are given 
by 


= e(vz, — 


—el(vz, 


where e = cy/g. 
Equation [1] assumes frictionless flow. To account for 
frictional pressure drop, two assumptions were made: 
1 Frictional drop is proportional to the square of the line 
flow. 
2. Friction may be lumped at the end of each pipe section. 
_ The form of Equation [1] as used for digital programming will 
_ be shown later. 
The characteristics of the mixed flow centrifugal propel- 
lant pumps are given by 


Da — Ps = — | vp | 


The values of Kg and y, are determined by solving the equa- 
tion at two values of pa — p; and v, for the same 6, where pa — 
ps and v, are taken from the head capacity curves of the 
pump. Other type pumps can be represented by modifying 
the right hand side of Equation [2]. The pump efficiency is 
represented by the equation 


Q 4 - 


Pp Pp 


The values of K,, and K,,’ are determined in the same manner 
as Ky and y, where Q, and 7, are taken from the pump ef- 
ficiency curves. The pump torque for a single pump is given 
by 


where h = 27/60, conversion factor from rpm to radians/sec. 
In the above equations, Q, and v, are related by Q,= Apu, 

where A, is the pump inlet area. 
Substituting [2 and 3] into [4] 


Ko 6,? — 


Equations [2, 3, 4 and 5] apply for either the fuel or the 


Mp = 


3 Numbers in parentheses indicate References at end of paper. 
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oxidizer pump. Total torque required by the two pumps is 
given by yu, where 


where 7’, = transmission losses, 0/0, = r, and 6/8, = r, rep- 
resent the gear ratios of the fuel and oxidizer pumps, 
respectively. 

The accelerating torque 


For a short ‘ies interval, the torque can be written approxi- 
mately as 


It, 
—At) + ur) = — [w, — wit—at)] = 


? — 6:).. 


or h = 
k D 0.8 
where (Ci = 0.023 — (7) [16] 
D \pvAr 


Substituting in [12] and rearranging 


The heat fae. nt into the liquid” must come from heat 
capacity of the metal, hence 


Equation [10] holds while the pump is operating as a pump. 
When the flow reverses, y; is replaced by y which is esti- 
mated to be 3y,. The range of transient behavior considered 
includes positive shaft speed, positive head and positive and 
negative fluid flow. 

The pump heads taken from the head capacity curves in- 
clude the difference between the velocity heads at the dis- 
charge and suction points. Since the desired pressures are 
static pressures, an allowance must be made for the velocity 
head by subtracting the differences between the rated system 
pump head and the calculated rated pump head. 

The valve equation used for either propellant line is given by 


wy? = Av? (Pit 


\,, is the valve admittance factor and is usually available as 
a function of valve opening. 

The condition of the oxidizer (liquid oxygen) in the tank and 
feed lines is at, or slightly above, the saturation temperature 
at atmospheric pressure (normal boiling point). When the 
LOX valve is opened, liquid flows through the valve into the 
injector manifold. Since the injector is at ambient tempera- 
ture, the LOX will absorb heat from the injector and vaporize. 
The quality of the oxygen and the weight flow rate of the 
liquid into the combustion chamber had to be determined as 
a function of time from initial valve opening. 

Since the liquid is at its saturation temperature, all the heat 
will be used to vaporize a portion of the liquid. The heat 
transferred will be primarily by convection. 


q/An = h(Tm 
ae The forced convection heat transfer equation is 


Nu = 0.023(Re)-8(Pr)®3. . 


hD 
— = 0.023 
k v 


or 


where w,’ = 


Substituting [6] into Equation [7] iA! 4 
5 + Mp, + (Mery + = ? — 4)... [9] 
6? 62 2 62 d 
1 Ke, ry Vi Ke, po| Ke, = V1, Ko, — 
2h f Q t Qos Q 


Continuity requires that 
[20] 


where / and g denote liquid and gaseous conditions. 
The quality of LOX, X, is defined as the fraction by mass 
of liquid in the mixture 


In Equation [17] q is also equal to h,w,., where h; is the heat of 
vaporization of oxygen 


> 


— Wort) = hy(to — = (Tn — To 
Aun 
hy 


and substituting from Equation [19] 


Ay 


From the mass balance, the weight of oxygen flowing into the 
injector must equal the weight flowing out plus the rate of 
change of gaseous and liquid oxygen contained in the injector, 
i.e. 


dW, 
dt dt 


LOX flow rate through injector. 
If the perfect gas law is assumed, P,;V = (W,./M)RT, the 


Wo = Wo’ + — + — [24] 


JET PROPULSION 


tot 


wh 


ti 
tir 
fr¢ 


di 


F 


: 
ant 
V 
dt oF 
. lo 
ma 
flo 
OF thls CQuation yleias 
se 
= t 
Tiny = To + (Tm — Tode 
Ap 
‘ 
22) 
_ 
Dy) 0.3 
f D pvAr 


[20] 


nass 


[21] 
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[22] 


total change isexpressed it 
RW, \ dV. 1 dW, 
1 g 25 
also from [20 and 21 ee 
4 
26 
and 
V, = volume of injector manifold (V,;) — volume of LOX in 
injector manifold = V,; — [27] 


Equations [23 through 27], after being expressed in finite 
difference form are sufficient to express the injector pressure 
P,, in terms of the valve flow t»,. 

rom the assumed incompressibility of the fuel (in the 
manifold), after the injector is full, the flow out must equal the 
flow in, if no leakage flow is assumed during filling. 


O t<t 
= ; [28] 
Wy ty 
where t, is the time at which the injector is filled, also 
tf 
f w = capacity of fuel manifold by weight. .... [29] 
o 


Equations [28 and 29] allow the solution for t,. Air and fuel 
vapor compression during the filling interval were neglected, 
since experimental evidence indicates this factor to be of 
secondary importance. 

After the manifolds are filled, the injector equation is 


The injector admittance \,; is constant for all flows. 

Propellants injected into the combustion chamber after 
time (t — 7) will not have burned at time ¢ due to a combus- 
tion time lag 7. The total amount of propellants burned 
from t’ = 0 tot’ = tis given by 


t 
f, = 
differentiating 


Wy = at [32 
From Crocco and Cheng (2) 
dr n 
Using Equation [32] 7 


n 
Woy = (Pec) Penn) | 


At Ve dp 


= Poa + RT. [34] 
but 
Wiltt—r) = (Writ 
+ tot) E De = Pann) | = 
c 
At Ve dpe 


where W1) or W(;—7) indicates time at which variable is 
evaluated. 

The aforementioned components of the mathematical 
model are typical and identical in form to the models used for 
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the gas generator loop. For this loop, however, a number of 
additional equations were necessary, describing the quasi- 
steady condition of the hot gases at the turbine outlet. 


Programming Details and Computing Logic 


The programming of the mathematical model was done in 
a general form. Not only can system parameters be changed 
at will, but components of the system may be replaced with 
other components of different characteristics. To some ex- 
tent, rocket engine models may be synthesized with available 
component models. 


Techniques 


All system components, such as valves, injectors and 
pumps, were numbered and given a type annotation. The 
equations describing the components and the equations con- 
necting the component to the upstream and downstream com- 
ponent were then programmed. 

Since the connecting equations are all line equations, any 
line component can be linked to another line component when 
the type and the next upstream component to which it is to 
be connected are specified. This scheme allows line con- 
figurations to be changed at will. 

An important digital feature of the compressible fluid flow 
equation is the fact that the equations use past time values for 
the extreme ends of the upstream and downstream lines. This 
allows each component to be computed without the necessity 
of knowing present time values at the other components. 

The thrust chamber and the gas generator combustion 
chambers are particularly amenable to this type of solution 
due to the combustion time delay. Present time pressures 
are thus computed from past time flows and mixture ratio. 

Solution of the turbopump loop required the most involved 
handling. Here the energy input to the system is transferred 
from the gas generator to the propellant lines, and involves 
simultaneous solution of the turbopump system and the 
propellant suction and discharge lines. In addition, the avail- 
able energy absorbed by the turbine is dependent cn the ex- 
haust conditions. In this section of the program, there are 
two separate iterations within an iteration. 


Modification of Line Equations 


The line Equation [1] is modified by the addition of two 
friction terms 


Pre — Pry, = (thz, tor, + tor, | 


where 


g Ay gA, 


is a constant converting velocity into mass flows, where R,’ 
is the resistance coefficient at section x, and R,,’ is the resist- 
ance coefficient at section 2x}. 

The use of the absolute mass flow term is to keep the fric- 
tion acting in the proper direction whenever the flow direc- 
tion reverses. 


Computing Time Interval 


The computing time interval and the line transmission time 
are determined as follows: 

Pressure disturbances travel through a medium with sonic 
velocity, and the transmission time of a line equals the length 
of the line divided by the speed of sound in the medium. 

In Equation [36] the transmission time determines time (,, 
since t — ¢, equals transmission time. 

When the transmission times of all the lines are deter- 
mined, the shortest time is the computing interval. For the 
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rest of the lines, the transmission time is divided by the com- 
puting interval, and the closest multiple integer is used as the 
number of computing intervals between end “ae of the 


line. 


Computer Simulation of a Pressure Fed System 

A rocket engine system, comprised of the elements shown 
in Fig. 2, was mathematically represented and programmed 
into the digital computer. The following equations repre- 


Computer Logie Line equations, similar in type to Equation [36]. 
The diagrammatic representation of a system link is shown 2 Valve equation, similar - type to Equation [11]. re 
below. 3 Manifold equations similar in type to [12 to 29]. a 
Line parameters At, e Line parameters At, e’ 
j 
Constants 
and equations 
Component |—— Component € 
PBis to h 
time ¢ — At Component V 7 je time ¢ — At 
The equations representing the computation logic for this 
schematic are 
shi 
Peis) = —e' — Wao(t—a1)] + Pao(t—ay — | | — | tao(t—at) | 
nt 2X2 X || PBo(t) sta 
Matrix of cal 
transfer function | WBo(2) || : 
= clo 
Ppo(t) = = + pyi(t—an) + | | + | | [39] sur 
4 Injector equations similar in type to Equation [30]. dy 
213 THRUST CHAMB- MW 5 Combustion chamber equations similar in type to Equa- flor 
tion [35]. 
Valve opening times and other pertinent initial data were 
found from an actual firing. 
THRUST CHAMBER 
é rat Start Transient 
a The functions plotted for comparison are chamber pressure 
(Po), fuel thrust chamber valve inlet pressure (P yo), fuel tank 
pressure (Py), oxidizer thrust chamber valve inlet pressure 
(P,,) and oxidizer tank pressure (P,;). The following points 
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of interest are highlighted in Figs. 3 and 4: 


al 


4 
| 
OPEN — 
cuosep- THRUST CHAMBER VALVES 
to 
9 4 
w 
NOTES 
? |. Dotted curves ore 
_4 from test dete 
---T~ 2. Solid ves or 
z / 
@ 
z 
2 os ay 
anal 
2 1 rock 
Re 
17] 
NON-DIMENSIONAL TIME 
Cons 
2 
Fig. 4 of start transient pressures 


JET PROPULSION 


- 
| 0) 
| m 
« 
~ 
5 
ee 
43 
piss 
~ 
EC 
ql 


1 In the actual engine, combustion begins at this point 
oxidizer and fuel begin leaking through the injector. The 
manifolds are not fully pressurized yet. 

2 The manifolds are filled at this point and are fully pres- 
surized, allowing appreciable flow to enter the chamber, re- 
sulting in a sharp rise in chamber pressure. The computer 
results do not show the slow rise in chamber pressure between 1 
and 2 because no injector fuel flow exists between these points. 

3 Initial fuel flow through the fuel thrust chamber valve, 
ndicated by a drop in valve inlet pressure. 

4 Initial oxidizer flow through the oxidizer thrust chamber 
valve. 

The computer results follow the test data in general ap- 
pearance and within about 10 per cent in magnitude. The 
nain differences are attributed to approximations of valve 
characteristics and line configuration. The effect of bends 
on pressure waves is not included in the mathematical model, 
and small changes in pipe cross-sectional area and wall thick- 

1ess are not taken into account. 


Shutdown Transient 


The points of interest for the shutdown comparison are 
shown in Fig. 5: 

1 Test data drop in the oxidizer valve inlet pressure is ap- 
varently due to the oxidizer valve lagging the fuel valve in 
starting to close. The result is an increase in oxidizer flow 
caused by the drop in chamber pressure and the increased 
flow reducing the valve inlet pressure. 

2 The valve inlet pressures increase during the valve 
closure due to the conversion of velocity head to static pres- 
sure and the resulting pressure-wave action. 

3 Indication of combustion instability due to feed-line 
dynamics. The computer indicated a loss of chamber fuel 
flow at this point. 


Summary 


Approximations in the mathematical model are those due 
to lumped resistances and the assumption of straight piping. 
The former does not materially affect the accuracy of the 
analysis. The latter, however, could be important ina system 
having numerous bends. Bends act as reflective elements, 
passing some of the pressure waves and reflecting the others. 
Thus, a length cf tubing with a bend at some intermediate 
point will show pressure oscillations of a line having a length 
up to the bend, modulated by a lower frequency oscillation 
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Fig. 5 Comparison of shutdown transient pressures 


due to the entire pipe length. This has been observed in a 
number of tests. 

In order to obtain conservative results with respect to pres- 
sure surge amplitudes, pumps were considered to pass pres- 
sure waves without attenuation. A contemplated improve- 
ment of the present digital program will incorporate a vari- 
able attenuation factor to be determined from test records. 

The present digital program requires approximately 30 
min of IBM 704 time for each engine starting transient solu- 
tion; a basic computing interval of a little over 1 millisee was 
used. A simplified program is under development which will 
neglect the transients in the pump discharge line; however, 
the transients in the longer suction lines will still be retained. 
This program is expected to require less than 10 min for solu- 
tion and will allow rapid evaluation of transients for numerous 
conditions, | 
2 
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Precision Measurement of Supersonic Rocket Sled 
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This continuation of a previous paper (Reference 1) 
analyzes the possibility of measuring the velocity of a 
rocket sled to an accuracy of 0.1 ft/sec (root mean square) 


Received Jan. 4, 1958. 

1 This paper represents work performed by the authors in their 
respective capacities as Member of the Technical Staff and 
Consultant of the Ramo-Wooldridge Corp., Los Angeles, Calif. 

2 Assistant Professor of Aeronautical Engineering. 
’ Professor of Instrumentation. 


DECEMBER 1958 


Velocity—Part Il 


over a 100 cps bandwidth. It is shown that a track coil 
system alone cannot achieve this goal using reasonable 
techniques. However, the addition of a sled-borne ac- 
celerometer of moderate accuracy suffices to insure attain- 
ment of the quoted velocity accuracy, if accelerometer and 
track coil data are combined in an optimum manner. 
Finally, an error analysis is performed on a practical com- 
putational procedure for combining accelerometer and 
track coil data. It is concluded that the desired accuracies 
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can be achieved with present day components using the 
SNORT track. 


Introduction 


ART I of this paper (1)4 was devoted to an expository 

discussion on the measurement of sled velocity to an 
accuracy of 0.1 ft/sec rms over a 100 cps bandwith. This 
part of the paper contains the analytical justification for the © 


results stated in Part I. a 
In section 2, we analyze the possibility of attaining the de- _ 


sired measurement of velocity with the existent track coil 
system. It is shown that any practicable system utilizing 
only discrete position data cannot achieve satisfactory re- 
sults, even when markers are spaced 2 ft apart with an ac- 
curacy of 0.024 in. rms, and time measurement is accurate to 
0.1 microsec. 

The third section is devoted to a study of velocity measure- 
ment through use of an optimum combination of position and 
acceleration data in a distortionless filter scheme. It is 
shown that velocity accuracy and bandwidth specifications 
can be met without imposing too stringent requirements on 
the measurement of position or acceleration. In fact, it is 
demonstrated that this approach is feasible within the present 
state of the art. 

A practical computational procedure for combining ac- 
celerometer and track coil data is presented in section 4. 
The equations which must be mechanized are given. A com- 
plete error analysis is performed, indicating that an accuracy 
of 0.1 ft/sec rms over a 100 cps bandwidth can be achieved 
by practical components and computation techniques. 


2 Differentiation of Track Position Data 


If track markers are used, the data obtained will consist of 
the time intervals required for the sled to pass from one 
marker to the next. Solely for simplicity in analysis, how- 
ever, we shall treat position as a function of time. We shall 
assume that the sled position is sampled at equal intervals 
of time. The velocity is then calculated by numerical dif- 
ferentiation of discrete data. (For a discussion of discrete 
differentiation, see (2).) 

The problem at hand may be stated as follows: Using a 
system of physical markers, is it possible to achieve a velocity 
measurement accuracy of 0.1 ft/sec with a bandwidth of 100 
eps? It should certainly be possible if the spacing between 
markers is sufficiently accurate and if the distance between 
them is extremely small. However, we shall show that the 
accuracy and bandwidth requirements cannot both be met 
by any system of markers now available or conceivable from 
the standpoint of time, manpower and the present state of 
engineering achievement. 

The block diagram of Fig. 1 shows how the velocity-meas- 
urement error arises. It is seen that there are two sources of 
this error. First, the discrete differentiation process can only 
approximate true (continuous) differentiation. Thus, except 
for restricted inputs (e.g., constant acceleration for certain 
numerical differentiation processes), there will be an error be- 
cause the operation performed on the position data is not a 
true differentiation. Second, the noise input indicated in Fig. 1 
also yields an error. 

We may now proceed directly from Fig. 1 to the power spec- 
trum of the velocity error. Let A(w) be the power spectrum 
of the sled acceleration. Since an integrator has a transfer 
function 1/iw, the sled position spectrum is A(w)/w (4). 
G(w) denotes the transfer function of the discrete differentia- 
tion process, and N(w) the spectrum of the position noise 


input. Since the transfer function of a true differentiator 

is iw, the power spectrum of the velocity-measurement error is 
le 

= Nw) + | iw — | 


‘ Numbers in parentheses indicate References at end of paper. 
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Differentiator 


Position 
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Ideal Time 
Differentiator Delay 
Fig. 1 Block diagram of error computational process 
from Fig. 1. 


Since A(w) and N(w) are fixed by the characteristics of the 
sled, track and marker system, the value of #(w) depends on 
the choice of G(w). Evidently, G(w) should be selected to 
minimize ®(w), subject to the limitation that only diser te 
samples of sled position are available with time interval h. 

The general subject of optimum discrete filters is sufficien' ly 
complex to merit further mention, but without a complete— 
and therefore prohibitively involved—analysis. The reacer 
is referred to (3), which treats in detail discrete filters that 
permit zero lag for certain classes of signals (polynomials) 
while minimizing the error due to noise. Because our data 
need not be available instantaneously, the zero lag require- 
ment is relaxed to optimum lag, which improves the filtering 
possibilities considerably. In other words, the present ve- 
locity of an object can be estimated more closely by a device 
which looks into past and future for time 7 than by a system 
which uses only data in the past over a time interval 27. 

For the purpose of this analysis, we shall construct a dis- 
crete differentiator (filter) which is optimum in a more re- 
stricted sense. We shall then show that, even with the most 
optimistic assumptions about the marker system, no possible 
discrete differentiator can meet the specified accuracy and 
bandwidth requirements. 

A discrete filter, using 2// + 1 samples at time intervals h, 
may be described by 


+M 


7=—-M 


where x(r) is the sled position at time 7, and v*(t) is the ve- 
locity-measurement output of the filter. 

If it can be assumed that N(w) is white noise, we find that 
for M = 5 (eleven sample points) the first term of Equation 
[1] is minimized by 


[3] 


if we require that there be no measurement error [®(w) = 0] 
when N(w) = 0, and when the sled position is representable 
by an arbitrary quadratic equation in ¢ over any time interval 
(t — 5h,¢t + 5h)5; that is, the acceleration is constant. 
Proceeding as explained in Part I (1), we find that the trans- 
form of the discrete differentiation (or filtering process) de- 
scribed by Equations [2 and 3] is 
a 5 
(w) 55h [4] 
for w < m/h. 
We evaluate N(w) on the assumption that the errors in 
marker position are independent of one another. Let tlie 


5 Without this constraint, it is obvious that the noise error is 
minimized by W; = 0 (allj), and therefore G(w) = 0. 
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standard deviation of such errors be ¢,. There is also an error 
in measuring the time between successive markers. The 
actual timing operation consists of counting cycles of a fre- 
quency standard into a register, the counting being initiated 
and terminated by a pulse emitted the instant a marker is 
passed. We assume a counting frequency f, and note that the 
error in timing any interval may be incorrect by as much as 
1, ie., 6 = 1/f. However, it is equally probable that the 
time error is any lesser figure in the interval (—1/f, +1/f). 
The distribution of timing errors is then rectangular, with the 
mean square deviation 1/3f%. Converting this into an equiva- 
lent marker location error gives }(v/f)?, since v/f is the dis- 
tance traveled by the sled per cycle of the frequency standard. 
The total mean square equivalent error in marker location is 
therefore 


lyf \? 


Because the position sampling interval is h, and the errors 
ol successive samples (marker positions) are independent of 
one another, the form of N(w) can be taken to be that of white 
noise—perfectly flat—except that it is cut off at |w| = m/h 
and is zero for all larger |w|. But since its mean square value 
is o%, it follows that 


and therefore 
4 
| 

h me 
[7] 

lol > 

0, h 


Lastly, we must assume a value for A(w). In a typical 
case, the acceleration spectrum might be covered by a white 
noise spectrum of sufficient magnitude. A value of A(w) 
which is somewhat pessimistic but of the proper order of 
magnitude for a number of applications is 

A(w) = Ao? = 1000 (ft/sec?)?/rad/sec....... [8] 

To reduce the measurement error, we should consider only 
those frequencies below some cutoff frequency w,, such as the 
100 cps requirement. This is exactly equivalent to passing 
the measurement error through an ideally flat filter with a 
sharp cutoff at w,. Since &(w) > 0 for every w, this can only 
reduce the error. Further, the error increases very rapidly if 
w,/2r has a value comparable to 1/2h due to the frequency 
foldover effect of the sampling. We must have w, < m/h 
then, so that expressions containing sin hw can be approxi- 
mated by the leading terms of their power series expressions. 
By this series of steps we obtain successively from [1] 


55h 


WwW — 


j= 


and by approximating the sines as described above 


2 2 2 
30 


With a cutoff w, the mean square velocity measurement 
error becomes 


which in view of Equations [8 and 10] reduces to 
€,? = w.%h(0.1060 + 2930h*)............ [12] 
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To express this in terms of sled velocity and distance be- 
tween markers we note that this distance Az is given by 


[13] 


so that the final expression for the mean square error is | 


3A. 3 
0.10004 + 2930 (~) | 
v v 


We now estimate the rms error under the assumption that a 
maximum effort is made on a crash basis to instrument a 
22,000 ft track. If Ar = 2 ft, and o, = 0.002 ft, 11,000 
markers have to be surveyed and placed along the track with 
a mean deviation of only 0.024 in. Also, take f = 10mc. The 
cutoff frequency is taken as w, = 2007 (100 cps). Then, for 
v = 1000 ft/sec, 


[14] 


Equation [1], as well as those subsequent to it, have two 
separate terms, corresponding to the two sources of error de- 
scribed earlier. The first of these is associated only with the 
noise input (physically, marker location and timing errors), 
and may be called noise error. The second is due solely to the 
difference between approximate (discrete) differentiation and 
true differentiation; we shall denote this by the term ‘‘fre- 
quency response error.’ Let these types of errors be desig- 
nated by (e,?),, and (€,”),,, respectively. We have, then, from 
Equation [14] 


[16] 
and 
It is evident that changing the weighting coefficients W; in 
M 
v*(t) = + jh)............ [18] 


will alter the value of rms error. Consequently, the question 
arises whether, by increasing M and/or varying the W,’s, a 
smaller error than Equation [15] may be obtained. 

Consider first an increase in M (the above implies M = 5) 
so as to employ a greater number of samples. Since we have 
used a perfectly flat filter to a cutoff w,, we already have an 
infinite memory filter,® i.e., every sample obtained along the 
track is used. Thus, if such a filter is represented by 


= > Yat + ih)........ [19] 
there follows a | 
i=-—o@ jJ=-M 
from which 
v(t) = Viet +hh)............ [21] 
M 
[22] 


Now a change in V; implies a change in Y, W, or both. If 
Y is varied (without, however, affecting w, which is specified) 
the filter must depart from flatness in the region0 < w < w,. 
Suppose the variation is such as to compensate for the 


6 A filter of this type and the meaning of infinite memory with 
respect to such a filter may be found in Valley and Wallman (4), 
p. 722, and, more generally, on pp. 721-727. 
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Fig. 2 Block diagram of velocity measurement error when 
velocity is measured by two separate, independent devices 


drop in the frequency response due to W; and the limited 
value of M@. This would be a lead or differentiation filter 
which, while improving the frequency response error, would 
increase the noise error which is already excessive. Any other 
change in Y might decrease the noise error at the expense of 
the frequency response error, or possibly even increase both 
errors. 

Consider also the changes in W. But W has already been 
optimized with respect to the noise error, and therefore any 
change whatever in W must increase the noise error. 

Thus we have proved that any change in the discrete filter- 
ing coefficients must either enlarge the frequency response 
error, the noise error, or both. Since either error is in itself 
excessive (see Eqs. [16 and 17]), it follows that there is no 
system of weighting coefficients which can reduce the total 
error to 0.1 ft/sec on an rms basis with the parameters used. 

It should be noted again that the parameters—that is, the 
marker spacing of 2 ft, mean deviation in marker placement 
of 0.024 in., and counting frequency of 10 mc—are all marginal 
with respect to the techniques available and the time re- 
quired. Therefore, any marker system using feasible parame- 
ters will have still greater measurement errors than those 
quoted and considered excessive. 

The problem may also be approached from an alternative 
viewpoint, namely, what bandwidth can be achieved with the 
above parameters, given that the rms error is to be 0.1 ft/sec 
or less. Using Equation [14], it is seen that this cutoff fre- 
quency is w, = 59.5 rad/sec. 

It must be concluded, therefore, that a system of markers 
along the track, unaided by other means of velocity measure- 
ment, presents little promise of fulfilling the requirements set 


forth in Part I. 
3 Optimum Combination of Position and 


Acceleration Data 


We wish to find a lower bound for rms velocity measurement 
error when linear, non-time-varying filtering methods are 
employed to treat independently position and acceleration 
data to obtain velocity. The lower bound so calculated 
serves the dual purpose of: 

1 Determining whether it is at all possible to attain the 
desired accuracy with existent track coil and accelerometers. 

2 Evaluating specific filter scheme accuracy relative to the 
best accuracy which can be achieved by any method. 

Starting from more general considerations, we assume that 
the same time-varying quantity is to be measured by two 
separate, independent devices. The outputs of these de- 
vices are passed through filters of responses Y;(w) and Y2(w). 
Defining the error as being the difference between the meas- 
ured and actual quantity, we see that the block diagram of 
Fig. 2 is applicable. 

Suppose that the signal and both errors are stationary and 
mutually uncorrelated. In particular, let the input (signal) 
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spectrum be ¢,(w), and the spectra of € and € be ¢:(w) and 
¢2(w), respectively. Then, according to Fig. 2, the output 
(error) spectrum ®(w) is 


Pw) = [¥i(w)|? G(w) + ¥x(w)|? + 


— ¥i(w) — Y2(w)|? ow)... ... 
and the mean square error e? is given by pee, 
@= f B(w)dw. 4) 


In general, ¢1, $2 and @, are determined by the measur ng 
devices. The problem of minimizing e? must then be attacl.ed 
by choosing Y; and Y2 in an optimum fashion. This is exactly 
the Wiener filter problem. However, an arbitrarily long ‘ag 
can be permitted in the process of data reduction, so that filer 
realizability is not a problem. 

Turning to Equation [23], we see that ®(w) (and therefire 
€?) is minimized by Y, and Yo, which are real and non-neva- 
tive. Thus, we need only differentiate ®(w) with respect to 
the magnitudes of Y; and Y», set the derivatives to zero, sd 
solve for the magnitudes. Whenever reasonably accur: te 
instrumentation is employed, ¢,(w) >> ¢:(w) and ¢,(w) > 
¢2(w) so that the procedure described yields approximat«ly 


= 


gi(w) + dx(w) 


and 
y 
= 


The corresponding mean square error is 
gi(w) + 

The analysis will now be applied to velocity measurement 
of a rocket sled by means of markers spaced equally along the 
track in time, aided by an accelerometer mounted on the sled 
itself. 

Let ¢;(w) be taken as the error spectral density resulting 
from measurements by the markers alone. Ordinarily, ¢;(w) 
would be expected to contain terms resulting from aliasing 
effects, imperfect differentiation,’ and marker spacing and 
timing errors. 

Aliasing effects in the marker data will be due to sled mo- 
tion frequency components at frequencies higher than half 
the sampling frequency. If the time interval between suc- 
cessive markers is At, this is equivalent to position samples 
being taken with sampling angular frequency 27/At. Alias- 
ing of this data can be eliminated if the position samples con- 
tain only frequencies below z/At, that is, if it is possible to 
eliminate higher frequency components in some manner. 
These frequency components may actually be eliminated as 
follows: Pass the accelerometer output through a unity gain 
high-pass filter with low frequency cutoff m/At, integrate 
twice to obtain position, sample the integrator output at the 
same instant the position is sampled by the marker systems 
and subtract this result from the marker system samples. 
In this way all the aliasing error can be eliminated, except for 
a small residual entirely due to accelerometer errors at the 
higher frequencies. The latter effect is small (especially since 
the high-frequency accelerometer errors are integrated twice) 
and may be disregarded in this analysis. 

We treat next the error due to imperfect differentiation 
again assuming an interval At between successive track 
markers. The position samples may be passed through an 


[26] 


7 Ordinary differentiation cannot be applied to a signal con- 
sisting of sampled data. 
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interpolation filter flat to an angular frequency of w/ At, and 
cutting off sharply at that frequency. The resulting con- 
tinuous signal can be differentiated in ordinary fashion; aside 
from aliasing and spacing errors, the differentiation will be 
perfect up to angular frequency 7/At. In the calculations 
which follow, it will be assumed that outputs of all frequencies 
are available to the filter Y\(w). This is not strictly correct 
in light of the above. However, it will be seen that the re- 
sponse of Y,(w) cuts off well below w = w/At, and that only a 
truly insignificant portion of the total output is due to the 
track coil data at frequencies above m/ At. 

Marker spacing and timing errors may both be regarded 
as errors in position samples. In practice, the greater part of 
this error is attributable to sled bending and flexing, so that 
errors may be assumed independent from marker to marker. 
The effect is then the same as a band limited white noise 
position error, whose spectra density is Ato,?/m up a to fre- 
quency where is the mean deviation of marker posi- 
tion. When the position data are differentiated as indicated 
in the preceding paragraph, the resulting error spectral 
density of the velocity calculated from the marker position 
samples is 


ws 
ow) = [27] 
0, | otherwise 


For the power spectrum ¢2(w) of the integrated output of 
the sled-borne instrumentation accelerometer, we consider a 
zro offset, a scale factor error p and a noise spectrum R(w) 
generated on the acceleration output within the accelerometer. 

The zero offset contributes to the power spectrum as a delta 
function at the origin. By Equation [25], Y2(0) = 0, and there 
is no net error due to the offset. 

The scale factor error may be virtually eliminated by the 
first correction method proposed in section 4 of this paper. 
Consequently, the validity of these calculations is not compro- 
mised by omitting consideration of the extremely small in- 
accuracy which remains after scale factor compensation. 

Lastly, the accelerometer noise error is taken to have a 
spectral density of R(w), so that the velocity which is derived 
by integrating the accelerometer data has an error spectral 
density of 


If the optimum filters Y,(w) and Y2(w) are used, the total 
spectral density is computed to be 


In the absence of further information, it is reasonable to 
suppose that the acceleration spectrum is white. We assume 
then, that 


R(w) = [30] 
With this substitution, the mean square error is determined 
by straightforward integration. We have 
= 5 [31] 


2/2 


where @,” is the mean square velocity error. 

In order to evaluate the mean square velocity error, one 
chooses values of N? and ¢,? which are well within the range 
of achievement, and require no development time and effort. 

Marke > avai ft intervals to an accuracy 

arkers are av ailable at 100 ft inte rvals to an accuracy of 
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0.01 ft rms. Then, if v represents sled velocity in ft/sec 


100 
Ait = — sec 
= 10~4 (ft)? 
ss N? = 0.0318 (ft/sec?)?/rad/sec (see footnote 8) 


and the rms velocity measurement is 


v 8 

where v is in ft/sec. It is seen that the rms velocity error is 
below 0.1 ft/sec for any velocity in excess of 15 ft/sec. Ata 
speed of 200 ft/sec, this error is only 0.073 ft/sec. By the 
time 2000 ft/sec velocity is reached, the rms error has dropped 
to 0.054 ft/sec. Furthermore, these figures apply to an arbi- 
trarily large bandwidth, since the upper limit of the integral 
of the error spectral density was taken to be infinity. 

It is possible to speak of the frequency response of Y,(w) and 
Y.(w) in the usual sense. For example, by substituting the 
expressions for ¢)(w) and ¢2(w) into Equation [25] for ¥i(#), 
we see that 


which is a low pass filter with a cutoff frequency® of 1.8 v'/4 
rad/sec and an attenuation of 24 db/octave. 

Since Y.(w) = 1 — Y,(w), it follows that Y2(w) is a high 
pass filter, with the same cutoff and attenuation characteristics. 

It was stated earlier that, while the filter Y\(w) on the 
track coil data exhibited response beyond half the sampling 
frequency, the response was so small that the sharp cutoff 
assumption on ¢;(w) (see Eq. [27]) was of little consequence 
in shaping Y,(w). Since half the sampling angular frequency 
is given by 7v/100, a substitution of that figure into Equation 
[34] shows that the response at that frequency is down by a 
factor of more than 10°/v', or about 10~4 at a velocity of 1000 
ft/sec. Response inaccuracy at higher frequencies is still 
smaller, related inversely as the fourth power of frequency. 

The filters Y,(w) and Y2(w) have to be closely matched, 
otherwise the third term of Equation [23] becomes pro- 
hibitively large. This consideration (together with the dis- 
crete nature of the marker data) suggests that a discrete 
(digital) filter of the form 


+N 


2 


is appropriate, with u and w representing the outputs of the 
marker system and sled-borne accelerometer, respectively, 
and with U;, W,, N and 6t chosen to give a discrete represen- 
tation of Y;(w) and Y.(w). 


4 A Practical Method of Obtaining Velocity 
From Position and Acceleration Data 


We shall now present a constructive method of velocity 
measurement. The computational procedure necessary to 
accomplish the measurement of velocity will be given ex- 
plicitly, so that the method should prove directly useful. 

An error analysis of the scheme described here constitutes 


8 This figure has been chosen to be consistent with the noise 
values used in section 4. The actual choice of this number is ex- 
plained in greater detail there. 

® The cutoff frequency is defined as that frequency for which 
the attenuation is 3 db. 
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a part of this section. In this fashion, the adequacy of the 
technique is established for a measurement accuracy of 0.1 
ft/sec rms over any bandwidth. 

The method of this section is not claimed to be the optimum 
method. It does have the virtues of computational simplicity 
and accuracy close to that of an optimum method. Indeed, 
without an exact knowledge of the sled acceleration spectrum 


and the sled-borne accelerometer error spectrum it is im- | 


possible to design a truly optimum filter. Furthermore, the 
work of section 3 indicates that, even under simple assump- 
tions on the two spectra, the computational mechanization of 
the filter scheme is rather complex. 

A brief description of the velocity measurement method 
under discussion will precede its analytical representation. 
If the sled-borne instrumentation accelerometer has bias and 
scale factor errors, these may be estimated by comparing the 
doubly integrated accelerometer output with track coil data 
from three coils along the track. Corrections based on the 
bias and scale factor error estimates may then be made to 
the accelerometer output. The correction just outlined will be 
called the “gross correction.”’ Its calculation is so quick and 
simple that it may serve to provide ‘quick look data” avail- 
able almost immediately after the end of a run. 

The gross correction may result in data adequate for many 
purposes; however, to achieve an accuracy of 0.1 ft/sec or 
better, further effort is required. A final correction is made 
by comparing the average velocity obtained by differencing 
successive track coils with the average velocity obtained from 
accelerometer data on which the gross correction process has 
already been performed. This entire process is illustrated in 
Fig. 3. 

The actual output of the accelerometer a*(t) may be ex- 
pressed as 


a*(t) = Ko + (1 + Kida(t) +a(t)......... [36] 


where Kp is the bias error, K; the scale factor error, a’(é) all 
other accelerometer errors (including particularly random 
noise), and a(t) the true sled acceleration. 

In order to perform the gross correction, it is necessary to 
estimate Ky and K,. If the sled starts at time zero and 
passes by the jth and kth track coils at times ¢; and ¢,, respec- 
tively, the doubly integrated (position) accelerometer output 
is 


r*(t;) = + Ki) 2(t;) + 2(t;)....... [37] 


and 
= Ko — + (1 + Ki) (te) + 2'(te)...... [38] 
at times ¢; and ¢,, where x(t) is the true sled position. 


The track coil passed at t; is thought to be at position 2; 
rather than x(t;) because of surveying errors, sled flexing and 


vibration. Hence, there is a track position error €; which is 
equal to 


The track error €; should also include an error which results 
because the time of passing a given track coil is not measured 
exactly, but only to a given increment. The SNORT track 
instrumentation is capable of measuring this time to the near- 
est microsecond, so that the resulting error is small compared 
to an estimated one sigma track error of 0.01 ft. 


Equations [37 and 38] yield 


be. 4) K,. 


4, 


(t)—x 
x(t) 


itl 
i 


x. t 
i isl 


Fig. 3 Block diagram of practical method for obtaining velocity 
from acceleration and position data 


= A(j, k) = = 


where 


In Equations [40 and 41] we know the values of 2x*(/), 
2;, t, t; and t, This enables us to construct estimates A, 
and K, for Ky and K,, respectively. We shall use 


= 2 

—— — 1 [44 


If the error in estimating Ko is defined as dK, = Ky — Ky, and 


similarly for K,, one may write 
dK, = 


— + Ki) + — 


AG, k) 


= 
= AG, k) — + Ki) + — 


We shall assume that a’(t) has zero mean (any constant non- 

zero mean is absorbed in Ko), as do track coil errors €; and «. 

Then the estimates Ky and K, are unbiased. 
At any time ¢ the true velocity of the sled is 


v%(t) — Kot — v'(t) 


from itl Our estimate i(t) of the velocity shall be 
v*(t) — Kot 
v(t) = [48] 


which is easily computed after Ky and K, are calculated 
from [43 and 44]. The error in this estimate of v(t) is 


(dK, )o(t) + (dKo)t + 


In our application, a so-called 1 per cent accelerometer or 
better is required. Then |K,| < 0.01, and likewise K, will be 


Ky = — + Ki) — 2(t;)] — — + Ki) — [40) 
(t.) — + Ki) — — t2[x*(t;) — + Ki) — —1........ [41] 7 
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small. On the basis of this argument it is legitimate to ap- 
proximate da(t) by 


di(t) = — (dR,)u(t) — (dKo)t — v'(t)......... [50] 


As stated previously, the gross correction described above 
may suffice for some purposes, and in any case proves useful 
for quick look data. Therefore, we shall calculate the rms 
error in the measured velocity after the gross correction has 
been applied. 

Since Ay and K, are unbiased estimates, and a’(t) has zero 
mean, the mean square value of da(t) will be identical with the 
variance.!° To perform the required calculation, di(t) is 
squared, and its expectation taken. We shall denote “expec- 
tation of’ by prefixing the letter F to the appropriate quan- 
tity. Now 


an error of less than 1.05 X 10~‘ (ft/sec)? in R,,», and an error 
of less than 10~° (ft)? in R,,. 
If the indicated calculations are carried out 
E(ds*(t)] = 0.14 + 0.0962 — 0.00252 /sec)?. . . [58] 
whenever ¢t < ¢;. Fort; < t < &, we have instead 
E[do?(t)] = 0.1¢ + 0.0356¢? + 0.0035¢* (ft/sec)?. . [59] 


In particular, if ¢ = t; = 10 sec, the rms value of do(t) is 


WV 
and 


V 


E{ = + CE[(dK,)?] + Ef [v'(t)]?} + + + . (51] 


The remainder of the work will be greatly eased if a few 
assumptions are made at this point. Since these assumptions 
will correspond closely with fact, no apology is necessary. 
The great preponderance of the track coil error is due to 
flexing and bending of the sled itself, so that €; and e are 
orthogonal for any 7 # k. The track coil error is also as- 
sumed to be uncorrelated with the accelerometer error. 
With regard to the latter, Ko, K, and a’(t) are all mutually 
(pairwise) uncorrelated. 

An additional supposition applicable only to the gross cor- 
rection analysis is that E(e;?) < E{ [x’(t;)]*}. The justifica- 
tion lies in the fact that the rms value of ¢; is about 0.01 ft, 
while x’(t;) has an rms magnitude of 8.2 ft in the representa- 
tive example which is to follow, and 2’t,) has an even larger 
rms value. Consequently, the calculation of the expectations 
in [50] will omit all considerations of track error. Since 
E{ [x’(t) }?} tends to grow as ¢*, the same approximation is not 
appropriate to the final correction in which we deal with time 
intervals of the order of 0.1 sec. 

Reference to Equations [45, 46 and 51] will show that terms 
such as E[zx’(t,)x’(te)] occur in the evaluation of do(t). To 
show how such terms are evaluated we define first the co- 
variance FE [a’(t)a’(te)] by 


t2) = 
and similarly 
R..(h, te) = [53] 
R,,A(h, ts) = [54] 
Given R,,(ti, te), these covariances may be computed as 
follows 


t te 


h 


te 
R,,(h, te) = f, T2)dt [57] 


As an example, we shall calculate dv(t) for a rocket sled 
traveling along a 20,000-ft track. We shall assume that the 
sled accelerates at 200 ft/sec? for 10 sec, and then decelerates 
at the same rate for the same length of time. Then, with z; = 
10,000 ft, and x, = 20,000 ft we shall have ¢; = 10 sec, and 
i, = 20 sec. The accelerometer has an assumed noise error 
a'(t) of 3.2 ft/sec? rms, flat from 0 to 50 cps; this implies a 
spectral density N? = 0.0318 (ft/sec?)?/rad/sec. We may 
thus write R,,(h, t2) = 0.0318 6(4; — te) as an approximation 
to white noise over the entire frequency spectrum. The dif- 
ference between the exact R,, and the one just given implies 


10 This statement is true if the approximation of [50] is valid. 
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when t = t = 20sec. The rms values of dKyand dK, are 
0.115 ft/sec? and 0.00167, respectively. Without the pro- 
posed correction, a 1 per cent instrument might have Ky = 
3 ft/sec? and K, = 0.01, so that errors would be in excess of 60 
ft/sec! 

Since only three track coils have been used to correct the 
accelerometer, it is evident that better results may be achieved 
through further corrections made by use of track posi- 
tion data obtained from other track coils. The procedure 
which is proposed here consists of correcting the singly inte- 
grated accelerometer output by the difference between average 
velocities between two track coils as determined from track 
coils and accelerometer over that interval. That is, in the 
time interval t; < t < t+; our new estimate of sled velocity 
is v(t), where 


0) = + + (62 


in which the difference between average velocities A; is given 
by 


“az Ai = zi) wee | [63] 


and At; is defined by 


The error in our estimate is then 


Git — dX(ti+1) — dX(t;) 
do(t) = Ai, + di(t) [65] 


It will be noted that aliasing does not appear in any of the 
error terms of [65] regardless of the frequency components 
of the rocket sled motion. This is to be expected, because the 
track coil data are used only to determine the average velocity 
between two track coils as shown by [62 and 63]. 

For the purposes of our error analysis we shall assume that 
the ith and (i + 1) track coils are not identical to either the 
jth or the kth (on which Ko and K, have been estimated). 
Since there are 200 coils along the track, this assumption will 
generally hold. If the assumption is not true, F { [do(t)?} will 
contain terms in addition to those that we shall calculate. 
However, these terms are small for practical ¢;, so that our 
error estimates will apply along the entire track. 

The mean square error of velocity measurement is most con- 
veniently expressed as 


2E 
(At;)? At; 


tiga 
E{do(t)]*} = f + Et [do(t)]*} 


E{[dx(tis1) — dx(t;)]?} 
(At;)? 
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Here da(t) is given by [50]in terms of dKy, dK, and v'(t). In 
turn, dK, and dK, are expressed as functions of z’(t;) and 
z’(t,) in [45 and 46]. Making the firstindicated setofsub- 


— At) 


(At;)? 
= tins + bi , 
2 4 [o(t) — — ) + (t — —>— )\ — 
bashed ti 

+ if f dn Ai J. — 2R,,ft, »|t + R,,(t, t)..... [67] 
where for brevity of notation, we write the average velocity = °°; 
over the interval (t;, t:+1) as 


In some cases it may be sufficient to obtain an upper bound 
for E[dv(t)]?. The following inequalities are useful for this 
purpose 


Curve Computed from Eq 39 with 


= 300 ft/sec* 
| v(t) = [70] Sled Velocity, ft/sec 
: é : Fig. 4 Upper bound on rms velocity measurement error vs. 
where max is the maximum expected acceleration. Then sled velocity 
1 tina 


Equation [70] may be obtained from [66] by repeated use of method presented here promises sufficient accuracy to ap- 

the Schwarz and Minkowski inequalities (5). proach the specifications. Indeed, the figures presented here 
A numerical example will indicate the order of accuracy refer to upper bounds based on pessimistic accuracy estimates, 

attainable by the velocity measurement scheme under discus- so that greater precision should be realized in practice. 

sion. We shall again assume that the sled travels down a 20,- 

000-ft track, first accelerating at a rate of 200 ft/sec? for a 

period of 10 sec, and then decelerating at the same rate for References 

the een" period of time. The rms position erwleage the track 1 Beutler, F. J. and Rauch, L. L., “Precision Measurement of 

coils will be taken as 0.01 ft, and the acceleration noise error Supersonic Rocket Sled Velocity, Part I,” Jer Proputston, vol. 

in the accelerometer as 3.2 ft/sec? rms uniformly over a 50 27, 1957, pp. 1021-1024. 

cps bandwidth. We shall again base our calculation on white © 2 Milne, W. E., “Numerical Calculus,” Princeton University 

noise of the same spectral density. ‘The accelerometer may also Press, 1949, chapter IV. 

have bias and scale factor error; these do not enter into our 3 Johnson, K. R., “Optimum, Linear, Discrete Filtering of 

equation as per the discussion on the first correction. The Signals Containing a Nonrandom Component,” JRE Transac- 


maximum acceleration, dmax, is 300 ft/sec? for the purpose of or Information Theory, vol. IT-2, no. 2, June 1956, pp. 


f 4 Valley, G. E., Jr. and Wallman, H., “Vacuum Tube 

ig. 4 gives an upper bound or rms ve ocity measurement Amplifiers,’’ McGraw-Hill (Radiation Laboratory Series), New 

error by the method of this section. Equation [71] has been York, 1948, pp. 721-727. 

used to find the numerical values for the graph. 5 Graves, L. M., “Theory of Functions of a Real Variable,” 
The error analysis demonstrates that the measurement McGraw-Hill, New York, 2nd ed., 1956, pp. 233-234. 
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Some thirty upper atmosphere rocket soundings have 
been made in the Arctic so far in the International Geo- 
physical Year to measure temperature, pressure, density 
and winds. Aerobees in the various models and Nike- 
Cajuns have been the vehicles. The experiments have 
included the freely-falling sphere for measuring densities, 
the rocket-grenade experiment for measuring pressures 
and densities in the supersonic flow, yielding the ambient 
conditions by theoretical conversion. Data from a number 
of these firings have been reduced and compared to those 
taken at other latitudes. The preliminary analysis indi- 
cates different temperature, pressure and density distri- 
butions than elsewhere. Also, these parameters seem to 
be more sensitive to the time of measurement than at 
lower latitudes. 


Introduction 


HE meteorological parameters of temperature, pressure, 

density and winds have been and are being measured in 
the Arctic during the International Geophysical Year by a 
variety of techniques developed over the previous ten years 
by various groups in this country that have been carrying out 
the upper atmosphere rocket research program. The scope 
of the program has been described elsewhere (11).* Tables 
1 and 2 are summaries of the meteorological firings. 


Table 1 Meteorological rocket soundings (Arctic) 
Experiment No. firing Location Time 
sphere 8 North Atlantic 

and 
Churchill = W,D 
grenade 10 Churchill = 
aerodynamic 6 Churchill 
gage 3 Churchill WS 
D,N 
USSR = 6 Franz Joseph 
Land W,D 
W = winter, S = summer, D = day, N = night 


The IGY is only a year old at present; most of the rocket 
research groups are still in the middle of their firing schedules; 
the results, in terms of reduced data that are available at this 
time, are quite limited. The scientists themselves have not 
had the opportunity to publish their information. There- 
fore, this paper is intended to review the techniques that are 
being used in the Arctic to measure the meteorological param- 
eters and to present a few of the preliminary results that the 
working scientists have so kindly made available. Other 
published data have been included when pertinent. 

In the planning of the IGY rocket program for Fort 
Churchill, a guiding principle was that proven rocket-borne 
techniques would be used. Each of the methods used, 


1 Presented at the ARS Semi-Annual Meeting, Los Angeles, 
Calif., June 9-12, 1958. 

3 Chief, Space Research Instrumentation, Applied Physics 
Division. Member ARS. 

* Numbers in parentheses indicate References at end of paper. 
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Soundings in the Arctic 


whether sphere, aerodynamic, grenade or gage, has been tested 
at White Sands Proving Ground and has yielded significant 
data (la,b). Thus, there is added the significance of having 
measurements made with the same method in different geo- 
graphic locations. 


Descriptions of Methods 


There follow quite simplified descriptions of the various 
experiments for measuring the meteorological parameters of 
temperature, pressure, density and winds at rocket altitudes. 
References (1) give more thorough discussions. 


The Freely-Falling Sphere Experiment (2) 


This method of measuring the density of the atmosphere 
has been developed by the University of Michigan group 
under L. M. Jones, first under Signai Corps sponsorship using 


Table 2 Meteorological rocket soundings (Arctic) 
Alt. 
IGY no.! Date? Location Experiment km 
AM6.08 20Oct.— North 
through— through— Atlantic sphere (D) 90 
AM7.12 20Nov.56D 
AM6.31 200QOct.56D Churchill aero(T) 113 
AM 2.21 23 Oct. 56N aero (T) 145 
SM1.01 12Nov.56N grenade (T, W) 68 
NN3.12 17Nov.56D gage (P,D) 210 
21 July57N grenade (T,W) 92 
8M1.03 23 July 57N grenade (T,W) 87 
NN3.13 29 July 57D gage (P,D) 210 
F 
SM1.04 12 Aug.57D grenade (T,W) 74 
SM1.05 19Aug.57N grenade (T,W) 93 
SM2.06 25 Aug.57D grenade(T,W) 130 
AM4.01 1 Sept.57D eS aero (T) 160 
SM1.07 11 Dec.57N grenade (T,W) 85 
SM1.08 14Dec.57D . grenade (T,W) 97 
AM6.02 25Jan.58D a aero (T) 160 
SM1.09 27Jan.58N = grenade (T,W) 98 
SM2.10 27Jan.58D eS grenade (T,W) 152 
sphere (D) 
AM6.36 27 Jan.58D aero (T) 130 
AM6.03 29Jan.58N sphere (D) 160 
AM6.37 24Feb.58N aero, T) 144 
NN3.14 24 Feb.58N gage (P,D) 206 
F 
AM6.05 4 March58 sphere (D) 170 
AM6.38 24 March58 = aero (T) 137 
1TGY no.: A code number indicating the sponsoring 
agency, the instrumenting agency and the type and num- 
ber of the vehicle used. 
2The D or N after the year indicates day or night 
respectively. ‘‘Aero’’ is the aerodynamic technique of 
N. Spencer, Univ. of Mich. ‘‘Sphere”’ is the freely-falling 
sphere technique of L. M. Jones of the Univ. of Mich. 
“Grenade” is the rocket-grenade experiment of W. G. 
Stroud of the Army Signal Corps. ‘Gage’ is the multi- 
ple-gage technique of H. LaGow of the Naval Research 
Laboratory. P, T, D and W are pressure, temperature, 
density and winds, respectively. 


nal Re: velopment Laboratory, F 
saboratory, Fort Monmouth, N. J. 
I } 
| 
| i 
va 
i 
| i 
q 
q 
q 
| 
7 
= 
). 
= 
| 
| 
; 
i 
i 


 weme ACCELEROMETER 
Fig. 


Aerobees and later under Air Force sponsorship using Nike- 
Cajuns. The present procedure consists of carrying aloft a 
7-in. diam, 9-lb sphere, ejecting it from the Nike-boosted 
Cajun, and measuring its drag by means of an internal acceler- 
ometer as it falls through the atmosphere. 

The basic equations are 


where 
D = density of air 
v = velocity of sphere 
Cp = drag coefficient 
ANTENNA 


F=M(g-a) 
DRAG=3Dv'C,A 


BOBBIN- 


SR- SOUND RANGING 

D -DOVAP 

L -AEROBEE LAUNCHER 
G -GRENADE EXPLOSION 


y 


ALPHATRON 
NOSE GONE 
ALPHATRONS 


The deceleration is measured; the drag coefficients are 
known from scaled wind tunnel tests at low pressures; the 
velocity is obtained from the free-fall from peak, and the 
other known quantities permit the calculation of D. 

The altitude range over which the experiment is effective 
depends on the height from which the sphere is dropped, the 
ratio of area to mass of the sphere, and the sensitivity of the 
accelerometer. The present configuration yields data be- 
tween 10 and 90 km. The accuracy is about 10 per cent 
over this altitude range, the primary limitation being the 
uncertainty in the coefficient of drag over the range of Mach 
and Reynolds numbers. Pressures and temperatures may 
be derived from the computed densities by means of the 


hydrostatic equation. 
| 


The temperatures and winds of the upper atmosphere have 
been measured by sound-ranging on successive explosions 
created by carrying aloft and ejecting charges from Aerobees. 
The present experiment is conducted jointly by the Army 
Signal Corps group under W. G. Stroud and the University of 
Michigan group under L. M. Jones. The fundamental parani- 
eters measured are the positions and times of the explosions 
in space, and the times and angles of arrival of the successive 
sound waves at the ground. The experimentally determined 
altitude limit for the method is about 90 km above which 
sufficient sound energy cannot be coupled into the atmosphere. 

The average temperatures in the layers between grenades 
are calculated from the relationship 


The Rocket-Grenade Experiment (3,4) 


if RT x (T)'/2 [3] 
where 
C = the velocity of sound in the layer ‘- = 
y = ratio of specific heats of air 
R = gas constant . 
7 = temperature a 
M = average molecular weight of air | 
K = constant over the altitude range | 


The average wind speeds and directions in the layers between 
grenades are determined by measuring the differences be- 
tween the true positions of the explosions and their apparent 
positions as given by the angles of arrival of the sound waves 
at the ground. The accuracy of the measurements is limited 
by that with which the times of arrivals of the low frequency 
sound waves at the individual microphones can be measured. 
The average temperatures have an accuracy of +5 C, the 
winds, +10 m/sec and +18 deg in azimuth as determined 
in previous firings. 


= 
The Aerodynamic Experiment (1,5) 


By means of the Taylor-Maccoll (6) theory concerning 
supersonic flow around right-circular cones and the Kopal 
tables (7) for such cones in yaw, it is possible to determine 
ambient temperatures from a supersonic rocket by measuring 
the ram pressures and cone sidewall pressures by appropriate 
gages. Spencer of the University of Michigan has applied 
alphatrons for such pressure measurements using both Aero- 
bees and Nike-Cajuns. Pressures at the nose tip and roughly 
0.8 of the distance back along the right-circular cone are 
measured. 

Thus, the pressure P, at a distance x back along the cone 
is given by 


where 7’ is the ambient temperature, v the speed of the missile, 
a and @ its angles of attack and roll, and m and + the molecular 
weight and ratio of specific heats, respectively. For gages at 
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different distances, x and y 


P, f (z, v, a, 6, m, Y) 


= 


For any rocket flight, all these quantities are known or deter- 
mined point by point except 7’, which may be computed 
point by point along the trajectory. The pressures may be 
derived from these temperatures using the hydrostatic equa- 
tion. The temperatures obtained over the region 40 to 90 
km are accurate to +10C. 


The Gage Methods (8,9) 


LaGow and his colleagues of the Naval Research Labora- 
tory have applied the earlier work (8) of NRL to Aerobee-Hi 
borne measurements in the Arctic. These measurements of 
pressure and density (from which the temperature may be 
derived) are carried out as follows: A family of pressure 
gages is scattered over the surface of the rocket measuring 
stagnation, ambient and cone pressures. Bellows, Pirani 
and Philips gages and the Havens-cycle gage measure the 
pressures in different pressure ranges. See Table 3. Roll 


Table 03 of sensitivity of various pressure-density 
gages 


= Range 
Gage A = mm Hg 
Bellows 760 to 20 
Pirani (conductivity Zto3s X 
Havens-cycle 10% to 10-5 
Philips (ionization) 103 to 10-6 
Alphatron 10% to 107? 


modulation appearing on the side cone gages also permits 
the determination of the pressure at higher altitudes (1077 
mm). The distribution of the gages along the walls of the 
vehicle is critical, since ambient pressures are measured di- 
rectly by various gages located about 7 calibers behind the 
stagnation probe. These Aerobees are sealed throughout 
to prevent outgassing and permit the extension of the meas- 
urements to the limits of gage sensitivity. 

The density is determined from the Rayleigh supersonic 
pitot tube formula which takes the form 


0.144 P, — 0.066 P, 
D= | [6] 
J 
where 
P, = stagnation pressure 
P, = ambient pressure 
V = total velocity 1 


This multiple-gage approach has made possible measureme nts 
to much higher altitudes than the lan techniques. Its 
accuracies, as given by NRL, are 


Altitude, km Error, per cent 


Pressure data 7 23-80 5-15 
110-150 20-50 
30-80 2-10 
110-300 20-200 
Other Methods 


Two other methods of measuring the density at rocket 
altitudes should be mentioned. Aerobee flight NN 3.17 
by NRL, which was instrumented with time-of-flight mass 
spectrometers for gas and ion composition (10), obtained 
density data because the rocket’s altitude stayed within 54 
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Fig. 4 The gage method 


deg of the trajectory during ascent, and the spectrometer at 
the tip yielded density data between 110 and 170 km. These 
are winter night data. Also, it has been observed in the 
grenade experiment that at the time of each explosion a 
period disturbance appears on the Dovap cycle record. 
This disturbance is apparently a measure of the velocity of 
propagation of the shock wave from the explosion which is a 
function of the ambient density—thus the density at each 
grenade altitude may be determined. Bartman of the Uni- 
versity of Michigan is at present working out this technique. 
A specific problem in the comparison of the rocket data is 
the time variations. Simultaneous measurements with the 
different experiments is highly desirable but difficult, because 
of the limited space on the vehicles and the limited facilities 
for launching and tracking on range. However, on Jan. 27, 
1958 at noon, an Aerobee carrying the rocket-grenade ex- 
periment and a small sphere was fired. Twenty-nine minutes 
later a Nike-Cajun carrying the aerodynamic experiment 
was successfully fired. A comparison of these results from 
three techniques will be particularly interesting and valid. 


The USSR Method 


The USSR is also conducting Arctic and, in addition, 
Antarctic meteorological rocket sounding programs (14) using 
their meteorological sounding rocket (15,16) to measure tem- 
peratures and pressures up to about 110 km. Between 
November 1957 and March 1958, six rockets were launched 
in Franz Joseph Land and four near Mirnyy in the Antarctic. 

As described in the literature, the instrumentation consists 

of thermal (Pirani?) and membrane manometers for measur- 


ae ing pressure, and an electric resistance thermometer for 


temperatures. A unique feature of their technique is the use 
of a parachute throughout the entire downward leg of the 
.: trajectory to orient and slow the separated instrumentation 
section in its fall. A telemeter is used to recover the data 
which apparently are now being processed. 

Previously published USSR results for the upper-middle 
latitudes (16) show general agreement with the U. S. liter- 
ature in this field. 


Results 


Of the twenty-seven firings, including the pre-IGY group, 
that have been conducted so far, data from only about ten 
have been reduced or even partially reduced. The respon- 


_ sible scientists have had only limited opportunities to publish 


their results. The data presented in the following paragraphs 
must be considered preliminary. Except for the grenade ex- 
periment which is the author’s responsibility, the inclusion 
of data in this paper is due to the cooperation of the respon- 
sible scientist. 
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The Sphere Experiment latitudes up to 80 km altitude have been obtained by the Uni- 
In aseries of five Nike-Cajun firings in the North Atlantic in ot r 
November 1956 rt of th IGY test ieee shown in Fig. 5. An initial conclusion is that in this region ‘ 
dist th of the atmosphere the densities are lower than the White 
Sands measurements which are the basis for the ARDC mode] 
‘ atmosphere. Also plotted on the graph in the region 90 to 
10 ———— 100 km are several points from Aerobee NN 3.12, fired at \ 
DENSITIES the same latitude at about the same time (Nov. 1956). The 
data are consistent. 
— ARDC MODEL These sphere samplings taken at latitudes of 49, 58 and 66 
107 moog — ooAVER. SPHERE deg N when compared with the sphere data from Wallops é 
+ X-X NN 3.12 = Island (38 deg N) and White Sands (32 deg N) suggest « X 
+ strong latitude variation in the densities at 50 km of about " 
> 2 per cent per deg of latitude. This is greater than we might , 
tl 
The Grenade Experiment 
sie ys The ten planned Aerobee-grenade firings of the joint Signal a 
-— = Corps-University of Michigan (Jones) program at Churchi'! tt 
= = have been successfully completed. Five winter firings and | 
Pg TAN five summer firings were made. Some were daytime, some th 
. \ nighttime; on Jan. 27, 1958 firings were made at midnight E 
> ad = = and at noon, to measure the diurnal variation of temper- a 
‘a = atures and winds up to 90 km. a! 
a x The reduction of the grenade data involves the determin- - 
z Y ation of the positions of each of the grenade explosions from in 
= io? = the Dovap records and the calculation of the angles of of 
=" arrival of the sound waves. The calculation of the temper- 
. atures and winds (12) from these elements has been pro- 
v grammed on an IBM 650 computer so that the final steps of 
io the reduction go quickly. But the Dovap reduction for th 
aoa =~ positions is a matter of hand counting of the Doppler cycles, fo! 
_ a = some 100,000 of them per firing, which takes several months én 
mp x Data from firings SM 1.01 (Nov. 1956) and SM 1.02 (Jul) fre 
7 wv 1957) have been reduced, and the preliminary results are Nj 
0 10 20 30 40 50 60 70 80 90 100 110 120 shown in Figs. 6 (temperatures) and 7 (winds). We are not : 
ALTITUDE (Km) 
Fig. 5 Densities—sphere 
sid 
100 ZA per 
suc 
TENT. ST'D init 
90 a SMI.OI is | 
o SMI.02 | the 
0 tha 
80 a RAWIN int 
x tha 
60}+~1 
Q = the 
ay 
» the 
= dat: 
- a : 
30 
¥ ‘he The 
Pan 
20 20 N 3.12 
A= 50M/ altit 
L* 10M/S of t 
5M/S 
| = Jucy | 
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Fig.6 Grenade data—temperatures Fig. 7 Grenade data—winds as tl 
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The temperature data show a difference between the sum- 102 a 12 10-2 
mer and winter distributions throughout the entire altitude 
region. At 50 km, a maximum in the temperature distribu- 4 ++ eee ASCENT 
tion, the spread is about 20 C, a significant amount. The LL TN xxx DESCENT : 
winter measurement is characterized by low temperatures | ESTIMATED AVERAGE 10-3 
below 50 km and large scatter from layer to layer. The H 
summer sounding is, by contrast, quite smooth and shows an w OF PRESSURE CURVE, 
extremely low temperature at 80 km of 165 K. This is the Vv ASS DENSITY CORRECTED FOR 
lowest temperature we have measured in the atmosphere. 10-4 \ | ASSUMED WIND OF 80 

The temperature data obtained from radiosondes flown a Ss=a=a WASEC. FROM SOUTH ‘a 
at the times of the rocket soundings are also plotted in the E eS Te =e et saa : 
figure. The soundings overlap in the 30 km altitude region, S AAAS Tq + 
and the data match quite well. : 

The grenade winds, shown in Fig. 7, have the character- 
istic summer-winter difference we have observed in the lower a tet WS tt : 
latitude measurements (4). The winds are weak and from re } me 
the east in the summer, very strong and from the west in the ry ttt pe 
winter. We have every reason to believe the data at this | | 
time, even the 150 m/sec velocity at 58 km. However, ° Ss 5g} 
the analyses are preliminary. +7: 

An additional finding of the Churchill firings is the fact } mt ats Mott 
that above about 90 km the experiment does not work. | | | | 
Even with the 4-lb explosives, sufficient energy cannot be Hoseas EEEy 
coupled into the atmosphere to permit sound arrivals at the t OWN: = a > 
ground. This is consistent with Schrodinger’s early work 
on the absorption of sound in the atmosphere, since the mean- | | Eonar 08 
freepath at these altitudes is about equal to the wave lengths 80 90 100 110 120 130 40 150 160 170 180 190 20 20 220 
of the sound waves. ALTITUDE (KILOMETERS) 

7 Fig. 8 NWN 3.12 densities 
The Aerodynamic Experiment 4 

As do all other experiments except the sphere experiment, 
the aerodynamic method needs detailed trajectory data be- 4 
fore the alphatron measurements of the soundings can be 0 \ 
converted to atmospheric parameters. At present no data 
from the six aerodynamic experiments on both Aerobees and i 
Nike-Cajuns have been made available. SE 
The Gage Experiments ROCKET PANEL 

Again, since the accurate trajectory data required for the E 4 
complete analysis of measured parameters have taken con- ‘ — 
siderable time and effort, the data reduction of the gage ex- i h 7 


periments is incomplete and preliminary; however, three 
successful firings have been carried out at Churchill. An 
initial reduction of data from NN 3.12 fired in November 1956 
is plotted on Fig. 8, densities vs. altitude. Again throughout 
the entire region of the measurement, the densities are lower 
than those previously measured at White Sands (13). This 
is true up to about 200 km. The group of points in the lower 
right-hand corner are single measurements taken at the peaks a 
of the other rocket trajectories when the velocities were mini- 
mal. Summer day (SD) densities are a factor of 20 higher io7E 
than those of winter night (WN), and winter day (WD) 
a factor of 10 higher than winter night. The confirmation of 
these results will provide an extremely interesting insight into ; 
the dynamics of the high atmosphere. The winter night = 
data are from NN 3.17. 

Another measurement of interest was the ion-composition 
spectrometer firing of November 1956 (NN 3.17) which, be- 100 120 140 160 180 200 220 
cause of favorable altitude throughout the upper leg of the ALTITUDE (KM ABOVE SEA LEVEL) 
trajectory, yielded density data between 100 and 200 km. Fig. 9 NRL-48 (NN 3.17) densities 
These data are plotted on Fig. 9 with a comparison with the 
Panel data (13) and the average ascent-descent data of NN - a 
3.12 (Fig. 8). If real, the day-night difference in densities at 


altitudes above 160 km is consistent with our present picture at least as long as the preparations. However, a Russian 
of the radiation processes in the ionosphere. geophysical rocket, carrying ionization and magnetic gages in 


DENSITY (GMS/M°) 
@ 
= 
= & 
| 


? 


240 


bas = OB flight to 473 km, made measurements at 260 km on both 
uber The USSR Results ascent and descent. Results: Pressure at 260 km was 
a wan a 107 mm Hg. This could be compared to the ARDC model 


We have seen no detailed results from the Russian meteoro- atmosphere in which at 260 km the pressure is 3 X 107° 
logical soundings and suspect that the problem is the same mm Hg. This suggests, as do the satellite data, that the 
as that of the American groups. The data processing takes densities at these altitudes above 200 km is higher than the 
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previous extrapolations of rocket soundings had led us to ex- 
pect. 
Despite the tentative character of the results that have 
been presented here, the various figures have been drawn 
comparing these recent high latitude measurements with the 
ARDC or Panel model atmospheres which were based largely 


Some Comparisons of Results 


on the rocket data taken in middle latitudes. That there — 


would be differences between the Arctic and the middle 
latitude data could be suspected; the winds in the atmosphere 
certainly suggest it. Thus, at high latitudes the density is 
apparently lower than the ARDC model up to about 200 km, 
above which altitude the densities are higher, as suggested 
by some of the satellite densities. A difficulty with the 
satellite measurements is that they are averages over changing 
geographic regions of the atmosphere; our recent data suggest 
considerable time and space variations in these parameters. 


Conclusions 


It would be premature to attempt to draw any far-reaching 
conclusions from the limited data now available. However, 
the number of successful firings that have been conducted 
when combined with those still to be fired within the IGY 
should yield an amount of data of sufficient quality to give a 
good understanding of the world-wide circulation patterns 
and the dynamic processes of the upper atmosphere. By the 
end of the IGY meteorological rocket soundings will have 
been made 60 deg N (Fort Churchill), 50 deg N (Russia), 
32 deg N (White Sands), 12 deg N (Guam), and from ship- 
board at 50 to 65 deg N (North Atlantic) and 60 to 70 deg S 
(Antarctic). The analysis and interpretation of the data 


present an exciting prospect. _ 
> 
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Nomenclature 
6 = boundary layer thickness & 
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distance from leading edge of flat plate to point of flame 
stabilization 


Re = Reynolds number 

7, = free stream velocity an 
p = free stream density ne 
Vg = blowout velocity 4. 
a, b = parameters in the equation Vg = a + bRe pate 
m = mass flux per unit width in boundary layer ad 
u = local velocity in boundary layer 
y = distance from flat plate _ 
@ = equivalence ratio 


Introduction 


HE stabilization of premixed flames in high velocity 
streams has received much attention in recent years. 
Nevertheless, there are a number of aspects of flame stabiliza- 
tion which deserve further study and consideration. One of 
these is the influence of the upstream boundary layer thickness 


Eprror’s Note: The Technical Notes and Technical Comments sections of Jet PROPULSION are open to short manuscripts describing new 
developments or offering comments on papers previously published. Such manuscripts are published without editorial review, usually 
within two months of the date of receipt. Requirements as to style are the same as for regular contributions (see masthead page). 
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on the performance of a flameholder. In previous work, Put- 
nam (3)? studied the importance of laminar boundary layers 
on the stability limits of premixed propane-air flames stabi- 
lized on axial rod flameholders. Putnam controlled the 
boundary layer thickness by providing suction grooves along 
the rod and varied the boundary layer thickness by a factor 
of two. He found that as long as the boundary layer and the 
main flow are nonturbulent, the thickness of the former has 
little bearing on the limits of stabilization. Gross (2) studied 
stabilization on the trailing edge of a very thin plate aligned 
parallel to a laminar flow. He determined that with a lam- 
inar boundary layer, a flame may be stabilized on a flat plate 
without separation and in the absence of turbulent eddies. 
In his analysis Gross showed that the velocity gradient may 
be used conveniently in the correlation of the data. The 
work of both Gross and Putnam was concerned with laminar 
boundary layers, and their applicability to the turbulent 
regime requires further study. 

Ziemer and Cambel (5) experimentally investigated the 
stabilization of propane-air flames in the boundary layer of a 
heated flat plate maintained at constant temperature. Their 
investigation was concerned with true boundary layer burn- 
ing in which a mechanism of recirculation of hot gases was 
not necessary for stabilization. They demonstrated that the 
point of stabilization of a flame in a laminar boundary layer 
corresponds to a position at which the local velocity in the 
boundary layer equals the flame propagation velocity. 

With a laminar boundary layer, there will be some point 
in the boundary layer at which the velocity is equal to the 
burning velocity of the combustible mixture, and thus a flame 
may be stabilized. However, the existence of such a condi- 
tion is not a sufficient condition for flame stabilization, but 
the quenching distance must be considered also. Flame 
stabilization may be expected if the burning velocity equals 
the local velocity at a distance from the solid boundary at 
least equal to the quenching distance. However, it is well 
known that for flame stabilization in highly turbulent streams 
a zone of recirculation is necessary. 

In this investigation the effect of the approach boundary 
layer on flame stabilization was studied. In particular, a 
Reynolds number range of 1.5 X 10° to 8.65 X 10° was con- 
sidered. It was noticed that decreasing the boundary layer 
thickness at the flameholder decreased the maximum blowoff 
velocity and the range of the stability curves while shifting 
them in the direction of richer mixtures. 

For a turbulent boundary layer, the thickness is given by 
Schlichting (4) by the following equation 


0.37L 
where the Reynolds number is defined as v 
VoL 


It follows from Equations [1 and 2] that the boundary layer | 
thickness increases as L4/, 


Description of Apparatus and Experiments 


The tests were conducted in a combustion tunnel of conven- 
tional design having a rectangular cross-section test chamber. 
The flame was stabilized at the edge of a sudden expansion 
in the chamber walls (1). 1n accordance with Equation [1], 
the boundary layer thickness at the point of stabilization was 
varied by changing the characteristic length L preceding the 
sudden expansion. The flameholder geometry is illustrated 
in Fig. 1. In making the boundary layer calculations, no cor- 
rections were made for the laminar boundary layer present in 
the initial stage of boundary layer development. 

The blowout data for the three approach lengths used are 
shown in Fig. 1. The curves indicate that the flameholder 


* Numbers in parentheses indicate References at end of paper. 


DECEMBER 1958 


} 
| 
x 
fo} 
Flameholder 
210}+— 
190 
>= 170 
q 
2 
j 
/ 
110 -— + — 
x 
06 0.8 1.0 12 L4 6 


Equivalence Rotio @ 


Fig. 1 Blowout velocity vs. equivalence ratio for different 
approach lengths 
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Fig. 2 Blowout velocity vs. Reynolds number for constant 
equivalence ratio 


performance is improved as the thickness of the boundary 
layer at the point of stabilization is increased. The blowout 
data were correlated in the following manner. From Fig. 1, 
a cross plot of blowout velocity vs. Reynold’s number at con- 
stant equivalence ratio was made. This family of curves is 
shown in Fig. 2. The curves of Fig. 2 can be described 
analytically as 


Since the family of curves in Fig. 2 has the equivalence 
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Fig. 3 Curves of parameters a and b/a vs. equivalence ratio 


ratio as parameter, it follows that the following functional 
expressions may be written 


Values of a and 6 were determined from Fig. 2, and analyti- 
cal expressions for a and 6 in terms of equivalence ratio were 
obtained. Fig. 3 illustrates the degree of scatter of the ex- 
perimental values of a and 6 from the analytical curves ob- 
tained. The ratio b/a is plotted in Fig. 3 for later use. 

Since the curves of Fig. 1 indicated that thickening the 
boundary layer at the point of stabilization widened the 
stability limits, it seemed reasonable to expect some degree of 
correlation between blowout velocity and mass flow rate in the 
boundary layer. Such a correlation was obtained as follows: 

The mass flow rate per unit width in the boundary layer is 
given by the expression 


Since the boundary layer at the sudden expansion is turbulent, 
a velocity distribution obeying the + law is assumed. Thus 


u y 1/7 


In integrating Equation [7], the fluid may be assumed in- 
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Fig. 4 Analytical blowout curves with experimental data 


compressible, because of the relatively low free stream veloci- 
ties preceding the flameholder. Combining Equations [1, 2 
and 9] yields 

m = [10] 


Substituting the value of uw for a stoichiometric air-propane 
mixture at 70 F gives 


m = (10.42)(3.31 X {11] 
or 
Re = 671 
Combining Equations [3 and 11] yields 
Va = c(l + dm®!4).. [12] 
where 
c = c(h) = —357G? + — 364 ....... [13] 
d = d(d) = 10-3(1.855 — 1.29) ....... . [14] 


Using Equations [12, 13 and 14], blowout curves were cal- 
culated for the three characteristic lengths used in the experi- 
ments. These curves are shown in Fig. 4 along with the ex- 
perimental points for comparison. Fig. 4 serves to define 
the maximum blowout points which were not obtainable ex- 
perimentally due to an inadequate air supply. This correla- 
tion of data indicates that for the recessed wall flameholder 
in the operating-range covered by these experiments, the 
parameters determining the blowout velocity are the equiva- 
lence ratio and the mass flow rate in the boundary layer at the 
sudden expansion. 

The length of the recirculation zone was observed by NaC! 
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injection and was found to be essentially constant above a 
certain minimum Reynolds number, about 2.5 X 105. Ata 
fixed value of @, the transport of active chemical species and 
thermal energy from the recirculation zone to the mixing zone 
must increase, up to a certain point, as the velocity of the 
main stream increases. Since the data correlation indicated 
that for a fixed value of ¢, the boundary layer mass flow rate 
| increased with an increasing blowout velocity, it appears that 
a thickened boundary layer improves flameholder perform- 
ance by increasing the time rate of recirculation of hot gases. 
This reasoning can be made more plausible by a physical argu- 
1 ment which follows. 

A necessary condition for recirculation to occur is the ex- 
istence of an adverse pressure gradient. Recirculation then 
occurs because the boundary layer fluid near the wall does not 
possess sufficient momentum to penetrate the increasing pres- 
sure field. Hence, fluid particles with relatively low momen- 
tum are decelerated, reach a zero velocity in the downstream 
direction, and are ultimately accelerated in the upstream di- 
rection. Thus as the boundary layer at the sudden expansion 
is thickened, more fluid becomes available with insufficient 
momentum to penetrate the adverse pressure gradient. This 
results in an increase in time rate of recirculation of hot gases. 
An increased rate of recirculation means an increased rate 
of energy and species flux to the mixing zone resulting in im- 
| provement of flameholder performance. 


Acknowledgments 


This paper is based on the MLS. thesis of the author. The 
author is indebted to the General Electric Company for a 
fellowship grant and to Northwestern University for financial 
aid through Faculty Research Project FR P103-54. Special 
thanks are due to Professor Ali Bulent Cambel, teacher and 
friend, for suggesting the topic and giving needed advice and 
encouragement. 


References 


1. Huellmantel, L. W., Ziemer, R. W. and Cambel, A. B., 
“Stabilization of Premixed Propane-Air Flames in Recessed 
Ducts,” Jet Proputsion, vol. 27, Jan. 1957, p. 31. 

2 Gross, R. A., “Boundary Layer Flame Stabilization,” JET 
} PROPULSION, June 1955, p. 288. 

3 Putnam, A. A., “The Effect of Boundary Layer Thickness 
on Flame Stability,” Fuel, July 1954, p. 355. 

] 4 Schlichting, H., “Boundary Layer Theory-Turbulent 
Flows,’”’ U. 8. NACA T.M. 1218, April 1949. 

e 5 Ziemer, R. W. and Cambel, A. B., “Flame Stabilization in 
the Boundary Layer of Heated Plates,’’ Jer Proputsion, vol. 
28, no. 9, Sept. 1958, pp. 592-599. 

Visibility of Orbital Points 


W. J. BERGER! and J. R. RICUPITO? . 


; Conyair-Astronautics Division of General Dynamics Corp., 
San Diego, Calif. 


| Equations are developed for computing the values of 
time at which a vehicle, moving in a known ballistic orbit, 
appears visible to line-of-sight instrumentation at any 
e specified geodetic station. Transformations are found 
is for describing the vehicular position in orbital-plane co- 
= ordinates, in geocentric fixed star coordinates, and in local 
e spherical coordinates originating at the observing station. 


Introduction 


r ODERN activity in astronautics has created a need for 
C a simplified mathematical discussion concerning visi- 
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Fig. 1 Geometric configuration for satellite orbits — 
“grt 
7 


bility of orbital points. By “visibility of an orbital point” 
we mean that there exists a mathematical line joining a fixed 
station on Earth to an instantaneous position of an orbiting 
vehicle such that the line does not cut through the Earth. 
Neglecting refraction, this characterization of visibility ap- 
plies especially to all “line-of-sight”? instrumentation sys- 
tems, e.g., optical cameras, radar, Dovap, Azusa, Minitrack, 
etc. Questions of visibility occur not only in test flights of 
missiles but also in tracking of satellite vehicles and in de- 
tection of enemy ICBM’s. 


General Theory 


The general basis of the problem may be outlined as follows: 
In Fig. 1 let CC, represent the curved surface of the Earth; 
let S be an instrumentation station on the Earth, and let the 
curve M,, M, M; denote the path of a vehicle (missile or arti- 
ficial moon) moving in a gravitational orbit about the Earth. 
The plane tangent to the Earth at S cuts the orbit at points 
M, and M,, and we suppose that the vehicle passes M, before 
M2. 

Assume a geocentric equatorial Cartesian system with 
(U,V,W) axes having constant directions relative to the 
“fixed stars.”? The W axis may be taken coincident with 
the axis of rotation of the Earth. In the (U, V, W) system 
let the orbit of M be represented by the functions 


U = alt) V = Bt) [1] 


As an approximate picture, sufficient for our purpose, we can 
consider the orbit of to be a fixed track relative to the equa- 
torial (U, V, W) system, with the Earth spinning about the 
W axis at constant angular velocity 2. Because of gravita- 
tional perturbations and other disturbing effects, such as 
collisions of the vehicle with ions and micrometeorites, the 
elements of an orbit may vary slowly with time. If such 


3’ The “fixed stars,” an old astronomical term, denotes the 
stellar system which defines sidereal time. 

‘ Because of precession and nutation, the axis of the Earth is 
not strictly constant in direction, but these effects are negligible 
in our problem. 
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variations are known functions of time, they can be readily 
included in our theory, since they enter only algebraically, 
and we do not have occasion to evaluate their time derivatives. 
Any station S on the Earth has known coordinates as func- 
tions of time. 


where p and S,, are constants determined by che figure of the 
Earth. 

When a particular mathematical model is assumed to repre- 
sent the figure of the Earth, we can specify, as a known func- 
tion of time, the unit vector Vs normal to the Earth at S. 
We can construct, as a function of time, the vector P = 
SM joining S to any orbital point M regardless of whether 
or not this vector cuts through the Earth. 

Any value of time for which M is visible from S must satisfy 
the inequality 


[3] 


If P-/ Ns <0, then the point M isnot visible from S. Values 
of time ty, and ty, such that P- Ns = 0 determine the first 
and the last visible points 1, and M2. 

For ty, < t < ty, it will be of interest to transform (a, 8, 
y) into local coordinates originating at S. 


Specialized Kinematic Theory 


To obtain a simple mathematical method, we assume that 


the orbiting vehicle describes a Keplerian ellipse with the 
centroid of the Earth as one focus. Since the gravitational 
potential function of the Earth is not spherically symmetrical, 
the elliptical orbit will be subject to dynamic perturbations 
which we shall not consider. However, we retain the sphe- 
roidal shape of the Earth insofar as it influences the geometry 
of the problem. 
PoOsITION IN THE ORBITAL PLANE In Fig. 1, curve M,, M, M2 
now represents an ellipse with one focus at the centroid O 
of the Earth. Let X, Y and Z be orthogonal Cartesian axes 
originating at O with the Z axis normal to the plane of the 
orbit and the Y axis pointing toward the perigee P;. We 
suppose the orientation to be such that, if X, Y and Z are 
unit vectors along these axes, then the rotation of the X axis 
into the Y axis gives the sense of motion of the oribiting ve- 
hicle provided that the vector product X @ Y = Z. 

In classical mechanics, six “elements” are necessary to de- 
fine an elliptical orbit. These are 


=semimajoraxisoftheellipse j | 
e = eccentricity wih 
to = time of perigee passage 


and three angles (J, L, Q) which fix the orientation of the orbit. 
Usually these angles are defined as being 


I = inclination of the orbit. We take J positive when meas- 
ured from the equatorial plane toward the orbital plane 
on the descending side of the orbit. 

L = longitude of a node (ascending or descending). In Fig. 
1 we show L as positive measured westward in the 
equatorial plane from the U axis to the line OD which 
points to the descending node D. 

Q = angle between the major axis of the orbit and a nodal 
line. We take Q positive when measured from the 
X axis to the line OD of the descending node D. 


The position of the vehicle in the orbit is described by Kep- 
ler’s equations 


a(1 — e cos E) 


+*) tan 


ll 


6 
= 


t) = time of perigee passage 

t = any later time TD Pay 

T = periodic time 

E = eccentric anomaly 
r = |OM| = radial distance of the vehicle at time ha 

= “mean anomaly’’ 
A = “true anomaly” = angle described by OM in the orbit! 


plane during the time between the most recent perigee 
passage and the present position at time ¢. 


As in the text by Smart (4),® series expressions for FE, .1 
and hence r can be developed in terms of @, namely 


3 
sino sin 20 + “sin 30 + 


sin 36 + 
[6] 


These expansions (2) were originally obtained by Lagrange. 

Laplace proved that Equation [5] converges if e< 0.662743; 
a condition which happens to hold for orbits of most ballistic 
missiles and satellites. 

Adding 2/2 to A, we obtain y, the usual polar angle of 0./ 
measured in the orbital plane counterclockwise from the 
X axis. 

Hence 


= —rsinA 
=rcosA 


gives the position of the vehicle in the (X,Y,Z) system at any 
time t. 

Position IN THE (U, V, W) System _ By vectorial geom- 
etry (3), with (X, Y, Z) and (U, V, W) as unit orthogonal 


vectors along the respective axes, we find 


= YurYw 
x. = cosL cosQ — cos/ sin L sin Q 
xz, = —sin L cosQ — cos! cos L sin Q 
= sin] sinQ 
yu = cos LsinQ + cosT sin L cos Q 
yo = —sin L sinQ@ + cosI cos L cosQ 
Ye = —sinI cosQ 
z = sinJ cosL 
20 = cos! 


5 Traditionally the mean anomaly is denoted by the symbol 1/ 
which then conflicts with the symbol for mass. 
6 Numbers in parentheses indicate References at end of paper. 
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Since Equation [8] gives the components (x, y, z) of the 
position vector R(t) of the vehicle in the (X, Y, Z) system, 
we have by Equation [9] 


R = (2,y,2)| ¥ | = (2, Yu Yo Yo || V [10] 
2u 2w W 


CRITERION OF VISIBILITY 


If 

S = position vector of station S relative to O — 

Q = angular velocity of Earth rotation & 

p = distance from the W axis to S 

@ = geodetic latitude (1) of station S 

ts = time of perigee passage i 

Xo = longitude of the station S at to Ya io 


measured in the (U, V) plane from the U axis in same sense 
(clockwise) as L, then the angle 7 between U and the pro- 
jection of S in the (U, V) plane, measured counterclockwise 
from U is given by 


Hence 
+ S.W............. [12] 
[13] 
Sw 


and Vs, the unit normal vector at S, is given by 


Ns =N.O+ NV +NW............ [14] 
N, cos cos 
N, 1 cos @ ain {15] 
Nw sin @ 
By Equations [8, 12, 14], the criterion of visibility is 
[16] 


Position RELATIVE TO THE Station At the station S, 
let (Si, Se, S3) be an orthogonal coordinate system with the 
S; axis pointing upward (along Ws), and the S, and S2 axes 
pointing east and north, respectively. 
Let (.5;, Se, 53) be unit vectors along the (S;, S2, S3) axes. 
We have already defined 5, = Ns by Equations [14 and 
15]. Similarly 


5, = (—sin [17] 
= —(sin cos n)U — (sin @sin n)V + (cos @)W..... [18] 


Since MZ = R — Sis the position of the vehicle V/ relative 
to S, then the coordinates of M in the (S,, So, S3) system are 


(M-S,, M-5:, M-S3)...... [19] 


Other useful coordinates relative to S are the azimuth ag, 
elevation es, and range Rs of M given by 


M-S, 
M-S: 
€s = are tan [21] 
— 
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Calculation of Re-Entry Velocity Profile 


ALFRED A. ADLER! 


Bell Aircraft Corp., Buffalo, N. Y. 


A simple approximate equation, widely used to compute 
the velocity vs. altitude profile of a ballistic missile re- 
entering the Earth’s atmosphere, gives quite a good ap- 
proximation of the results sought. Modifications to the 
derivation of this equation are presented which include 
the effect of gravity and provide for a more accurate repre- 
sentation of the air density variation. These refinements, 
to a large extent, eliminate the deleterious effects of the 
assumptions made in the more simple equation without 
unduly increasing the complexity of the final relation. 
Two equations result, one covering the altitude range from 
initial re-entry down to 75,000 ft, the other providing for 
the remaining distance to sea level. 


Nomenclature 
V = velocity, ft/sec 
h = altitude, ft 
p = air density, lb sec?/ft* oO 
B = 1 In Ps 
hp 

Cp = vehicle drag coefficient ut, 
S = vehicle reference area, ft? ; 
W = vehicle weight, lb 
g = gravitational acceleration, ft/sec? 
y = flight path angle with respect to local horizonta. 
r = geocentric radius, ft 
= gsrs?, 14058.65 X 101? ft?/sec? 

72 
E’ = total energy per unit mass, = = ft?/sec? 

2 r 
Subscripts 
0 = initial condition 
s = sea level 
k = assumed constant value 


N COMPUTING the velocity vs. altitude profile of a 
ballistic vehicle re-entering the Earth’s atmosphere, the 
following equation is widely used 


where 
26; W sin ly| 
B = B = 41.5 X 107° 


This equation results from solving the general equation of 
motion in the velocity direction, assuming that the effect of 
gravity is negligibly small compared to the drag effect, and 
assuming that @ and a are constant. These assumptions are 
necessary to obtain a very simple solution. 

Although Equation [1] yields a good approximation to the 
results sought for certain uses, a comparison of these predic- 
tions with a machine-computed numerical solution to the 
equation of motion assuming only that @ remains constant, 
indicates that the assumptions of zero gravity and constant 6 
made in the derivation of Equation [1] are to be avoided if 
possible. Due to the lack of a gravity term, the velocity pre- 
dicted by Equation [1] decreases slowly with decreasing alti- 
tude for large values of altitude, whereas actually the velocity 
increases with decreasing altitude at high altitudes. Thus, 


~ Received June 10, 1958. 
1 Aeromechanics Section, Space Flight Division. Member ARS. 
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the velocity vs. altitude profile given by Equation [1] for a 


particular problem will depend on the assumed initial re-entry 


altitude. The assumption that 6 remains constant, though 
a simplifying assumption, is at variance with published data 
(1).2 It will be shown that it is possible to make modifica- 
tions to the derivation of Equation [1] which, to a large ex- 
tent, eliminate the deleterious effects of these two assumptions 
without increasing the complexity of the final relation unduly. 

It is reasonable to state that the velocity of a re-entering 
ballistic vehicle at any point in its descent is the sum of the 
initial re-entry velocity, the velocity gained since initial re- 
entry due to gravity, and the velocity lost since initial re- 
entry due to drag. The velocity gained due to gravity can b 
determined by considering the interchange between potentia! 
and kinetic energy in an inverse square central force field 
The velocity lost due to drag is essentially given by Equatioi: 
[i]. Thus 


= 
= 2 
rath 


This relation amounts to an analytical solution of the equatio: 
of motion assuming 6 and a are constant and that the gravit; 
and drag effects are linearly independent. The addition o 
the gravity dependence increases the accuracy of the solution 
considerably, and makes the final result no longer dependent 
on the choice of initial re-entry altitude. 

A plot of 6 vs. altitude from ARDC 1956 data (see Fig. 1 
indicates that while 8 = const is a good approximation abov: 
75,000 ft, below this altitude 6 is more nearly a linear function 
of altitude. If Equation [1] is rederived, holding a constant 
and letting 8 = 6, — m (75,000 — h) rather than a constant, 
and the results substituted in Equation [2] in place of the drag 
term, the relation to be used below 75,000 ft becomes 


— — 75,0008 


exp (see ¥ ve lerf (W/mh + Vb) — erf (75,000 


Vm + vo) 


From Fig. 1, = 13.416 X Vb = B,/2/m = 
1.0435. Substituting these values 


V= v2 + h + Vo x 


ag 
exp (8.1446a[erf(13.416 X 10-*h + 1.0435) —0.99626}) 


Thus, the velocity vs. altitude profile calculation should be 
made using Equation [2] over the altitude range from initial 
re-entry down to 75,000 ft, and using Equation [3] from 75,000 
ft to sea level. It is apparent that these equations may read- 
ily be adapted to the use of a piecewise constant a to allow 
for variable Cp>S/W sin 

A comparison of the differences between the velocity vs. 
altitude profiles computed by use of Equation [1] and by use 
of Equations [2 and 3] and the results of a machine-computed 
numerical solution of the equation of motion, holding a con- 


stant, is presented in Fig. 2 for two initial re-entry condi- 


tions. The results indicate that where greater accuracy is 
required without undue effort, Equations [2 and 3] may be 
used. 


2 Minzer, R. A. and Ripley, W. 8., “The ARDC Model 
—AFORC TN 1956,” Air Force Surveys in Geophysics no. 86, 
AFCRC TN-56-024, Dec. 1956. 
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S. L. BRAGG! and S. P. Q. BYWORTH? 
Rolls-Royce Ltd., Derby, England 


TP HE characteristic velocity and thrust coefficient values 
calculated for an ideal thrust chamber are based on the 
total pressure at the throat, P;. The normal test bed meas- 
urement is that of static pressure at the injector plate, p-. 
The following analysis demonstrates a simple method of 
calculating the ratio between the two pressures. 
The equation of motion of flow in the combustion chamber, 
which is assumed to be parallel, and of area Ag, is 


= MV 
where 

Z = axial momentum of the liquid propellants at entry | 

M = total propellant flow 

p = static pressure at combustion chamber exit 

V = velocity at combination chamber exit 

The actual characteristic velocity is defined by the relation 
: AP 
[2 


where A, is the throat area. 
Substituting for 17 in Equation [1] and rearranging 


Z AW 

Pe AcDe Acc Ps 
Now the dimensionless groups V/c* and p/p: are functions 

of the area ratio A,/A, and the physical characteristics of the 


exhausting gases. The figure shows the variation of the de- 
nominator function 


s 


(4242)" 


A,.c* 


with area ratio, for extreme values of the specific heat ratio- 
y = 1.0 and y = 1.4. Also shown are two curves for the ex- 
pansion of a typical mixture (2.35:1, oxygen: kerosene), assum- 
ing shifting (equilibrium) composition, and constant (frozen) 
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transfer may be defined as 


composition, respectively. The latter corresponds to y = 1.23 
approximately. 

It is apparent that the value of y, or the assumptions made 
as to composition changes during expansion, have compara- 
tively little effect on the relation between the demoninator 
function and the area ratio. Thus the relation is probably also 
unaffected by combustion inefficiency, or any other departure 
from ideal performance in the combustion chamber, and the 
mean curve can safely be used for practical analysis. 

The term Z/A-p¢ is usually small but not negligible com- 
pared to unity. If each propellant is injected under an injec- 
tion pressure differential AP, through ports of effective area 
A,, inclined at @ to the chamber axis 


[4] 


This equation, used in conjunction with Equation [3] and 
the graph, enables the throat total pressure at any operating 
condition to be calculated from the measured injector end 
static pressure, in a very simple way. 


An Empirical Method for Calculating 
Heat Transfer Rates in Resonating 
Gaseous Pipe Flow 


GERALD MORRELL! 


National Aeronautics and Space Administration, 
Cleveland, Ohio 


UBIANSKII (1)? has shown experimentally that a sound 

field at resonance can produce an increase in heat transfer 
rate from a heated surface to surrounding air provided the 
velocity perturbations are at least as large as the steady 
velocity. Male, Kerslake and Tischler (2) have reported 
large increases in wall erosion rates (hence, heat transfer 
rates) accompanying large amplitude resonance burning in 
cylindrical rocket combustors. Since exact solutions of the 
energy and flow equations even for simplified cases of linear 
perturbations are tedious and difficult (3), it appeared that an 
approximate method suitable for many engineering applica- 
tions would be worthwhile. 

In this presentation the following major assumptions have 
been made: 

1 The perturbations are sufficiently similar to shocks that 
the one-dimensional pressure-velocity relations for normal 
shock waves can be used. 

2 The usual Stanton number-Reynolds number-Prandtl 
number relation applies to both resonance and steady flow 
provided suitable average conditions behind the shock are 
selected. 

3 The effect of aerodynamic heating due to the oscillations 
can be treated as a correction factor to the quasi-steady state 
solution. 

The ratio of resonance heat transfer to steady-state heat 


where T' is absolute temperature, and subscripts s, 0 and w 


_ refer to average conditions behind the shock, steady condi- 


tions, and wall conditions, respectively. The latter condition 
jis assumed the same for both cases. The film coefficient h is 


given by 
DV cpt —0.6 
h 
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HEAT TRANSFER RATIO PER CYCLE, Ro ve 


PRESSURE RATIO, € 


Fig. 1 Typical curves of heat transfer ratio as a function of 
pressure ratio for several values of Mach number /, and tem- 
perature ratio 0 


where the symbols have their usual meaning. Assuming that 
the ideal gas law applies, that (u/k), molecular weight and 
heat capacity are invariant, and that up « 7, the heat transfer 
ratio becomes 


where p is pressure. 

The ratio of maximum pressure behind the shock to steady- 
state pressure is assumed to be nearly equal to the shock 
pressure ratio 


and the shock temperature ratio 7,/T; as given in (4) is as- 
sumed to be nearly equal to the ratio of maximum tempera- 


ture behind the shock and steady-state temperature me 


a) 
ag 1 T» 

where a = (y + 1)/(y — 1) and ¥ is the specific heat ratio. 

For one-dimensional flow, the particle velocity behind the 

shock (5) is given by 

(a — 1)(E — 

V(a@ + 1)(1 + @é) 

where @ is the speed of sound at temperature 7», the plus sign 

is for a wave moving in the same direction as the undisturbed 

velocity Vo, and the minus sign is for a wave moving in the 

opposite direction.- For a reflected wave, the heat transfer 

ratio per cycle will be given by the average of these two con- 

ditions. 


830 


Vo+t 


> 


For these calculations average conditions behind the shock 
were chosen 


2, = 0.52 + 1) 
Zo 


where x is velocity, pressure or temperature. The expres- 


sions for heat transfer ratio become 
0.5/4, + + 1)°*10.5(B + 1) 8] 
(B + — @) 


_ + 1] + 1)°*10.5(B + 1) — 8) 
(B + 1)°“(1 — 8) 


where the plus subscript is for the forward facing wave and 
the minus subscript for the reflected wave, and 


R- 


(a — 1)(€ 1) 
MoV + 1)(1 + @&) 
MoV(a + 1)(1 + 


and M, is the Mach number of the undisturbed flow. 

The correction for aerodynamic heating due to gas oscilla- 
tions is assumed to be given by Ostrach’s derivation (6) for an 
oscillating flat plate. The expression, in terms of a heat 
transfer ratio, may be written 


— VPr) 


and 


At + 2|A 1) 
16,m(1—8) 
1 
=1+ — 8 (A2 + +1 


where Pr is Prandtl number, R, is the universal gas constant 
and m is gas molecular weight. 


The overall heat transfer ratios are 
Ro+ = Rs R: 
and 
Ry- = R-R! 


and for a reflecting wave, the heat transfer ratio per cycle is 
given by 


+ 


Ro, per cycle = a 


For small values of Mo and large values of &, |A| for the two 
cases will be nearly the same, and it becomes unnecessary to 
perform the averaging process in view of the other approxima- 
tions involved. Fig. 1 is a plot of Ry per cycle as a function of 
for several values of and Mg, and for c, = 0.3 Btu/lb deg R, 
y = 1.25 (a = 9), m = 22, and Pr = 0.75. The latter values 
are typical of some rocket combustion gases. 

If the erosion rates of (2) are assumed to be proportional to 
heat transfer rates, these data may be used to check the 
calculation method. From the engine performance data the 
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values of My and T>) are calculated to be 0.125 and 3650 R, re- 
spectively. The values of cp, a, m and Pr may be assumed 
to be the same as those above. The softening point of the 
plastic used for a wall material is given in (7) as 713 deg R, so 
that @ is approximately 0.2. For longitudinal waves the only 
measured value of £ was about 2, and the experimental value 
of Ry was about 3; the calculated value of Ro is 3.37. For the 
traveling wave form of the transverse oscillations £ is about 
3.7, and the maximum experimental value of Ry is about 11. 
Assuming the same values for My and 6 (this neglects the dif- 
ference in direction for Vo and V,) and considering only the 
forward traveling wave, the calculated value of Rp is 9.86. 
This agreement indicates that the calculation method pre- 
sented can give useful engineering estimates of the heat 
transfer rates produced in a resonating flow. 
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2 
The Gasification of Solid Ammonium 

Nitrate! 


W. H. ANDERSEN,’ K. W. BILLS,’ A. O. DEKKER,‘ 
E. MISHUCK,*® G. MOE’ and R. D. SCHULTZ’ 


Aerojet-General Corp., Azusa, Calif. 


When the surface of solid ammonium nitrate is super- 
heated by a high flux heating technique the ammonium 
nitrate decomposes in a different manner than the bulk 
phase liquid. 


Introduction 


N ORDER to understand the burning of ammonium ni- 

trate (AN) composite propellants, information is required 
about the rates and the nature of the chemiéal reactions oc- 
curring in the surface layer of initially solid oxidizer which 
is exposed to the flame zone. A number of investigations 
have been directed toward an elucidation of the mechanism 
of the thermal decomposition of AN (1).8 However, in all 
instances, the studies were carried out with the molten salt. 
There is general agreement that N.O and H,0 are the chief 
products of the decomposition of the liquid, from 169.6 C 
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(melting point) to about 260 C, and at atmospheric pres- 
sure. At higher temperatures and pressures, side reactions 
occur in the gas phase between the small amounts of the dis- 
sociation products, HNO; and NH;. The extent to which 
these side reactions occur increases with the temperature of 
the liquid and with the ratio of the exposed surface to the 
volume of the melt. Kinetically, the liquid phase (bulk) 
decomposition is a first order reaction, having an activation 
energy of about 38 kcal/mole (2). Ther saction is exothermic 
to the extent of about 9 kcal/mole. 

This note describes some experiments in which the surface 
of solid AN is heated at such a high rate of heat input that 
the surface gasifies at a rate which is much greater than the 
bulk decomposition rate. Under these conditions the surface 
gasification reaction proceeds by a mechanism which involves 
a small activation energy (about 7 kcal/mole), is highly 
endothermic, and produces essentially only — vapor and 
NH; as the primary products. This surface reaction is the 
direct result of the fast rate of heating of the surface, and 
appears to involve the same phenomena of solid superheating 
reported by Schultz and Dekker for ammonium chloride (2). 


Experiment 


The linear gasification rate of solid AN was measured as a 
function of surface temperature (7',) with an improved ver- 
sion of the “hot plate” device described by Schultz and Dek- 
ker (2). This method has the important advantage of close 
simulation of conditions that exist during propellant 
burning. Briefly, the measurements were made by pressing 
strands of AN against a heated surface, and measuring the 
linear rates of regression of the surface of the samples xs a 
function of 7,. 7’, was measured by means of a thin, flattencd 
thermocouple which was interpos: d between the s ample sur- 
face and a thin sheet of mica covering the hot element. The 
apparatus used was a prototype of the hot plate pyrolysis 
apparatus which is described in detail elsewhere (3). The 
samples were pressed vertically against the surface of the 
heating element with a constant loading pressure of about 
750 gm/cm?. The measurements were made in an atmos- 
phere of nitrogen at pressures of 3, 14.7 and 600 psia. In the 
experiments, the samples gasified cleanly; little evidence for 
the accumulation of liquid or solid residues was observed. 

The AN strands were prepared for pyrolysis by pressing 
dry, cp AN powder (particle size about 10 u) in a steel die 
to a pressure of 17,000 psi. Density measurements (oil 
immersion method) indicated that average strand densities 
of about 1.69 gm/cm? were obtained. The strands thus were 
close to crystal density, 1.72 gm/cm*. 


Results and Discussion 


The experimental linear pyrolysis rates B for strands of 
AN are shown in Fig. 1. The dependence of the pyrolysis 
rate on surface temperature in the range 180 to 300 C can be 
represented to good approximation by the equation 


B = 120 exp( —7100/RT3;) em/sec......... {1] 


where 7’, is the surface temperature, K, and FR is the gas con- 
stant. With the available apparatus, the activation energy 
is believed to be reliable to only about +3 keal. However, 
this fact does not affect the main point of this note, namely, 
that the surface gasification of AN is entirely different from 
the bulk dee omposition reaction. 

The linear pyrolysis of AN is a highly endothermic process, 
as is evidenced by the fact that a very large amount of heat is 
removed from the hot plate during the pyrolysis of the 
sample. This, together with analyses of the pyrolysis prod- 
ucts,? indicated that the overall process occurring in a very 

* Semiquantitative infrared and mass spectrometric analyses 
show that the major pyrolysis product is sublimed ammonium 


nitrate; very small amounts of NO», NO, Es , H.O, NH; and per- 
haps O» are also present. Below 200 C, NO was also found. 
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Fig. 1 Linear pyrolysis rate of ammonium nitrate as a function 
of surface temperature 


thin layer of AN exposed to a heat source, such as a hot plate 
or a flame, is chiefly the dissociation reaction 


NH,NO; (s) > NH; (g) + HNO; (g).......-- [2] 


which has an endothermicity of about 40 kcal/mole. Sur- 
prisingly, the activation energy for the linear pyrolysis of 
AN is only about 7 kcal/mole, indicating that a stepwise 
surface process is probably involved. The situation thus 
appears to be similar to that in NH,Cl, where the dissociation 
reaction is likewise highly endothermic, whereas the rate of 
reaction is controlled by a low activation energy (2). 

Fig. 1 shows that below about 180 C, the pyrolysis rates 
drop sharply, and become essentially independent of tem- 
perature. Moreover, the rates appear independent of pres- 
sure in the range studied. These results imply that the sharp 
decrease in pyrolysis rate below about 180 C is related to a 
melting process rather than a sublimation process, as in the 
case of NH,Cl (2). 

The mechanism of the surface dissociation of AN is not 
known, but possibly involves the following sequence of steps 


NH,NO; (s) [3] 
NHNO; (liq) == NH; (ads) + HNO; (ads)....... [4] 
[5] 
HNO; (ads) = HNO, [6] 


where NH; (ads), HNO; (ads) represent molecules physically 
adsorbed on a microscopically thin (solid or liquid) AN sur- 
face laver. The heat of vaporization of pure liquid HNO; is 
of the same magnitude as the activation energy for pyrolysis, 
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and hence it is reasonable to assume that during the linear 
pyrolysis of AN the desorption of HNOs, reaction [6], may be 
the rate-controlling step. The heat of fusion of AN is 1460 
cal/mole, and it is thus unlikely that reaction [3] is rate 
controlling. 

The thermal decomposition of any material at a given 
temperature proceeds by the fastest mechanism compatible 
with the conditions. The preceding results and discussion 
indicate that the rate of the solid surface decomposition of 
AN is greater than that for the bulk decomposition. This 
has important consequences for the interpretation of com- 
bustion phenomena involving the surface thermal decom- 
position of AN. In particular, it has led to new models for 
interpretation of the burning (4, 5) and detonation (6) char- 
acteristics of AN. 

In conclusion it is suggested that, according to the observa- 
tions described above, the rate-controlling process in the 
linear pyrolysis of solid AN occurs in the surface layer only, 
and not below the surface. The surface reaction is highly 
endothermic. The possibility of partial decomposition in 
the subsurface layers is not ruled out, but this partial de- 
composition is not considered to be rate-controlling. 
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Pressure Distributions on Blunt-Nosed 
Cones in Low Density Hypersonic Flow' 


L. TALBOT,? S. A. SCHAAF? and F. C. HURLBUT‘ 


University of California, Berkeley, Calif. 


> Ne note describes an experimental investigation of sur- 
face pressures on a class of blunt bodies in low density 
hypersonic flow. The general object of this study was to in- 
vestigate the effect of viscous interaction on the surface 
pressure. 

The configuration which was studied was a family of cones 
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Fig. 1 Surface pressure distribution, 40 deg cone angle 


with spherically-blunted noses. These cones were similar to 
those investigated by Machell and O’Bryant (1).5 The latter 
experiments were performed at flow conditions M, ~ 5.8 
and 6 X 104 < Rerw < 1.6 X 105, where M,, is the free 
stream Mach number, and Rere is the Reynolds number 
based on free stream conditions and the nose radius. The 
present tests covered the range 3.6 < M, < 5.8, but were 
at considerably lower Reynolds numbers, in the range 2 X 
10? < Rerw < 5 X 10°. At these lower Reynolds numbers, 
some viscous interaction effects can be expected to be pres- 
ent, and comparisons between the present tests and the re- 
sults of (1) thus serve to indicate the magnitude of such 
effects. 

The model dimensions are shown in Figs. 1 to 5. All 
models were of base diameter 1.75 in. The data points indi- 
cate the location of the pressure orifices, in terms of x/R, 
with x measured from the stagnation point along the surface 
of the model. Most of the tests were conducted in the 
nominal Mach 6 nozzle, under two flow conditions: M, = 
5.8, Re,,/in. 8980, and M,, & 5.5, Re,/in. 4500. A few 
additional results were obtained in the nominal Mach 4 
nozzle. Typical flow conditions are listed in Table 1. Ac- 
curacy of the present experiments is estimated to be within 
+2 per cent in pressure, +5 per cent in free stream Reynolds 
number and +1 per cent in Mach number. The values of 
x/R are accurate to +2 per cent on all but the models with 
0.05 in. nose radius. For these models, the orifices on the 
spherical portions could be located to about +10 per cent 
accuracy in 

It is convenient to discuss the results in terms of two re- 
gions on the model surface, the stagnation region, and the 
“transition” region, which is the region of the surface in the 
vicinity of the change in curvature between the spherical 
nose and the conical afterbody. 


Stagnation Region 


Figs. 1 to 5 show that in the stagnation region the surface 
pressures follow quite closely the “modified-Newtonian” 
approximation 


(1] 


where C, = 2(p — p,.)/p.U..” is the local pressure coefficient, 
Cymax the value of the pressure coefficient at the stagnation 
point, and @ the angle between the local normal to the surface 
and the direction of the undisturbed free stream velocity 
U.,. This agreement with the modified-Newtonian values 
was also found by Machell and O’Bryant. The fact that 
modified-Newtonian theory predicts the surface pressures in 
the stagnation region even at the Reynolds numbers of the 
present tests is really not too surprising. It has been pointed 
out by several authors (see (2) for a review of the subject and 
a comprehensive bibliography) that the success of the 
modified-Newtonian result can be “explained” by the fact 


5 Numbers in parentheses indicate References at end of paper. 
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that the decrease in surface pressure, caused by centrifugal 
effects (which yields the Newtonian-plus-centrifugal formula 
C, = 2 cos? 6 — (2/8) sin? 6), is almost exactly balanced by 
the increase in surface pressure due to the fact that the shock 
wave curvature is less than the body curvature. Now, the 
pressure drop through the shock layer Ap=> (1/3) (p,,. U,,?sin? 8) 
is independent of the shock layer thickness 6, to first order 
in k = p/p, where subscript 2 refers to conditions 
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behind the normal shock. The same appears to be true for 
the counterbalancing pressure increment due to lack of con- 
centricity of shock and body, because the modified-Newtonian 
theory is accurate over a wide range of 6/R. 

The effect of the boundary layer in the stagnation region 
is to change the effective thickness of the shock layer by the 
é*, the boundary layer displacement thickness. Detailed 
calculations show that in the stagnation region, because of the 
strong favorable pressure gradient, 6* grows very slowly with 
x and is in fact practically constant. Hence, one might expect 
the modified-Newtonian result to hold accurately in hyper- 
sonic flow even when 6* is not negligible, perhaps even up 
to the point where 6* approaches the inviscid shock layer 
thickness 6, and the two layers merge. Now, under wind 
tunnel conditions at M =~ 6, 6*/—R ~1/¥V Rer., (8); and to 
first approximation the shock layer thickness is given by 
6/R ~k (2). Thus the two layers can be expected to be dis- 
tinct provided 6*/6 ~ 1/k-/ Rep.. is somewhat less than 
unity. In the present tests, 6*/6 < 0.4. 

The actual values of Cp... found in the tests agreed, within 
experimental error, with their theoretical inviscid flow values 
given by the Rayleigh pitot tube formulas. Viscous effects 
on impact pressures (3) begin to be important (for round- 
nosed bodies) for Reza ~ 100 and less, which is a little lower 
than the range of the present data. Noncontinuum effects 
(i.e., slip flow phenomena) were likewise totally indiscernible 
in the data. It is now well established that such effects, if 
present, are overshadowed by viscous interaction phenomena 
at a “moderately low” Reynolds number (perhaps even for 
Reynolds numbers as low as Rerw ~ 10). In view of the fact 
that in the experiments the viscous interaction effects them- 
selves were unimportant, at least in the stagnation region, the 
absence of noncontinuum phenomena is only to be expected. 


Transition Region 


In this region of the body surface, the pressure distribution 
departs from modified-Newtonian to approach its final down- 
stream value on the conical afterbody. The ideal inviscid 
flow value on the cone is the Taylor-Maccoll value, as given 
in the Kopal Tables (4). The Kopal value is greater than 
the modified-Newtonian value for the cone surface, although 
for small cone angles the two are close together. 

The inviscid pressure distributions for blunt cones with 
large cone angle, as found by Machell and O’Bryant, are in- 
dicated in the figures. In the vicinity of the change in 
curvature, the flow overexpands, and the pressure drops to 
nearly the modified-Newtonian value for the cone. A re- 
compression region follows, and the final pressure is the Kopal 
value. The amount of overexpansion becomes progressively 
less as the cone angle is decreased. (Of course, the Kopal 
and modified-Newtonian values are nearly equal for small 
cone angles.) For blunt-nosed slender bodies (5) the over- 
expansion vanishes altogether. 

The present data, in comparison with those of (1), exhibit 
less overexpansion in the transition region. Also, the lengths 
of the transition regions on the models, measured in units of 
x/R, are somewhat greater in the present tests. It seems 
probable that these differences are caused by viscous inter- 
action. One may note that, at the lower Reynolds numbers 
of the M., ~ 4 data, there is no overexpansion at all evident 
on Model 1. 

Far back on the cone, downstream of the transition zone, 
one expects that due to weak interaction viscous effects the 
surface pressure will be slightly higher than the inviscid flow 
value, as is observed on slender pointed cones (6, 7) and 
slender blunt-nosed bodies (2, 5). The present tests were 
not designed to provide the high accuracy of pressure meas- 
urement in the weak interaction region necessary for a careful 
study of this region, and, in any case, the shortness of the 
models would have rendered any such measurements highly 
suspect. But within the accuracy of measurement, the pres- 


sures downstream of the transition zone were found to be 
slightly higher than the ideal Kopal values, the differences 
being of magnitudes commensurate with the predictions of 
the weak interaction theories (7). Of course, for large angle 
cones, the viscous interaction pressure increase in the weak 
interaction region is a very small fraction of the Kopal pres- 
sure, and is not even noticeable on the scale to which the 
data have been plotted. 


= 
‘Table 1 Flow conditions of tests 
Peo Measured 
(u Hg) Comax Re. /in 
3.69 57.2 1.80 1020. | 
4.05 111 2660 
5.46 4500 
5.80 107 8980 
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Note on Interplanetary Navigation 
ROBERT M. L. BAKER Jr.' 


Aeronutronic Systems, Inc., Newport Beach, Calif. 


ECAUSE mass onboard the interplanetary vehicle would 

be so expensive and reliability so important, techniques 
should be found for minimizing the amount of in-flight ob- 
servation and computation required by making the maximum 
use of terrestrial precomputation of trajectories. Thus, ob- 
servations would, for example, take the form of comparing 
the predicted and observed positions of all observable solar 
system bodies against the fixed star background, while the 
onboard computations would deal only with the differences 
between the precomputed and the actual trajectories. Pre- 
liminary studies have shown that these procedures would re- 
sult in a substantial savings in equipment weight and com- 
plexity, and that such a self-contained navigational unit would 
probably be superior to the employment of Earth-bound 
computers and a data link. 

Basically, it would seem that the navigational require- 

ments for an interplanetary voyage would be twofold: First, 
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to determine the true orbit that the space vehicle is proceeding 

along, and, second, to compute the corrective thrusts needed 
to place the vehicle on an orbit that will accomplish the de- 
sired mission (implicit in the latter is a knowledge of the ve- 
hicle’s orientation); see (1, 2 and 3).2. These requirements 
need not, however, be carried out explicitly. In fact, it will 
probably never be necessary to determine exactly where the 
vehicle is, since the real question is how to get it to its objec- 
tive. By eliminating the middleman (i.e., by eliminating 
the determination of the vehicle’s actual orbit), a great sav- 
ings can be achieved in computer complexity and computa- 
tional time. Under the assumption that the vehicle must de- 
termine its own corrective maneuver, the navigator or navi- 
gational device will have certain raw data available, such as 
star and planetary positions (obtained by celestial observa- 
tions or by star-trackers) and some electronic information 
(perhaps radar near the Earth or the target planet and a 
LORAN type signal in midcourse), all information of com- 
parable accuracy being useful. Such rough data must be re- 
duced to terms from which a corrective thrust program can 
be calculated. Precomputation would permit these calcula- 
tions to be carried out mainly before the flight. 

In general, observational inputs (probably obtained at 
different times throughout a month or so of flight), which are 
denoted by 0, Oz, O3, . . . O;, must be processed by the on- 
board computer. If the orbit is a ballistic one, then the 
corrective thrust is specified by a, 8, y, T and t,, where a, 8, 
and ¥ are angles that define the thrust direction with respect 
to some reference system, 7 is the magnitude of the corrective 
thrust, and ¢, is the time when the corrective thrust should be 
applied. These five variables are, of course, implicitly re- 
lated to the O;’s. If the orbit is a low thrust one, then the 
thrust program modification, i.e., a(t), B(t), y(t) and T(t), 
must be specified as a function of the time ¢, and these vari- 
ables will likewise be implicitly related to the O,’s. 

Precomputation, based upon accurate astronomical con- 
stants and accurate stellar and planetary position data would 
provide a set of the values of the observational quantities 
O;, that would be observed from the space vehicle if it were 
on the ideal computed orbit (i.e., a “prediction ephemeris’’). 
Of course, the space vehicle will not actually move along 
the ideal trajectory because of various unpredictable errors 
in the thrust program, uncertainties in certain astronomical 
conversion factors (which, for example, reflect themselves 
in incorrect takeoff initial conditions), unaccounted for 
perturbations, malfunctions, etc., so that the observed quan- 
tities will differ from the computed values by small differ- 
ences, denoted by Aye = O; (actually observed) — 0; 
(computed) (4). The accurate precomputation could, how- 
corrective thrusts directly from these differences, e.g. 


Qa anA,-(0;) + Q12A,-(O2) 


a — £8,y, T and ¢,. Consequently, the actual orbit need 
never be explicitly determined. 

This Taylor-series or other slightly more complex represen- 
tation can be expanded toany number of terms, as required 
by the expected vehicle inaccuracy and by the nature of the 
mission. The evaluation of these series would then be carried 
out by as'mplified onboard computer. This can be designed 
to employ a large recirculating magnetic tape as its main 
memory with, perhaps, a second recirculating tape of smaller 
size which would contain the arithmetic registers and the 
logical structure of the computer. Such a computer would 
be exactly tailored to the requirements for processing ob- 
servational data on the basis of precomputed information. 
This device can be made with a minimum of flip-flops, thus 
increasing its reliability. (It should be noted in passing 
that the importance of reliability cannot be overrated.) 
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Preliminary investigations indicate that our lack of precise 
knowledge concerning the conversion of certain astronomical 
constants might make a completely predetermined trajectory 
inadequate for all except the most rudimentary of inter- 
planetary missions (5). Such a circumstance, however, 
makes the employment of a differential correction procedure 
even more attractive, since differential corrections could 
be simply interpreted as a form of “homing,” in which the 
differentials A,, indicate the aiming error that must be elimi- 
nated by corrective thrusts; e.g., after converting from a geo- 
centric to a heliocentric reference system in midcourse, much 
of the constant error can be automatically eliminated since 
the observations will be made in the heliocentric system. 
Even if such a homing type of navigational scheme is found 
to be useful for certain missions, the requirements for pre- 
cision orbits cannot be relaxed. In the case of such terminal 
guidance, one would have to consider the trade-off between 
better observational information obtained nearer to the ob- 
jective planet as opposed to less corrective thrust required 
further away. Essentially, precision orbital techniques, 
involving precomputed data, would provide for the most 
efficient reduction of the earlier, less refined observational 
information and hence would provide for a more economical 
employment o’ corrective thrusts. 

Of course, a midcourse correction will always be required, 
no matter how accurately the orbit was established, in order 
to account for the inclination between the Earth’s and the 
planet’s orbit plane. (This inclination effect would result in 
maximum errors of 3,960,000 miles for Venus and 4,560,000 
miles for Mars.) 

It has been found (6) that (unless Venus at the time of 
landing were passing through the line of nodes) the optimum 
manner in which to transfer from the Earth’s heliocentric orbit 
plane to that of Venus or Mars would be to carry out the trans- 
fer in midcourse (in fact, almost precisely in midcourse). It 
is evident that if such a maneuver were carried out at take- 
off, it would be impossible to design a doubly tangent Earth- 
planet orbit without passing over the celestial pole—a 
maneuver that would be extremely costly in view of the fact 
that no advantage could be taken of the Earth’s speed; in 
fact this speed would have to be cancelled out. Other than 
double tangent heliocentric orbits are found to be rather 
inefficient, and even if this inefficiency were discounted, the 
penalty paid to transfer into any Earth-sun-planet orbit 
initially can become exorbitant. Spherical trigonometry can 
be employed to demonstrate that the minimum inclination 
angle between the Earth-sun-rocket plane and the planet- 
sun-rocket plane* would occur when the radius vector from 
the sun to the rocket and from the sun to the predicted posi- 
tion of the planet at landing is 90 deg. Thus, excluding other 
effects, the inclination maneuver is made most efficiently at 
this point (and not at the line of nodes of the Earth-planet 
plane; see (1)). 

The author acknowledges the contributions made to this 
note through discussions held with Drs. Samuel Herrick, Eric 
Durand and L. G. Walters cf Aeronutronic Systems, Inc. 
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Ignition by Flow Over Hot Surfaces 
WELBY G. COURTNEY! 


Experiment Incorporated, Richmond, Va. 


LTHOUGH it has long been known that a combustible 

material could be ignited by a hot surface, an interpre- 
tation of ignition phenomena in terms of heat transfer and 
chemical reactions has been severely limited by turbulence 
(1).2 Ignition during laminar flow over a heated flat surface 
(or in the laminar boundary layer when the free stream is 
turbulent) was treated theoretically by Chambré (2) and 
Toong (3). However, the present experimental work indi- 
cates that small-scale turbulence at the hot surface may be 
difficult to avoid in ignition during flow over hot surfaces. 
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Experimental Procedure 


Premixed propane-air gas at ambient temperature and 
pressure was passed at a constant initial linear flow velocity 
(volume feed rate/tube area) up an electrically heated ni- 
chrome wire positioned vertically in a glass tube (Fig. 1). 
When a minimum voltage was applied to the wire, flame 
formed at the top of the wire and then rapidly moved down 
the tube until it stabilized on the wire at a point which de- 
pended upon the initial flow velocity, burning velocity of the 
combustible mixture, flame-holding behavior of the wire and 
the voltage’ on the wire. 

The minimum total rate of energy input required for igni- 
tion was calculated from the minimum voltage for ignition 
and the wire resistance and length. The minimum rate de- 
pended critically upon the placement of the wire and tube 
with the lowest minimum rate being measured when thx 
tube and wire were vertical and the wire was located in th 
center part of the tube. The ignition data reported below 
corresponded to the minimum conditions obtained in th 
present work. 


Results and Discussion 


Shadowgraph examination of the gas as it exited from th« 
tube indicated that the gas near the wire was in laminar flow 
when the wire was cold, became turbulent when the wire was 
warmed, and was highly turbulent when the wire was heated 
to a temperature sufficient to give ignition. Although the 
free stream was always laminar, turbulence always formed 
near the wire when the heat flux was sufficient to give ignition. 

The turbulence extended with a 2 mm radius into the gas 
stream. Although low shadowgraph sensitivity did not per- 
mit observation of the flow pattern immediately adjacent to 
the wire surface, Eckert and Soehnghen (4) reported turbu- 
lence “right up to the surface” during free convection with a 
vertical heated plate. Since the present heat fluxes were 
considerably larger than those used by Eckert and Soehnghen, 
the gas immediately at the heated wire in the experiment 
discussed in this paper was probably also in turbulent motion. 
Therefore, a theoretical treatment of ignition assuming 


laminar flow at the surface would not apply to the present 


~ tially independent of equivalence ratio‘ (Fig. 3). 


work and may in general be of limited usefulness in inter- 
preting ignition during flow over a hot surface. 

Empirically, the minimum total rate of energy input for 
ignition dE/dt increased with increasing initial flow velocity 
and increasing wire or tube diameter (Fig. 2) but was essen- 
The present 
limited results corresponded to 


dE/dt =aV +b 
where V is the initial linear flow velocity, a is a constant which 


_ depends upon the wire and tube diameters but is independent 


of equivalence ratio, and b is a constant which is independent 
of wire and tube diameters and equivalence ratio. The inde- 
pendence of d#/dt on gas composition is probably due to the 
heat capacity of the combustible gas being essentially inde- 
pendent of composition. Interpretation of the effects of wire 
and tube diameter is uncertain at present. The increase in 
dE/dt with increasing flow velocity is undoubtedly related 
to decreasing contact time, but this relationship is obscured 
by the presence of turbulence. Ignition presumably involves 
heating a volume element of the combustible gas to the tem- 
perature range where autoignition readily occurs, but the 
temperature history of the volume element as it passes in 
turbulent flow along the wire is difficult to describe. Further 
work is planned. 

The author is pleased to acknowledge the advice of Dr. L. 


_ E. Line Jr. and I. R. King and the assistance of J. T. Scheu- 


rich, who performed the experimental work. 


3 The stabilized flame could be moved up or down the wire by 
decreasing or increasing the voltage on the wire (i.e., by changing 
the wire temperature). 

4 Stoichiometric air/fuel ratio divided by actual air/fuel ratio. 
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A New Instrument for Measuring 
Atmospheric Density and Temperature 
at Satellite Altitudes 


A. J. DESSLER,' W. B. HANSON,? M. HERTZBERG,? 
D. D. McKIBBIN‘ and R. C. WRIGLEY® 


Missile Systems Division, Lockheed Aircraft Corp., 
Palo Alto, Calif. 


A simple rugged instrument is described which can 
measure, at satellite altitudes: Atmospheric density, at- 
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mospheric temperature and satellite or ballistic missile 
angle of pitch and yaw. When placed in a satellite, the 
present instrument is sensitive enough to operate up to 
altitudes of at least 300 miles. It is not affected by vehicle 
outgassing and will measure only the properties of the 
ambient atmosphere. 


Introduction 


NEW instrument which is capable of the direct measure- 
ment of atmospheric density and temperature at very 
high altitudes has been developed and tested in our Space 
Physics Laboratory. This instrument can measure the at- 
mospheric mass density continuously and hence can provide 
not only the time- and space-averaged density, but also its 
variation with altitude, latitude and local time. Thus, one 
successful satellite flight should yield a detailed description of 
the atmospheric density profile and its variations. 

Satellite tracking data must be averaged over many orbits 
in order to be useful in the derivation of atmospheric density 
data. Therefore, analyses of satellite trajectories have thus 
far yielded values of atmospheric density which are a long 
term space and time average of the true density. The sever- 
ity of this averaging process is indicated by rocket data now 
available.6 Rocket measurements at 125 miles altitude show 
variations in density of a factor of more than 10, depending on 


the time and place of measurement. “4 


The Instrument 


When a high velocity vehicle moves through the upper 
atmosphere, the ambient gas molecules produce an impact 
pressure on its forward side equal to kpv?, where & is a con- 
stant determined by the accomodation coefficient ~ 1, p is 
the atmospheric density, and v is the velocity of the vehicle. 
This impact pressure is measured by a very sensitive micro- 
phone wh.ch is the sensing element for the instrument. The 
impact pressure is a direct measure of the density drag, or, 
since the satellite velocity is independently known, the pres- 
sure is also a direct measure of the density. 

By mechanically chopping the beam of gas molecules 
striking the microphone, an a.c. signal is produced whose am- 
plitude is proportional to pv?. To achieve a narrow band- 
width and hence low noise, the microphone signal is detected 
at a phase sensitive detector with a signai from the chopper 
used as a reference. All of the electronic circuitry necessary 
for this instrument has been completely transistorized. 

If a small hole is cut in the front part of a satellite moving 
through the upper atmosphere, the beam of ambient gas par- 
ticles which passes through the hole will spread, due to the 
transverse thermal velocity components of the individual 
molecules. The amount of thermal spreading will be propor- 
tional to the square root of the absolute temperature. Thus, 
when the microphone is placed a short distance behind the 
sensing hole and moved transversely to map the intensity dis- 
tribution across the beam, a measure of the amount of thermal 
spreading may be obtained. This distribution of the impact 
pressure behind the hole is an absolute measure of the kinetic 
temperature of the ambient atmosphere if the molecular 
weight of the gas is known. Direct methods for measuring 
the ambient gas temperature at altitudes greater than about 
100 km are limited by radiation effects which are larger than 
the heat transfer to a detector by the ambient gas. 

The direction of a molecular beam which passes through 
an aperture in the vehicle’s skin is obviously determined by 
the vehicle’s direction of motion. A suitable combination of 
apertures and detectors can be used to measure this beam 
direction with respect to the body axes and hence provide 
pitch and yaw information. 


6 For a discussion of these data see ‘Atmospheric Density from 
Present Rocket and Satellite Measurements,” by R. A. Minzner, 
presented at the Symposium on Satellite Geophysical Studies, at 
the Washington meeting of the American Geophysical Union on 


May 6, 1958. 
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Calibration 
The sensitivity of the density gage has been measured 
in the laboratory with the aid of a simple molecular beam 


apparatus. The entire gage was placed in a large vacuum 
chamber which was maintained at a pressure of less than 10~4 


mm Hg. An aperture of area 0.1 cm?, which was located 19 
em from the microphone, separated the large vacuum chamber 
from the beam source chamber. The gas pressure in the 
source chamber was controlled by a needle valve and measured 
with a Phillips ionization gage. Argon gas or air at room 
temperature was used for the beam. The product pv? was 
adjusted to simulate the impact pressure which would be en- 
countered in a satellite. Thus, although p and v independ- 
ently do not duplicate satellite conditions, the product pv? 
will be correct. Even a change in the accommodation coeffi- 
cient from 0 to 1 cannot change the product kpv? by more than 
30 per cent. 

In a typical calibration run, a source pressure of 5 X 1078 
mm Hg was used. Allowing for spreading of the beam after 
passing through the 0.1 cm? hole, this source pressure cor- 
responds to a pv? at the microphone of about 6 x 1074 
dynes/em?. The voltage signal to noise ratio at the detec- 
tor was about 4 under these conditions, using a 3 eps bandpass. 
This is the value of pv? which would be experienced in a satel- 
lite moving through a gas of density 10-5 gm/cm*. Taking 
the lowest atmospheric density values of Jastrow,’ 10—-% 
gm/cm? corresponds to an altitude of 375 miles. The voltage 
signal to noise ratio at 300 miles should be better than 12. 
Hence, quite reliable density values should be obtained at this 


altitude. 
Conclusion 


The instrument described above provides a simple rugged 
detector which, with suitable modifications, can be used to 
measure as a continuous function of position: Atmospheric 
density, atmospheric kinetic temperature and vehicle angle 
of pitch and yaw. Since the operation of the detector is 
dependent on the impact of high velocity gas molecules on a 
sensitive microphone element, it will not be greatly affected by 
vehicle outgassing. Laboratory tests indicate that this detec- 
tor, when moving through the upper atmosphere at satellite 
velocities, should operate at altitudes up to at least 300 miles 

in its present form. Calculations show that improvements 
are possible which can extend this altitude capability con- 
siderably. 


Acknowledgment 


We wish to thank Dr. F. 8. Johnson for his efforts in getting 
this project started. We also wish to thank J. Drake, C. 
Searing and W. Page for their aid in the construction of the 
instrument. 


7 Jastrow, R., Bull. Am. Phys. Soc., vol. 3, no. 189, 1958. 


Some Comments on Generalized 
_ Trajectories for Free Falling 


bs of High Drag 


WILLIAM SQUIRE! wah 
Bell Aircraft Corp., Buffalo, N. Y. 
N A recent note (1)? Turnacliff and Hartnett have discussed 
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ponential density variation. This problem was apparently 
Received Aug. 20, 1958. 


! Aerodynamicist, Aerophy sics Department. 
? Numbers in parentheses indicate References at end of paper. 


838 


first treated by Munk (2) in 1944. In 1951 a number of notes 
(3 to 5) appeared on the general problem of a body wth a drag 
proportional to the velocity and velocity squared in an ex- 
ponential atmosphere. 

The treatment given below is essentially based on Munk’s 
analysis. The novel feature is a graphical construction which 
with the use of a variable scale* gives particular solutions 
very rapidly. 


Analysis 
If distance is measured from sea lev el the basic differential 
equation is 


where as in (1) — wi 
=2/a 
and 
mo 


A new variable z = Ky,e~* is introduced, transforming [la] 


into 
di? 


The solution is 


where C is an integration constant and 


is the exponential integral, several tables of which are avail- 


able (6, 7). Besides the expansion 
n= © 5” 
Ei(z) = In ye + [5] 
In y = 0.5772 
. . 
cited in (1), there is an asymptotic expansion i 8 


which is useful for large values of 2. 
Graphical Solutions 


Since o? is a linear function of C, if 5? is plotted for two values 
of C, the value for any other value of C can be found by linear 


_ interpolation or extrapolation which is conveniently per- 


formed with a variable scale. Furthermore since 


a plot against In Z is easily adapted to any value of Kp. 


The values for C = 100 and C = 20 are given in Table 1 
and plotted in Fig. 1. 

The values 20 and 100 taken for C are of course arbitrary, 
and it is a simple matter to make up curves for any ranges of 
interest to the user with the values of C selected for easy in- 
terpolation. 

The first example given in (1) corresponds to 6? = 31 at 


’Gerber Scientific Instrument Co., Hartford, Conn. 
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Table 1* Numerical values for trajectory analysis 
5? 
In Z (C = 100) (C = 20) Ei (2) 
—10.00 81.15 1.15 —9.4228 
—9.00 83.15 3.15 —9.4227 
—8.00 85.12 —7.4225 
—7.00 87.08 7 15 —6.4217 
—6.00 88.85 9.13 —5.4203 
—5.00 90.55 17.09 —4.4165 
—4.605 91.05 11.84 —4.0179 
—3.507 91.42 13.7 —2.8991 
—2.996 90.62 14.52 — 1.6228 
2.303 87.55 15.16 —0.8218 
—1.609 80.53 15.03 +0. 1047 
—0.916 67.16 13.54 +0.7699 
—0.511 55.73 11.82 +1.3474 
—0.223 46.14 10.20 +1.8951 
0.000 38.18 8.75 4.9542 
0. 262 28.74 6.93 2.7214 
0.531 19.70 5.08 2.7214 
0.693 14.87 4.05 4.9542 
1.099 5.97 1.99 3.9210 
1.386 2.55 1.09 19.631 
1.792 0.67 0.48 85.990 
2.079 0.33 0.30 44.038 


* The first six values were obtained by using £;(Z) = 
In yz g. For the others the values were taken from 
Jahnke-Emde (7). 
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£=8. The value of Ky is 1900 deg, so In Kp is 7.55; therefore, 
the value of In z is —0.45, and the value of C is found to be 
58. The rest of the curve is shown dotted in Fig. 1. 


Discussion 


Because of the very large range of densities and velocities 
under consideration at present, it is extremely doubtful that 
the assumption of a constant drag coefficient holds. With 
the present rapid method of computation, it is possible to 
assume a drag coefficient which is constant in a limited range 
of Mach and Reynolds numbers, and piece together a number 
of such solutions. 

In connection with the application made by (1) to oblique 
flight, it should be noted that the path will only be a straight 
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line, as they assume, when the drag is much greater than the 
force of gravity, i.e., when the exponential integral term can be 
neglected. 
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On the Classification of the Chemistry 
in Combustion Experiments 
GERALD ROSEN’ 

Hughes Aircraft Co., Culver City, Calif. 


Nomenclature 


A = constant “frequency factor’ of the Spalding approxima- 
tion, gr-cal/em‘-sec?-K 


Cy = constant specific heat of the mixture, cal/gr-K 

f(T) = volumetric production of the product species, gr/em?- 
sec 

H = C,(T, — T.), heat of reaction, cal/gr 

k(T) = thermal conductivity of mixture, ca! /em-sec-K 

My, = molecular weight of unburned gas, gr/mole 

= pressure, keal/em? 

R = universal gas constant, 1.986 keal/mole-K 3 

T = temperature, K 

T. = temperature of unburned gas, K - 

T, = temperature of burned gas, K 

T =(T — T,)/(T, — dimensionless temperature rise 

Introduction 


D. B. Spalding (1)? has recently considered reaction rate 
expressions for the laminar flame which take the form 


This approximate representation of the overall chemical 
kinetics provides a decisive simplification in the theoretical 
treatment of many combustion problems. The exponent m, 
the well-known effective order of the overall chemical reaction, 
can be determined empirically from the pressure sensitivity 
of the laminar burning velocity (2); it should be remarked 
that m, then, is not necessarily an integer. The exponent n, 
which we shall call the effective index of the overall chemical 
reaction, seems to be an equally appropriate parameter for a 
generic classification of combustion experiments. 

The principal reason for assigning {m,n} labels to a flame 
experiment is that these parameters completely determine the 
approximate dimensionless structure of the flame, as a conse- 
quence of [1] and the flame equations. A knowledge of the 
changes in the {m, n} labels, in response to changes in the 
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: eae of the flame, will provide theorists with a clue for 
understanding the associated changes in flame structure. 

_ The purpose of this note is twofold. First, to show how the 
index can be related to the more venerable parameters of 
chemistry. Second, to give up-to-date theoretical 


Fig. 1 


expressions for the empirical determination of {m, n} labels 
from laminar flame experiments. 


Theory 
According to flame theory (1, 3, 4), the laminar burning 


velocity 


K(T)f(T)dT ] 


_ with the dimensionless eigenvalue parameter \ somewhat less 
than 3 for practical reaction rate expressions. The traditional 
semitheoretical approach estimates the temperature de- 

_ pendent functions with the forms 


K(T) = 


E 
Ba - (7) exp [4] 


Comparing these expressions with the Spalding approximation 
ol ], we see that 


T= exp R [6] 


t will be noticed that [6] has been adjusted so that both 
sides become exactly equal at T = T, (7 = 1). We determine 
the index n by requiring that the derivatives of both sides be- 
come equal at J = T,. This gives 


H H E 


In the figures we have plotted both sides of [6] for the 


| 
The Spalding approximation for typical values of the index n 


typical case of g + h = m, css = maa and for several 
values of the index, n = 5/6 (E/RT;,). Through the entire 
range of physical interest, the semitheoretical and Spalding 
approximations correspond uniformly to within about 2 
percent. (See figures.) 

Now we substitute [1] into [2]. 
obtain 


[ n) Bk 
Ss, = n) exp (g th | 


Ame ) (- *) 
RT, CT, exp 


with the Euler product of Gamma functions 


T(m + 1)T(n + 1) 
T(m + n + 2) 
Thus, if the burning velocity of a combustible mixture is 


measured at different pressures with the same flame tempera- 
ture, we have the influence coefficient 


S,,) 


Using [5 and 7], we 


B(m, n) = 


On the other hand, suppose the burning velocity of a com- 
bustible mixture is measured at neighboring flame tempera- 
tures with the same pressure, by varying the initial tempera- 
ture and density. Since the influence coefficient receives pre- 
dominant contributions from the terms within the braces in 
[8], we have 


{m, labels appear as 
(In Sy) 


LUI 


per 
accur 
tudin 
study 
the a 
ultim 
nome 


Th 
of det 
ing n 
lengtl 
drical 
mixtu 
levels. 
and 
servo 
250 Ik 
amplit 
appro: 
veloci' 
Detail 
descri 
their 


Sp 
nautie: 
5 Nu 


DEcE 


1.0 10 
Co 
6Nn(t-1) 
=) J (In 
Z 
0.2 04 0.6 08 Lo It 
limi 
tudi 
ont 
men 
wide 
assor 
used 
tran: 
publ 
& T \ 
P\n E 
k, exp( — — i 
Ree 
Memb 
3 Ste 


References 


1 Spalding, D. B., “Theory of Laminar Flame Propagation,” 
Combustion and Flame, vol. 1, 1957, p. 287. 

2 Gray, P., Lee, J. C. et al., “Propagation and Stability of 
the Decomposition Flame of Hydrazine,” Sixth Symposium 
(Int.) on Combustion, Combustion Institute, 1957, p. 255. 

3 Rosen, G., ‘An Action Principle for the Laminar Flame,”’ 
Seventh Symposium (Int.) on Combustion, Combustion Insti- 
tute, 1958 (in publication). 

4 Rosen, G., “Generalization of the Laminar Flame Action 
Principle for Arrhenius-type Rate Expressions,’’ Combustion and 
Flame, vol. 3, Jan. 1959 (in publication). 


On the Importance of the Sensitive Time 
Lag in Longitudinal High-Frequency 
Rocket Combustion Instability 


LUIGI CROCCO,! JERRY GREY? and DAVID T. HARRJE® 
Princeton University, Princeton, N. J. 


It was determined experimentally that there is an upper 
limit to the chamber length at which each mode of longi- 
tudinal high-frequency pressure oscillations will occur, 
and that this limit was accurately predicted for the funda- 
mental mode by Crocco’s sensitive time lag theory over a 
wide range of mixture ratios. The technique of using 
longitudinal stability boundaries for simple and direct ex- 
perimental determination of the sensitive time lag and 
associated parameters was established as valid and may be 
used to replace the complicated and difficult chamber 
transfer function measurements described in previous 
publications. 


Introduction 


N EXTENSIVE experimental study was performed on a 
small uncooled rocket motor using liquid propellants (95 
per cent ethy] alcohol and liquid oxygen) in order to establish 
accurately the nature of the stability boundaries for longi- 
tudinal high-frequency combustion pressure oscillations. This 
study comprised a small part of an overall program‘ aimed at 
the attainment of sufficient systematic knowledge to provide 
ultimate elimination of rocket combustion instability phe- 
nomena. 


Experimental Data 


The experiments discussed in this note consisted basically 
of determining the effects on high-frequency stability of vary- 
ing mixture ratio and combustion chamber length. The 
length changes were achieved by adding and removing cylin- 
drical sections of the experimental chamber (1),° and the 
mixture ratio was changed by adjustment of injector pressure 
levels. Steady-state chamber pressure (nominally 300 psia) 
and mixture ratio were maintained constant on each run by a 
servo control system (2). Thrust level was approximately 
250 lb. The experimental data included measurements of 
amplitude and frequency of combustion pressure oscillations, 
approximate axial distribution of steady-state chamber gas 
velocity, and the usual steady-state performance parameters. 
Details of the experimental apparatus and methods have been 
described in a number of previous publications (e.g., (1,2) and 
their bibliographies) and need not be repeated here. 

A brief summary of the experimental results is presented in 
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Fig. 1, which plots the observed stability boundaries as func- 


tions of mixture ratio and chamber length. The experi- 
mental definition of these boundaries is demonstrated in the 
three-dimensional plot of Fig. 2 (which includes, for clarity, 
only the fundamental mode), showing the measured ampli- 
tude of pressure oscillations as a function of length and mix- 
ture ratio. Fig. 2 and photos of the observed wave shape 
demonstrated clearly that shock-type pressure waves occurred 
quite close to the stability boundaries. 
Discussion 

The experimental results of Figs. 1 and 2 corroborate 
directly the Crocco theory’s well documented qualitative 
prediction (e.g. (1, 2) and Ref. 1, 2 of (4) etc.) that there 
exists an upper limit of chamber length beyond which any 
given longitudinal mode of oscillations cannot develop, as 
well as the lower limit which has been observed by a number 
of other experimenters (e.g., Berman, Tischler, Zucrow and 
their respective collaborators). Note that this existence of an 
upper limit clearly invalidates the oscillation-producing 
mechanism advanced by Zucrow and Osborn (3) as discussed 
by Crocco (4) since this mechanism cannot account for the 
observed cessation of oscillations with increasing chamber 
length (see Figs. 1 and 2). 

Besides being the only mechanism which predicts qualita- 
tively the upper longitudinal stability boundary, however, 
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Fig. 2 Three-dimensional plot of the experimental fundamental 

mode longitudinal stability boundary, showing the definition ob- 

tained. Nominal chamber pressure is 300 psia. (Second and 
third harmonic boundaries are omitted for clarity) 
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Fig. 3 Values of sensitive time lag 7 and interaction index n ob- 
tained by using experimental data from the lower first mode 
stability boundary 


the theory based on the sensitive time lag concept as first 
outlined by Crocco is unique in its capability for quantitative 
prediction of stability behavior. Previous publications giving 
results of the Princeton program have demonstrated that the 
theory has predicted quite well the effects on longitudinal 
stability boundaries of nozzle configuration and chamber 
pressure (e.g. (1) and Ref. 5 of (4), etc.) using measurements 
of the time lag and its associated parameters obtained for a 
few cases from chamber transfer function determinations (2). 
A more general, simpler and more direct experimental proof 
of the quantitative accuracy of the sensitive time lag theory 
as developed in detail in Ref. 2 of (4) will now be given, 
using the experimental data which appear here in Fig. 1. 

It is demonstrated in the cited reference that for the gen- 
eral case of longitudinal pressure oscillations in a rocket 
motor,® the theoretical conditions for neutral oscillations (i.e., 
the stability boundaries) may be expressed by 


n = N(w, a, u(x), L, k) 


T a, y, u(x), L,é,k,m) 


¥ 


sensitive time lag 

interaction index relating time lag to the chamber 
parameters (primarily pressure and temperature). 
It represents in a synthetic way the degree of sensi- 
tivity of 7 to variations in the chamber conditions 

observed frequency at which neutral oscillations occur 
(i.e., at the stability boundary) 

complex nozzle admittance parameter, calculated by 
the theory of (5) and conclusively verified by ex- 
periments described in (6) and later papers 


° Fortunately, it no longer appears absolutely necessary to 
justify the application of a linearized analysis to so nonlinear a 
process as the shock-type pressure oscillations observed. The 
use of linearized equations to predict stability boundaries, at 
which no finite oscillations exist, is a well-known technique ap- 
plied in the instability analyses of practically every field of dy- 
namics. Further, it has been demonstrated experimentally for 
this rocket motor configuration that longitudinal shock waves, 
artificially produced by cartridges, had a negligible effect on the 
stability boundaries measured in the undisturbed rocket chamber 
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Fig. 4 Comparison of theoretical prediction with experimental 

observation of the upper fundamental mode stability boundary 

The theoretical prediction is based on values of r and n from Fig. 
3, obtained from measurements at the lower boundary 


Y = mear value of specific heat ratio of combustion gases 

u(x) = axial chamber gas velocity distribution 

x = axial space coordinate 

L = chamber cylindrical length 

c = mean speed of sound in chamber 

m = mode of longitudinal oscillation (e.g., m = 1 corre- 
sponds to fundamental mode, m = 2 to second 
harmonic, etc.) 

k = estimated liquid droplet drag coefficient 


Details of Equations [1] and the definitions of some of the 
above parameters are not essential to the present discussion, 
and in the interest of brevity, since they are thoroughly docu- 
mented in the references cited, no further comment will be 
made here. It suffices to note that all quantities on the right- 
hand sides of Equations [1] were readily measured or inferred 
experimentally at the stability boundaries shown in Fig. 1. 

First, Equations [1] were used to calculate the critical 
values of 7 and n, using experimental data measured at the 
chamber lengths corresponding to the lower stability boundary 
(fundamental mode) of Fig. 1 at a number of the observed 
mixture ratio points. (Calculations to cover the full mixture- 
ratio range of Fig. 1 are now in process.) These values’ of n 
and the sensitive time lag 7 are shown in Fig. 3. 

Note that the values which appear in Fig. 3 do not in 
themselves constitute any check on the theory, since they are 
merely calculated from the measured lower stability boundaries 
of Fig. 1 and the theoretical Equations [1]. However, if 
the theory used to obtain the values of 7 and n shown in Fig. 
3 is correct, and if 7 and n are, as the theory postulates, a 
valid mechanism for description of the rocket combustion 
processes, then at each experimental mixture ratio, Equations 
{1] solved for L and w using the 7 and n values of Fig. 3 
should be capable of predicting quantitatively the chamber 
length and frequency at which the upper stability limit should 
appear. 

Fig. 4 compares the theory’s predictions of the upper first 
mode stability boundary with the experimental findings, re- 
plotted from Fig. 1. Quantitative calculations for the higher 
modes, whose existence is predicted by the theory, are now in 
process and will be presented in the complete report. The de- 
gree of correlation shown in Fig. 4 provides satisfactory 
verification of the sensitive time lag theory. 


7 It should be pointed out that good agreement had been found 
between these calculations and previous measurements of 7 and 
at one mixture ratio by the more complicated method involving 
chamber transfer function measurement, as indicated in earlier 

ublications, e.g., (1, 2). This, in effect, had constituted a 
imited but nevertheless quantitative verification of the theory. 
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Conclusions 


1 It was observed conclusively by experiments on oxygen- 
alcohol hardware that there exists an upper limit to the 
chamber length above which a mode of longitudinal high- 
frequency pressure oscillation cannot occur in a rocket motor, 
as predicted by Crocco’s sensitive time lag theory. This 
upper limit cannot be explained by any other mechanism 
advanced to date. 

2 The sensitive time lag theory as developed in Ref. 2 of 
(4) has been used to predict quantitatively the upper stability 
boundary (i.e., chamber length) of the fundamental mode of 
longitudinal high-frequency oscillation over a wide range of 
mixture ratios. These results constitute a considerable ex- 
tension of the previously established experimental verifica- 
tion for the lower fundamental mode stability boundary at one 
mixture ratio. 

3 The longitudinal stabilitv-boundary method of deter- 
mining the combustion parameters 7 and n, used to establish 
the above results, provides a remarkably simple and accurate 
determination of these important high-frequency stability 
criteria. This technique is expected to replace the former 
difficult and complicated transfer-function method described 
in previous publications. ao 


Extension of Results 


The above results were obtained for a single chamber 
pressure, propellant combination and injector type. Simi- 
lar experiments for different chamber pressures, widely 
different injector characteristics, and several combinations 


of fuels and oxidizers have been in process for some time. 
This rather complete survey of longitudinal rocket instability, 
of which the results given in this paper form a small part, will 
be included in a more comprehensive report at an early time. 
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Technical Comment 


Comments on the Zucrow-Osborn Paper 
on Combustion Oscillations 


LUIGI CROCCO! 


HE idea that the phenomenon of high-frequency com- 
bustion instability in liquid propellant rockets should be 
related to the effect of pressure waves on the local rates of 
combustion was first introduced in a paper which I published 
in this JouRNAL in 1952 (1).2- By means of the analytical 
treatment of a simple combustion model I was able to show, in- 
deed, that the appearance of instability is determined by (a) a 
space condition, that is, the regions of maximum energy re- 
lease must be sufficiently far away from pressure nodes, and 
(b) a time condition, that is, the ratio between a characteristic 
combustion time, called the sensitive time lag, and the period 
of the oscillations must lie in certain ranges. Observe that 
the sensitive time lag is not the same as the total time lag (on 
which the theory of low frequency rocket combustion instability 
is centered), and actually it must be much shorter in order to 
justify the high frequencies encountered. These ideas, ap- 
plied later to more complicated models, are summarized in a 
monograph (2) which also includes the experimental con- 
firmations available at the time. Many additional publications 
by various authors, both in the United States and abroad, 
have shown that in general the mechanism on which these 
theories are based has been accepted as correct. 
A different mechanism, which does not involve the idea of 
a time lag, has been advanced in a recent publication (3) by 
Zucrow and Osborn, based on the interpretation of experimen- 


Received Nov. 3, 1958. 

' Robert H. Goddard Professor of Jet Propulsion. 
Member ARS. 

* Numbers in Parentheses indicate References at end of paper. 


Fellow 


Princeton University, Princeton, N. J. 


tal results obtained through a technique similar to that al- 
ready used in recent years at Princeton (4). The main dif- 
ference in the experiments at Purdue consists of the use of 
premixed gaseous propellants, and promises to furnish in- 
teresting comparative results by eliminating the processes 
of atomization, evaporation and mixing always present when 
the propellants are liquid and unmixed. The conclusions 
reached by Zucrow and Osborn can be summarized as fol- 
lows: 

1 The unstable burning observed is due to the effect of 
pressure waves in increasing the local burning rates. 

2 The tendency toward instability is greater when the 
amount of propellant injected during one period of oscillation 
is larger, and therefore increases with decreasing frequency. 

3 As a result of argument 2, and in agreement with ex- 
perimental observations, a critical combustion chamber length 
exists below which the frequencies are too high for an appre- 
ciable pressure oscillation to be produced, and above which 
the always decreasing frequency provides the conditions for 
unstable operation with a pressure amplitude that steadily 
increases with increasing chamber length. 

4 The same argument can be applied to transversal types 
of instability to explain the aptitude of large rocket motors of 
small aspect ratio to exhibit this type of instability. 

5 The shock-fronted pressure waves traveling in the com- 
bustion chambers are a type of detonation wave. 

Analyzing these statements in order, we see that: 

1 This suggestion is exactly the same as the one advanced 
in my 1952 paper (1), and since then reformulated by 
various authors. 

2 This argument does not stand a more critical examina- 
tion, since in the appearance and maintenance of the unstable 
process what counts is not the absolute value of the excess 
energy released during one cycle as a consequence of the 
modified combustion rates, but the balance between this ex- 
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cess energy and the excess energy absorbed during the whole 
cycle by dissipative processes and by the nozzle. With de- 
creasing frequency both the released and the absorbed excess 
energies increase, and hence the balance is not particularly 
sensitive to the frequency level. In other words, it is not the 
excess energy during one cycle which determines the possibil- 
ity of unstable combustion, but rather the excess energy per 
unit time, or the excess power. 

3 This statement is contradicted by the results of some 
careful experiments carried out at Princeton. The first in- 
dications from these experiments have already been released, 
e.g., (4, 5 and 6). More complete results will soon be pub- 
lished in a full report now in preparation, but a preliminary 
release was believed necessary at this stage and is the object of 
a separate Technical Note in this same issue (7). It is clearly 
shown by these results that it is not true that only a lower 
critical length exists, above which oscillations become more 
and more severe, but that for each mode of longitudinal os- 
cillation there exists a range of lengths outside which the 
operation is stable, while inside it is unstable. Thus, by suf- 
ficiently increasing the chamber length, the operation be- 
comes stable again for a given mode (though it may still be 
unstable for higher modes). This result invalidates the 
mechanism suggested by Zucrow and Osborn (if a proof is 
still needed of its fundamental incorrectness). It also is a 
clear confirmation of the existence of a characteristic combus- 
tion time and of its fundamental importance in determining 
whether a rocket will operate under stable or unstable con- 
ditions. The analysis of the data actually allows the deter- 
mination of this characteristic combustion time (which we 
have called the sensitive time lag, or simply the time lag), 

and of other important characteristic quantities related to 

combustion, and gives a proof of the fundamental correctness 
of the combustion model which is the basis for the existing 

theoretical treatments of rocket combustion instability (2). 

This is shown in the Technical Note in this issue (7). 

4 Systematic quantitative experimental information about 
transversal instability is not yet available., The Princeton 
group and others in this country are actively engaged in a 
effort to provide data of this kind. But even before these re- 
sults become available, I believe I may predict that no matter 
what is the actual mechanism responsible for instability, the 
time lag will again play a fundamental role, and that its ratio 
to the period of the oscillations will have to lie in well deter- 
mined ranges if instability is to occur. The fact that in large 
rocket motors unstable operation is more likely to appear is 
probably related to the fact that in a chamber of no matter 
what shape there exists an infinite series of frequencies of 
natural oscillation corresponding to various modes and start- 
ing from a well determined minimum value. For small 
chambers, this minimum frequency may be such that the 
corresponding period of oscillation is too short compared to 
the sensitive time lag, and for no oscillation mode will the 
chamber be in a condition conductive to unstable operation. 
If, however, one of the dimensions of the chamber is increased, 
lower minimum frequencies and longer oscillation periods 
may be produced, thus upsetting the balance with the time 

_ Jag (if, as is likely to happen, this has not changed in the same 
_ proportion) and favoring the appearance of instability. If 
_ the length of the chamber is increased, the instability, when 
produced, is of the longitudinal type. If the diameter is in- 
-ereased, transversal instability is the result. By further in- 
creasing the dimensions, other modes may always produce 
_ frequencies in the range proper for the appearance of instabil- 
ity, so that it is quite possible that stability cannot be found 
- again by just increasing the dimensions, but that, in order to 


obtain stable combustion with large chambers, the overall 
energy balance must be modified through a change in the 
space distribution of combustion or an increase of damping. 
This may be particularly necessary in the case of transversa!| 
instability, since an unpublished theoretical study shows that 
the damping effect of the nozzle is absent for transversal os- 
cillations. In principle, there would also be the possibility of 
obtaining stability through a suitable increase of the time 
lag. But in practice, this may be either impossible (if th. 
sensitive time lag is essentially determined by the chemica| 
properties of the propellants, as it certainly is for premixed! 
gaseous propellants) or undesirable (because it may result ii 
the necessity of a larger chamber to maintain the performance: 
level). 

5 The concept that the waves propagating in the chamber 
are a kind of detonation wave is somewhat misleading. I: 
fact, in a rocket chamber the shock wave (when it exists 
merely accelerates the rates of a combustion process whic): 
would take place even in the absence of the shock wave, whil« 
in a detonation wave the chemical change is initiated by th: 
passage of the shock wave, on which it essentially depends 
Moreover, the velocities and amplitudes of the shock waves 
are much below the values they would achieve for an ordinar) 
detonation wave (which, by the way, are very well defined fo: 
the case of premixed gaseous propellant as used in the Purdu: 
experiments). Finally, to introduce the concept of this type 
of detonation wave as essential to the phenomenon of insta- 
bility would make a different mechanism necessary for th« 
case when the waves are of a sinusoidal nature, as they are in 
many instances of fully developed instability, and always at 
the onset of a self-started instability (see, for instance, Fig. 5 
of (3)). 

It is surprising that a highly competent group such as that 
headed by Professor Zucrow at Purdue University should ad- 
vance a new mechanism to explain the combustion instability 
in rockets, without any reference to previous formulations ex- 
tensively used in the published literature. Perusal, even 
limited, of this literature would probably have indicated to 
the authors the incorrectness of their conclusions. a 
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Why wait for the answers? 


... Observe and evaluate phenomena as they 
occur—with new Kodak Linagraph Direct Print Paper! 

This new photorecording paper gives immediate readout in 
moving-mirror galvanometer oscillographs—with no photochemical 
developing. Gives sharp, legible traces over a wide range of 
recording frequencies—from 0 cps to 3000 cps. 

Linagraph Direct Print records can be used in recording 
system readers, accept pen and pencil notations readily. They 
have adequate permanence and may be stabilized 
for extended archival storage. 

Now available in 5” x 200’, 6” x 100’, 7” x 200’, and 12” x 200’ rolls. 
Other sizes on request. For complete details, write: 


EASTMAN KODAK COMPANY _ 


Graphic Reproduction Division 
Rochester 4, N. Y. 
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Object locating (2,842,760). 
J. L. MeLueas, State College, Pa., 
assignor to Haller, Raymond and Brown, 


ne. 
i Means for selectively energizing a cath- 
ode-ray tube control grid with the output 
of a radar receiver or the output of a 
facsimile pickup head. 
Moving target indicator radar system 
(2,842,761). J. M. Downs, Glen Cove, 
N. Y., assignor to Sperry Rand Corp. 
Detector connected to a receiver with 
means to delay reflected signals at least 
one cycle of pulse repetition frequency. 
Subtracting and adding means integrate 
delayed and undelayed versions of the 
signals. 
Collision warning radar (2,842,764). N. 
L. Harvey, Eggertsville, N. Y., assignor 


Eprror’s Norte: Patents listed above 
were selected from the Official Gazette of 
the U.S. Patent Office. Printed copies of 
patents may be obtained from the Com- 
missioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 


Inc. 
Ranging apparatus with bandpass filter 
tuned to reject all frequencies outside 
the band corresponding to the Doppler 
shift in the echo-signals frequencies re- 
ceived. 


Infrared transmitting mirror (2,852,980). 
H. Schroder, Munich, Germany. 

Mirror for reflecting visible rays from a 
light source and transmitting the heat 
rays. Reflector consists of layers of a 
transparent material nonabsorbent of 
infrared rays, and alternate layers of 
material having a high index of refraction. 


Swingaway support for missiles (2,- 
852,981). C. A. Caya, Los Angeles, Calif. 

U-shaped suspension member on a 
bomb rack engaging a hood on the we 
of the missile and permitting the hoo 
to swing from an upwardly extending 
position to a folded position within the 
missile upon release. 


Aerial carry and release mechanism 
(2,852,982). C. W. Musser, Philadelphia, 
Pa., assignor to the U.S. Army. 

Hooks in holding grooves in a pair of 


supports released by firing a cartridge- 
actuated initiator which moves a piston 
in a cylinder. Mechanical linkage from 
the piston disengages the hooks attached 
to the store to be released. 

Control apparatus (2,853,255). A. P 
Rasmussen and F. J. Huddleston, Millers- 
ville, Md., assignors to Westinghouse 
Electric Corp. 

System for controlling a control surface 
of an aircraft. A decoupling link is 
placed between a manual controller and a 
servo controller. Means responsive to 
the power supply operates the decoupling 
link from one position to the other. 
Aircraft ejection seat (2,853,258). O. F. 
Polleys, Windsor, Conn., assignor to 
Chance Vought Aircraft, Inc. 

Capsule, with a body and cover, spaced 
behind the pilot’s seat. The seat may be 
moved into the capsule body, completely 
covered and ejected from the aircraft. 
Motion measuring system (2,853,287). 
C. S. Draper and C. Emmerich, 
Cambridge, Mass., assignors to Research 
Corp. 

Viscous damped mass element mounted 
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our illustrated brochure to see how our facilities can be put to 
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Wrap-around 
bumper for 
a continent 


Like a huge ‘“‘bumper” wrapped around the North 
American continent and reaching down along both 
the Atlantic and Pacific shores, the North Ameri- 
can Air Defense Command (NORAD) has been 
created for operational control of air defense 
units of the Army, Navy and Air Force of the U.S. 
and the RCAF Air Defense Command of Canada. 
Its field includes the vast area between the 
southern border of the United States and the 


Headquarters—NORAD—Colorado Springs 


northernmost limits of Canada and Alaska. Under 
the functional control of NORAD will be BMEWS 
(Ballistic Missile Early Warning System) and 
SAGE (Semi-Automatic Ground Control Environ- 
ment) for the defense of specified sectors. In addi- 
tion to its responsibility as prime contractor for 
BMEWS, the Radio Corporation of America is 
working on other important electronic assign- 
ments for NORAD. 
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ENGINEERS 
& SCIENTISTS 


Here is your opportunity to grow 
with a young, expanding subsidiary 
of the Ford Motor Company. Out- 
standing career opportunities are 
open in Aeronutronic’s new RE- 
SEARCH CENTER, overlooking the 
Pacific at Newport Beach, and the 
facility in dale, California. 
You will have all the advantages of a 
stimulating mental environment, 
working with advanced om 
in a new facility, located where you 
=~ enjoy California living at its 
inest. 


PhD and MS RESEARCH SPECIAL- 
ISTS with 5 to 7 years’ experience in heat 
transter, fluid mechanics, thermodynamics, 
combustion and chemical kinetics, and 
thermoelasticity. To work on theoretical 
and experimental programs related to re- 
entry technology and advanced rocket 
propulsion. Specific assignments are open 
in re-entry eg design, high temperature 
materi undary layer heat trans- 
fer with chemical reaction, thermal stress 
analysis, and high temperature thermo- 


ynamics. 
PROPULSION ENGINEERS with 5 years’ 
experience in liquid and solid rocket design 
and test. Familiarity with heat transfer 
problems in engines desirable. To work on 
program of wide scope in R & D of advanced 
concepts in rocket one components, and 
for missile project work 
ADVANCED AERODYNAMIC FACIL- 
ITY DESIGNER. Advanced degree de- 
sired. To supervise work in design and in 
instrumentation of advanced aerodynamic 
test facilities such as shock tubes, shock 
tunnels, plasma-jets, and hyper-velocity 


guns. 
STRUCTURAL ANALYSIS SECTION 
SUPERVISOR with 8 to 10 years’ experi- 
ence, including supervision, in the missile 
field. Graduate degree for design and 
required. Will be required to 
and methods, thermal stress, dynam- 
is, etc. to hypersonic vehicles, 
e-entry bodies, and ce vehicles. 
FLIGHT TEST & STRUMENTATION 
INGINEERS with 5 to 10 years’ experi- 
cone in laboratory and flight test instru- 
ill develop tech- 
niques utilizing advanced instrumentation 
associated with ce vehicles. 
THEORETIC. 
Advanced degree and at least 5 years’ ex- 
perience in high-speed aerodynamics. 
Knowledge of viscid and inviscid gas flows 
required. To work on program leading to 
advanced missile configurations. Work in- 
volves analysis of the re-entry of hypersonic 
missiles and space craft for determining 
optimum 
DYNAMICIST. degree, ap- 
plied mathematics background, and ex- 
perience in missile stability analysis de- 
sirable. Work involves re-entry dynamics 
of advanced vehicles and dynamic analysis 
of space craft. 
ENGINEER or PHYSICIST. With ex- 
perience in the use of scientific instru- 
ments for making nye measurement. 
Work related to flight test and facility 
instrumentation. Advanced degree de- 
sired with minimum of 3 years of related 
experience. 


Qualified applicants are invited to send 
resumes and inquiries to Mr. L. R. Stapel. 


AERONUTRONIC SYSTEMS, 
INC. 


A subsidiary of Ford Motor Company 


1234 Air Way, Bidg. 19, Glendale, Calif. 
CHapman 5-6651 


AERODYNAMICIST. 
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for deflections in a single degree of free- 
dom. A generator signals any deflection 
from a neutral position, and a force 
generator is energized in a direction to 
reduce the deflection of the element. 

Target position indicator (2,853,701). 
J. Freedman, R. O. Schlegelmilch and H. 
Sherman, Lexington, Mass., assignors to 
the U.S. Air Force. 

Transparent screen, with translucent 
backing, on which recorded data are 
projected. Plotter’s marks on the screen, 
together with projected recorded data, 
are sharper and clearer than the pro- 
jected data. A rapid process camera is 
set to be sensitive only to the more 
visible plotter’s marks and _ insensitive 
to the projected image. 


Moving target cancellation circuit (2,- 
853,702). M. C. Johnson and R. F. 
Ahrens, Eatontown, N. J., assignors to 
the U.S. Air Force. 

Modulation of the beam of a cathode- 
ray tube in accordance with an alternating 
voltage, input and output circuits and a 
transformer with two primary windings 
coupled to a secondary winding. Voltages 
induced in the secondary winding by 
the input circuit currents cancel, and the 
output circuit is coupled to the secondary 
winding. 


Video train bracketed by time-spaced 
control pulses (2,853,703). H. T. Hayes, 
Long Beach, Calif., assignor to Gil- 
fillan Bros., Inc. 

Generator for developing range marks 
and modulating them in accordance 
with voltages representing the scanning 
movements of a radiated antenna beam. 
Cursor voltages derived electronically 
establish a predetermined glide path 
ag line in relation to runway course 
ine. 


Radio direction finders (2,853,704). A. J. 
Ortusi and A. Robert, Paris, France, 
assignors to Compagnie Generale de 
Telographie Sans Fil Corp. 

System for determining the angular 
position of a first point with respect to a 
second point. UHF energy is emitted 
over a single carrier wave from one of the 
points. Energy is alternately concen- 
trated at the second point. Current 
from the two points is measured, and the 
ratio between the currents establishes 
the distance relationship between the 
points. 


Jet engine thrust control (2,853,851). 
M. E. Chandler, New Britain, Conn., 
assignor to Pratt & Whitney Co. 

Engine mounted in movable relation- 
ship to the aircraft and moved in propor- 
tion to the reaction thrust, so the air- 
craft is propelled in a definite selected 
speed throughout the speed range of the 
aircraft. 


Boundary layer control for aerodynamic 
ducts (2,853,852). A. G. Bodine Jr., 
Van Nuys, Calif. 

Control of the aerodynamic-acoustic 
vibration originating in a layer of air next 
to the wall surface of a jet engine duct 
through which air travels at sonic speed. 
Frequencies in the region are controlled 
by a sound wave attenuator mounted 
outside the duct, a part extending through 
the wall and communicating with the 
air layer in the vibration area. 


Rocket (2,853,946). A. C. Loedding, 
Princeton, N. J., assignor to Unexcelled 
Chemical Corp. 

Lightweight shell with elongated propel- 
lant charges each having a polygonal 
cross section and operating at a pressure 
of 500 psi. Spaces between charges are 
filled with an elastic inhibitor to form a 
solid unit between the shell and charges. 


to the U.S. Navy. 


Acceleration pressurized bipropellant 
liquid fuel rocket (2,850,975). BD. 
Teague Jr., Alpine, N. J., assignor to 
Bendix Aviation Corp. 

Thrust to the rocket is initiated by a 
booster charge using either hypergolic o1 
nonhypergolic propellants, avoiding the 
need of a staging system in the injector 
system. 

Thrust modifying device (2,847,822). 
G. F. Hausmann, Conn. 
assignor to United Aircraft Corp. 

Streamlined fins upstream of the nozzle 
of a jet powerplant, the rear portions 
movable for diverting the gases trans- 
versely of the duct axis. 

Reversible thrust nozzle construction 
(2,847,823). T. L. Brewer, Ridgewood, 
N. J., assignor to Curtiss-Wright Corp. 

Tail pipe extension for jet engines with 
vanes shaped to close the side walls and 
pivotally movable to direct the exhaust 
flow outward. 

Wing tip jets (2,848,181). W. L. Landers, 
Smithburg, Md., assignor to Fairchild 
Engine and Airplane Co. 

Pair of auxiliary jet engines mounted at 
the wing tips providing supplemental 
thrust to the airplane. Failure of one 
of the wing tip engines automatically 
cuts off the fuel supply to the engine on 
the opposite wing tip. 

Ultrasonic apparatus for the non-destruc- 
tive evaluation of structural bonds 
(2,851,876). J. S. Arnold, Palo Alto, 


Calif., assignor to the U.S. Air Force. 

Ultrasonic frequency modulated cycli- 
cally over a narrow band of frequencies 
outside the resonant frequencies of the 
structure being tested. 


Rocket motor with recrvstallized silicon 
carbide throat insert (2,849,860). E. C. 
Lowe, Niagara Falls, Ontario, Canada, 
assignor to Norton Co. 

Silicon carbide used for throat insert 
because of its high resistance to the non- 
oxidizing flame of a rocket motor reaction 
blast. The material can wthstand tem- 
peratures as high as 2250 deg and is 
extremely erosion resistant. 

Pressure vessels (2, ey 133). E. M. 
Ramberg, Glen Head, N. Y. 

Cylindrical container fluids, con- 
structed of a solid assembly of non- 
metallic fibrous strands of threads and 
helical wrappings. A solidified bonding 
material impregnates the threads, and a 
layer of metal foil is interposed between 
some of the wrappings. 

Temperature responsive fuel control for 
gas turbine (2,848,868,). R. W. Jensen, 
Encino, Calif., assignor to the Garrett 


orp. 

Needle valve for varying the flow of the 
fuel supply through a restricted inlet of 
one of two parallel flow paths, in response 
to temperature changes of the combustion 
air. 
Blast release detent (2,848,925). 
Hood, Sherman Oaks, Calif., 


assignor 


Rocket launcher casing and_ tube 
adapted to contain a rocket having tail 
elements and a retaining element at the 
rear end. Blast from an ignited rocket 
releases a detent from the retaining 
flange engaging position. 
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OF A SERIES 


Unique Oldsmobile-developed two-s age automatic 
choke is a major step forward in improving 
automobile operating economy. 


One of the important carburetor developments during the 
past few years was the automatic choke, a device that 
allows the automobile to be started in cold weather, and 
then keeps it running until the engine is sufficiently 
warmed up to sustain itself. Every automatic choke has 
two separate functions: 1) choking, which enriches the 
fuel-air mixture for starting, and 2) the idle speed con- 
trol, which keeps the engine from stalling once it is 
started. In the past, and on all present carburetors except 
those used on the 1959 Oldsmobile, these two functions 
have operated simultaneously with the result that the en- 
gine ran on a rich mixture for the same length of time 
that the fast idle was “on”. This resulted in excess fuel 
consumption. 


With the introduction of the 1959 Oldsmobile, the two 


functions have been separated with a new and exclusive 


OLDSMOBILE > 


TWO STEPS TO NEW FUEL ECONOMY 


two-stage automatic choke developed by Oldsmobile- 
engineers. An ingenious system of over-running levers 
allows the choke fly to open 75% sooner than previ- 
ously required. The fast idle, however, remains “on” for 
the full warm-up period 
so the engine will not 
stall. This early elimina- 
tion of the choking 
function represents a 
considerable fuel saving 
in cold weather when 
numerous short trips 
are made. 


At Oldsmobile the Inquiring Mind is always at work, 
finding new and better ways to design, engineer and build 
finer automobiles for the most discriminating of buyers— 
the Oldsmobile owner. Discover the difference for yourself 
by visiting your local Oldsmobile Quality Dealer and taking 
a demonstration ride in a 1959 Oldsmobile. 


OLDSMOBILE DIVISION, GENERAL MOTORS CORPORATION 


Pioneer in Progressive Engineering | 
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Thermodynamics of 


One-Component 
. Systems, by William N. Lacey and 
Bruce H. Sage, Academic Press, Inc., 
New York, 1957, xi + 376 pp. $8.00. 
Reviewed by Joun F. LEE 
North Carolina State College 


The author states in the preface that 
this is a textbook intended to “help the 
_ student in his transit from science studies 
_to engineering, from idealized thermo- 
_ dynamics to the combinations of thermo- 

_ dynamics with mechanics needed for 
_ dealing with everyday problems of the 
engineer.’’ However, tue book fails in 
itS avowed purpose as a textbook and 
appears more appropriate as a review 
volume for an engineer or scientist who 
has already had a thorough and complete 
course in the basic thermodynamic 

_ principles and their applications. Part I 
of the book is devoted to a summary in 

rather a catalog form of the more impor- 
basic thermodynamic _ principles. 
“Part II is devoted to applications in the 
restricted area of the traditional flow 

_ processes covered in most thermodynamics 

_ textbooks. The applications are covered 
most cursory fashion. The coverage is 
not only limited in scope, but is even 
more severely limited in depth. The 
basic principles of thermodynamics are 
stated concisely, accurately and clearly. 

The interpretation, the subtleties and 

_ the implications of the basic principles 
are left to the creative imagination or 
ingenuity of the reader. There is consid- 
_ erable merit in this approach, but one 
- eannot avoid the feeling that the authors 

_ have previously carried it to an extreme. 
_ One might expect that in Part II some 
_ Of the deficiencies mentioned previously 
presenting the thermodynamics princi- 
_ ples in Part I might be alleviated. How- 
_ ever, the applications covered in Part II 
are linked only inferentially to the princi- 
_ ples discussed in Part I. It is difficult to 
- see how the material presented in this 
- book could serve as a basis for either 
scientific or engineering action on the 
_ part of a student or practicing engineer 
Neverthless, despite the shortcomings 
mentioned in this review, the volume 
_ should serve as a useful topical summary 
_ of thermodynamics for students or engi- 
neers whose interests fall within the 

_ scope of the book. 


by R. Hanbury Brown and A. C. B. 
Lovell, John Wiley & Sons, Inc., New 
York, 1958, 207 pp. $6.50. 
Reviewed by J. G. Botton 
California Institute of Technology 


_ This book gives an excellent account of 
the first twelve years work in the new 
science of Radio Astronomy. The authors 
“are two men who have directed the efforts 
of the Jodrell Bank Experimental Station 
of Manchester University, which is one of 

- the world’s foremost radio observatories. 


It is written chiefly from an experimen- 
talist’s point of view and can be under- 
stood by anyone who can follow articles in, 
say, the Scientific American. The book 
is copiously illustrated with explanatory 
diagrams and photographs, and only a 
bare minimum of mathematics is used. 

The book is divided into 11 chapters. 
The first two contain a brief résumé of the 
necessary astronomical background and 
the factors governing the reception of 
extraterrestrial radio waves through the 
Earth’s upper and lower atmosphere. 
The third chapter gives an account of 
some of the specialized antenna and 
receiver techniques used in radio astron- 
omy. It is, however, principally confined 
to the meter wave lengths; that is 
the range in which the authors have been 
accustomed to work. It is somewhat 
surprising that the authors make no 
mention of the new high-frequency 
devices, such as the maser and the para- 
metric amplifier, which promise such 
great increases in receiver sensitivity for 
the future. 

There are then seven chapters on the 
results of various phases of radio astro- 
nomical investigation. These chapters are 
on galactic and extragalactic radio emis- 
sions (including the “radio stars’’), the 
hydrogen line, the scintillation of radio 
stars, solar radio waves, meteors, radio 
and the aurora borealis, and finally radio 
investigations of the moon planets and the 
Earth satellites. The last three are 
concerned with radar-astronomy in which 
echo techniques are used, the others with 
reception of radio waves emitted by 
extraterrestrial objects. The better chap- 
ters are those on galactic and extra- 
galactic radio emission, scintillation and 
meteors. These, rather logically, are the 
subjects with which the authors have 
been intimately concerned. The chapters 
on the 21-cm hydrogen line and the sun 
seem relatively brief compared with 
the amount of experimental work which 
has been done in these two fields. The 
chapter on the sun, for instance, contains 
only a brief reference to the work by 
Wild and his associates on the dynamic 
spectrum of solar disturbances. There 
is also no reference in the scintillation 
chapter to Wild’s classification of scintil- 
lations using the same dynamic spectro- 
graph. 

As a matter of correction, in the chapter 
on the hydrogen line and that on radio 
stars, it is mentioned that there is a 
discrepancy between the distance of the 
Cassiopeia radio-star as determined opti- 
cally and by means of the 21-cm ab- 
sorption spectrum. Since the book 
was written, Dr. Walter Baade, of the 
Mount Wilson and Palomar Observa- 
tories, has made a new determination of 
the distance based on a longer period of 
observation of the expansion of the 
visible remnant. His recent determination 
of the distance is in excellent agreement 
with the rather easier radio measurement. 


Ali Bulent Cambel, Northwestern University, Associate Edito: 


The final chapter of the book is on the 
giant 250-ft steerable telescope at Jodrel! 
Bank. In collecting area, this instrument 
is larger by a factor of 10 than any othe: 
similar instrument in operation or unde: 
construction. Its successful completior 
represents a triumph in engineering skil 
on the part of the designer and in fore- 
thought for its originators. The re- 
viewer can remember gazing down into 
the foundation pit for the central pivot 
at a time when most of the world’s radic 
telescopes would have fitted into that pit. 
Extremely valuable work has been done 
by this instrument in both radar and 
radio tracking of the Sputniks and 
American Earth satellites. 

The minor criticisms of this review 
detract only slightly from the value o 
this book. In general it is very wel 
written and is sure to appeal to a wide- 
spread audience. There are few, i! 
any, typographical errors. The upper 
diagram in Fig. 33 might amuse those 
acquainted with pen recorders. This 
surely is a case of over-idealization. 


Jet Propulsion, by Walter J. Hesse, 
Pitman Publishing Corp., New York, 
1958, 585 pp. 

Reviewed by P. Roy CHoupHURY 
University of Southern California 


Because of the increasing student 
interest in the field of jet propulsion, this 
book is welcome as an undergraduate 
text. 

During the last academic year, this 
reviewer had the difficult job of selecting 
an undergraduate textbook for a course in 
Rockets and Thermal Jets, which covered, 
among other topics, not only the funda- 
mentals of jet propulsion systems but also 
the fundamentals of gasdynamics. Because 
of the meager supply of books covering 
these subjects in one volume, the re- 
viewer feels that the present text has made 
a timely arrival in its field. Besides the 
basic fundamentals of operation of rockets 
and air-breathing engines, the book covers 
a review of thermodynamics, compressible 
flow, turbomachinery, charts of standard 
atmospheric data, tables of compressible 
isentropic flow and normal shock tables. 

“Jet Propulsion” is most appropriate as 
a text for an undergraduate course. It 
can be covered in two quarters of the 
senior year. However, the graduate 
students in the field would have had a 
major part of the material in other re- 
quired courses. Nonetheless, non-mechan- 
ical or -aeronautical graduate students 
(e.g., civil engineering, etc.) can prof- 
itably use this as a survey course for one 
quarter. 

In the review of the basic physical 
laws, the author writes the equation of 
motion with friction and calls it Bernoulli’s 
equation. According to most authors, 
Bernoulli’s equation is an integrated form 
of Euler’s equation and does not contain 
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AERO-THERMODYNAMICISTS 
EXPLORE HIGH-SPEED RE-ENTRY 


requiring attention to interaction between aerodynamic and thermodynamic 

mI phenomena. Typical of these is the problem of high-speed atmospheric re-entry. 


A report to Engineers = er Expanding research and development activities have coincided with acceleration 
on top priority programs like our Polaris IRBM. At the same time, positions 


4) 7s Advanced weapon system technology has brought to the forefront problem areas 


and Scientists from ® for qualified engineers and scientists have opened up that are unequalled in +e 
Lockheed Missile Systems —_ responsibility or in opportunities for moving ahead. ‘ 
where expanding missile ‘ Positions in aero-thermodynamics include such areas as: aerodynamic ga e 
: characteristics of missiles at high Mach numbers; missile and weapon system 8 
programs insure more design analysis; boundary layer and heat transfer analyses in hypersonic flow 7 
promising careers pas fields; and calculation of transient structural and equipment temperatures a4 
: resulting from aerodynamic heating and radiation. 


— a In addition, openings exist at all levels in Gas Dynamics, Structures, Propulsion, 
Test Planning and Analysis, Test Operations, Information Processing, Electronics, 
and Systems Integration. For these and other positions, qualified engineers 

and scientists are invited to write Research and Development Staff, 


- Dept. 2512, 962 W. El Camino Real, Sunnyvale, California. 


 kLochheed 
¥ / 


SUNNYVALE, PALO ALTO, VAN NUYS, SANTA CRUZ, SANTA MARIA, CALIFORNIA 
CAPE CANAVERAL, FLORIDA ¢ ALAMOGORDO. NEW MEXICO — 


MISSILE SYSTEMS DIVISION 


Maurice Tucker, Aero-T hermodynamics Department Manager, right, 
discusses combined aero-thermodynamic re-entry body tests being 
conducted in Division's new “hot-shot” wind tunnel. Others are 

Dr. Jerome L. Fox, Assistant Department Manager, Thermodynamics, 
and Robert L. Nelson, Assistant Department Manager, Aerodynamics. 
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any term due to friction. Perhaps, for 
the benefit of the students, it would have 
been better to conform to the common 
terminology. In another instance, in 
order to illustrate that in a cyclic process, 
energy in the form of heat can not be 
completely converted to continuous useful 
work, the author states “.. .heat energy 
is a low grade type energy... .’’ Although 
this statement, in the author’s mind, may 
describe the situation adequately, its 
place in a textbook of this sort is not very 
appropriate. 

In the section on afterburners, a brief 
discussion of flame stabilization by dif- 
ferent types of flameholders would not 
have been out of place. Lastly, in view 
of the ever-increasing importance of 
solid propellants, the chapter on rockets 
should definitely have included a discus- 
sion of the internal ballistics of rockets 
of this type. 


Aircraft and Missile Propulsion, Volume I, 
Thermodynamics of Fluid Flow and 
Application to Propulsion Engines, 
by M. J. Zucrow, John Wiley & Sons, 
New York, 1958, xiv + 538 pp. $11.50. 

Reviewed by A. E. Funs 
Northwestern University 


Ten years ago, the first edition of 
Professor Zucrow’s book “Principles of 
Jet Propulsion and Gas Turbines’ was 
published. In the intervening period 
there has been an immense growth in the 
field of jet propulsion. The second edi- 
tion, also, has grown and comprises three 
volumes, the first of which is reviewed 
here. Volume II deals with the cycle 
analysis and performance calculations of 
ramjets, rocket engines and the gas 
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Results of Static Test Experiment with 
Hot Water Rocket Model, by E. Shafer 
and W. Michely, Forschungsinstitut fur 
Physik der Strahlantriebe E. V. Mitteilun- 
gen 11, Aug. 1957, 47 pp. (in German). 

A Rocket Drive for Long Range Bomb- 
ers, by Eugen Sanger and Irene Bredt, 
Forschungsinstitut fur Physik der Strahl- 
antriebe E. V. Mitteilungen 13, Oct. 1957, 
266 pp. (Reprint of ZWB Untersuchun- 
gen und Mitteilungen 3538, Aug. 1944.) 
(In German.) 

Propulsion Systems for Space Flight, 
by R. B. Dillaway, Aeron. Engng. Rev., 
vol. 17, April 1958, pp. 42-49, 52. 

Rocket Pumps Have Reached High 
Efficiency, by Kurt R. Stehling, Aviation 

Age, vol. 29, pp. 32-33, 35, 37, 41-42. 

Taming the Supersonic Turbojet, by 

Edward A. Simonis, Aviation Age, vol. 29, 
April 1958, pp. 60-69. 


_ Eprror’s Note: Contributions from Pro- 


fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine, Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
_ Minnesota, are gratefully acknowledged. 


turbine power plant as applied to the 
turboprop and turbojet. Volume III 
discusses the components of these engines. 

In the preface the author states that 
the prime objective of the second edition 
is to furnish the student and the practicing 
engineer with an understanding of the 
fundamental principles governing the 
functioning and operating characteristics 
of the engines employed for propelling 
high speed aircraft and missiles. To 
achieve this objective, the five chapters of 
Volume I are devoted to a review of the 
fundamental principles, the general charac- 
teristics of propulsion systems, the ther- 
modynamics of compressible fluid flow, 
flow through nozzles and flow through 
diffusers. Throughout Volume I, the 
author follows the policy of deriving 


general results and then eliminating 
those terms irrelevant to a_ specific 
application. 


Compared to the first edition, this book 
has considerably more examples with 
detailed solutions integrated into the 
text. The engineer who wishes to enhance 
his knowledge through self study and the 
professor who may want to use Volume I 
as a classroom text will be pleased to 
know that there are a generous number of 
exercises. All the tables have been 
grouped in the appendix. Some of the 
notation has also been changed, adding 
to the clarity of the presentation. The 
extensive lists of references, for which the 
first edition was noted, have been brought 
up to date. 

Since Volume I deals mainly with the 
fundamentals which are as valid today as 
they were when originally discovered by 
Newton, Euler, Mach, Rayleigh and 
others, it is not possible to evaluate how 


successfully the author has “modernized”’ 
the text. It is apparent, however, that 
the author’s viewpoint includes propulsion 
as applied to missiles as well as aircraft. 

The value of this book arises not so 
much from the precision with which the 
fundamental facts concerning thermo- 
dynamics and fluid mechanics have becn 
stated, but from the completeness of tlie 
material included and the variety of 
applications treated. This book, like 
the first edition, will be a useful referen: e 
for engineers working in the field, «s 
well as a good introductory text fir 
newcomers to jet propulsion 


Books Received | 


Separation and Purification of M:- 
terials, by R. Hammond, Philosophic::| 
Library, New York, 1958, 327 pp. $10. 

Dispersion of Materials, by R. Han.- 
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Report on the Ex- 
perimental Investigation of the RPI Wave 
Engine, by H. J. Lopez and R. E. Jensen, 


Rensselaer Polytech. Inst., Dept. Aeron. 
Engng., Tech. Rep. AE 5801, Feb. 1958, 
23 pp. 

Booster Propulsion for Space Vehicles, 
by R. S. Wentink, Inst. Aeron. Sci., Pre- 
print 828, Jan. 1958, 29 pp. 2 figs. 

An Optimum Design for a Semi In- 
finite Rocket Wall Containing a Circum- 
ferential Keyway, by Bernard W. Shaf- 
fer, Ira Cochin and Morton Mantus, /nst. 
Aeron. Sci., Preprint 767, Jan. 1958, 26 
pp., 9 figs. 

Photographic Study of Rotary Scream- 
ing and Other Oscillations in a Rocket 
Engine, by Theodore Male, William R. 
Kerslake and Adelbert O. Tischler, NACA 
Res. Mem. E54A29, May 1954, 37 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstr.1 25, p. 20, 
3/18/58.) 

Propellant Vaporization as a Criterion 
for Rocket Engine Design; Relation Be- 
tween Percentage of Propellant Vaporized 
and Engine Performance, by Marcus F. 
Heidmann and Richard J. Priem, N ACA 
TN 4219, March 1958, 19 pp. 

Application of a High-temperature Static 


Strain Gage to the Measurement of Ther- 
mal Stresses in a Turbine Stator Vane, by 
R. H. Kemp, C. R. Morse and M. H. 
Hirschberg, NACA TN 4215, March 
1958, 36 pp. 

High-altitude Performance of 9.5-inch- 
diameter Tubular Experimental Combus- 
tor with Fuel Staging, by Wilfred E. Scull, 
NACA Res. Mem. £54A06, March 1954, 
55 pp. (Declassified from Confidential by 
authority of NACA Res. Abstr. 125, p. 
20, 3/18/58.) 

Component and Over-all Performance 
Evaluation of an Axial-flow Turbojet En- 
gine over a Range of Engine-inlet Rey- 
nolds Numbers, by Curtis L. Walker, 8S. C. 
Huntley and W. M. Braithwaite, NACA 
Res. Mem. ¥52B08, July 1952, 42 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstr. 125, p. 19, 
3/18/58.) 

Performance of an Impulse-type Super- 
sonic Compressor with Stators, by John F. 
Klapproth, Guy N. Ullman and Edward 
R. Tysl, NACA Res. Mem. E52B22, 
April 1952, 22 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 125, p. 19, 3/18/58.) 
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Du Pont hme to 


When you ignite a rocket propellant . 


(1) Speed rocket-fuelignition 
(2) Provide remote detonations 


. actuate a valve ... release 


a switch... you often rest the success or failure of the flight on some 
form of detonation. Now Du Pont research brings you two new prod- 
ucts that reduce the chances of failure almost to the vanishing point 


. and give you faster and better controlled detonation. 
aye 


DU PONT “MDF” MILD 7 
DETONATING FUSE 


. carries non-destructive detonation 
safely past delicate instruments to 
perform remote functions 

“MDF” is a small diameter flexible 
metal tubing, available in practically 
any length you desire, and contain- 
ing a detonating composition as its 


core. 


SEALING 


USES OF “MDF” include: shearing or severing 
covers or diaphragms; activating various 
valves, releases, arming devices, etc.; simul- 
taneously initiating multiple explosive charges. 

Five strengths are available (1, 2, 
5, 10 and 20 grains of PETN per 
foot), providing a mild explosive 
power sufficient to carry a detonation 
from one point to another, but usu- 
ally not enough to ignite or initiate 
other explosives or to burst violently 
out of its tubing. 

Thus “MDF” may be used to 
transfer a detonation through or 
around delicate instruments. ; 


Sample Kits Available « Sample kits con- 
taining either ‘“‘MDF”’ or ‘‘Pyrocore,”’ plus 
various primers and accessories, are avail- 
able from Du Pont. Kit prices and technical 
bulletins can be obtained by writing to E. I. 
du Pont de Nemours & Co. (Inc.), Explosives 
Department, Wilmington 98, Delaware. 


Strength, velocity and tempera- 
ture resistance can be varied accord- 
ing to your needs, and other types 
of tubing can be provided, including 
plastics, wires and textiles to over- 
come problems of abrasion, corro- 
sion, damaging effect and noise con- 
trol. 

You can plan precise detonation 
speeds with “MDF” because velocity 
for a given core loading is accurate 
within + 1%. 

And you can depend on each det- 
onation coming off as planned. In 
recent tests, thousands and _ thou- 
sands of feet of ““MDF’” were fired, 
without a single failure. Approach- 
ing 200,000 feet, the tests were dis- 
continued because they failed to pro- 
duce any failures on which to base 
reliability figures. 

Handling of “MDF”’ is safe and 
convenient compared to ordinary 
commercial explosives, because only 
small quantities of explosives are in- 
volved and even these are insensi- 
tive to ordinary shock. The tubing 
is quite flexible and lends itself eas- 
ily to stringing, wrapping, threading. 


DU PONT “PYROCORE” IGNITER 


. provides faster, more uniform igni- 
tion of rocket propellants 


“Pyrocore,” too, is a flexible metal 


EXPLOSIVES DEPARTMENT 


Better Things for Better Living . 


tubing available in any length. Its 
core is a detonating-ignition compo- 
sition. 

Propellant ignition may be ef- 
fected along its length at velocities 
ranging from 12,000 to 21,000 feet 
per second, depending upon the type 
and amount of core composition you 
specify. 

In laboratory tests, ““Pyrocore” cut 
ignition time 99.5% when compared 
with a standard primer. Total igni- 
tion of a 2114” length of cannon 
primer was achieved in less than 14 
millisecond, as against 50 millisec- 
onds required with a squib primer. 


TYPICAL APPLICATIONS OF “PYROCORE” 
IGNITER. You can string, thread, weave it 
through or around the propellant. 


The stringlike form of “Pyrocore” 
gives you new freedom in designing 
rocket motors. For example, you can 
thread “Pyrocore” through jelly roll 
or basket type igniter assemblies, or 
string it around the propellant 
grains. Or you can weave coils or 
loops of ‘““Pyrocore” directly into the 
propellant grain either during or 
after grain formation. 

Further, ‘““Pyrocore” has very low 
brisance (shattering effect) and can 
convey a non-electric stimulus safe- 
ly past expensive equipment and 
bulkheads. 

Initiator assemblies may be lo- 
cated outside the propellant cham- 
ber, eliminating a major source of 
residue in the reaction zone. The 
“‘Pyrocore” itself leaves very little 
residue. The portion of “Pyrocore” 
not actually used for ignition may be 
easily insulated with any one of a 
variety of plastic compounds. 
“Pyrocore” is a Du Pont trade mark. 
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sonic Impeller with Axial Discharge, by 
Linwood C. Wright and Karl Kovach, 
NACA Res. Mem. E53B09, April 1953, 
37 pp. (Declassified from Confidential by 
authority of NACA Res. Abstr. 125, p. 
20, 3/18/58.) 

Altitude Performance Investigation of 
Two Flame-holder and Fuel-system Con- 
figurations in Short Afterburner, by S. C. 
Huntley and H. D. Wilsted, NACA Res. 
Mem. E52B25, May 1952, 41 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 125, p. 19, 3/18/58.) 


Altitude Starting Characteristics of an 
Afterburner with Autoignition and Hot- 
streak Ignition, by P. E. Renas, R. W. 
Harvey Sr. and E. T. Jansen, NACA Res. 
Mem. E53B02, April 1953, 25 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 125, p. 19, 3/18/58.) 


Turboprop-engine Design Considera- 
tions, I: Effect of Mode of Engine Opera- 
tion on Performance of Turboprop Engine 
with Current Compressor Pressure Ratio 
by Elmer H. Davison, NACA Res. Mem. 
E54D19, May 1955, 34 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 125, p. 21, 3/18/58.) 


Turboprop-engine Design Considera- 
tions, II: Design Requirements and Per- 
formance of Turboprop Engines with a 
Single-spool High-pressure-ratio Com- 
pressor, by Elmer H. Davison and Mar- 
garet C. Stalla, NACA Res. Mem. E55B18, 
May 1955, 32 pp. (Declassified from Con- 
fidential by authority of NACA Res. 
Abstr. 125, p. 22, 3/18/58.) 


Correlation of Turbine-blade-element 
Losses Based on Wake Momentum Thick- 
ness with Diffusion Parameter for a 
Series of Subsonic Turbine Blades in Two- 
dimensional Cascade and for Four Tran- 
sonic Turbine Rotors, by Robert Y. Wong 
and Warner L. Stewart, NACA Res. 
Mem. E55B08, April 1955, 31 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 125, p. 22, 3/18/58.) 

Photographic Investigation of Air-flow 
Patterns in Transparent One-sixth Sector 
of Annular Turbojet-engine Combustor 
with Axial-slot-type Air Admission, by 
Charles C. Graves and J. Dean Gernon, 
NACA Res. Mem. E54128a, Dec. 1954, 
24 pp. (Declassified from Confidential by 
authority of NACA. Res. Abstr. 125, p. 
21, 3/18/58.) 

Effect of Nozzle Secondary Flows on 
Turbine Performance as Indicated by 
Exit Surveys of a Rotor, by Warren J. 
Whitney, evend A. Buckner Jr. and 
Daniel E. Monroe, NACA Res. Mem. 
E54B03, April 1954, 11 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 125, p. 20, 3/18/58.) 

The Transonic Flow Field of an Axial 
Compressor Blade Row, by James E. 
McCune, Inst. Aeron. Sci., Preprint 792, 
Jan. 1958, 26 pp., 6 figs. 

Use of Effective Momentum Thickness 
in Describing Turbine Rotor-blade Losses, 
by Warner L. Stewart, Warren J. Whitney 


and James W. Miser, NACA Res. Mem. 


_ E56B29, May 1956, 26 pp. 


(Declassified 


_ from Confidential by authority of NACA 
_ Res. Abstr. 125, p. 23, 3/18/58.) 


Engine Afterburners, by 


Low-velocity Turning as a Means of 
Minimizing Boundary-layer Accumula- 
tions Resulting from Secondary Flows 
within Turbine Stators, by Warner L. 
Stewart and Robert Y. Wong, NACA Res. 
Mem. E54B16, May 1954, 18 pp. (De- 
classified from Confidential by authority 
of NACA Res. Abstr. 125, p. 21, 3/18/58.) 

Performance Characteristics of Several 


Short Annular Diffusers for Turbojet 
William E. 


Mallett and James L. Harp Jr., NACA 
Res. Mem. E54B09, May 1954, 31 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstr. 125, p. 21, 
3/18/58.) 

Effect of Ambient Conditions on the 
Performance of a Pressure-jet Power- 
plant for a Helicopter, by Richard P. 
Krebs, NACA Res. Mem. E56B21, June 
1956, 49 pp. (Declassified from Confi- 
dential by authority of NACA Res. Abstr. 
125, p. 22, 3/18/58.) 

Exploratory Investigation of a Helicop- 
ter Pressure-jet System on the Langley 
Helicopter Test Tower, by Robert A. 
Makofski and James P. Shivers, NACA 
Res. Mem. L56B17, July 1956, 38 pp. 
(Declassified from Confidential by author- 


ity of NACA Res. Abstr. 125, p. 26, 3/18/ — 
58.) 


Ram-jets, by R. P. Probert, J. Roy. 
Aeron. Soc., vol. 62, March 1958, pp. 151— 
173. 

Altitude Performance Investigation of 
Two Single-annular Type Combustors 
and the Prototype J40-WE-8 Turbojet 
Engine Combustor with Various Combus- 
tor Inlet-air-pressure Profiles, by Adam 
E. Sobolewski, Robert R. Miller and John 
E. McAulay, NACA Res. Mem. E52J07, 
May 1953, 46 pp. (Declassified from Con- 
fidential by authority of NACA Res. Abstr. 
124, pp. 10, 2/10/58.) 

Compressor Turbine Matching of Two- 
spool Turbo-props, by A. W. Morley, 
Aeron. Quart., vol. 9, Feb. 1958, pp. 17-33. 

The High-density Turboprop: The 
Vanguard Examined, by Roy Allen, 
Aeronautics, vol. 38, March 1958, pp. 47- 
51. 

Propulsion, by Jerry Grey, Aviation 
Age, vol. 28, March 1958, pp. 36-43. 

Acoustic, Thrust, and Drag Characteris- 
tics of Several Full-scale Noise Suppres- 
sors for Turbojet Engines, by Carl C. 
Ciepluch, Warren J. North, Willard D. 
Coles and Robert J. Antl, NACA TN 
4261, April 1958, 48 pp. 

Transonic Drag of Several Jet-noise 
Suppressors, by Warren J. North, NACA 
TN 4269, April 1958, 34 pp. 

Discharge Coefficients for Combustor- 
liner Air-entry holes, II: Flush Rectan- 
gular Holes, Step Louvers, and Scoops, 
by Ralph T. Dittrich, NACA TN 3924, 
April 1958, 56 pp. 

Effect of Prior Air Force Overtempera- 
ture Operation on Life of J47 Buckets 
Evaluated in a Sea-level Cyclic Engine 
Test, by Robert A. Signorelli, James R. 
Johnston and Floyd B. Barrett, NACA 
TN 4263, April 1958, 41 pp. 

Maximum Theoretical Tangential Ve- 
locity Component Possible from Straight- 
back Converging and Converging-diverg- 
ing Stators at Super-critical Pressure 
Ratios, by Thomas P. Moffitt, VACA TN 
4271, April 1958, 21 pp. 

Equivalent Performance Parameters for 
Turboblowers and Compressors, by Hunt 
Davis, Trans. ASME, vol. 80, Jan. 1958, 
pp. 108-116. 


A Method for Calculating Vibration Fre- 
quency and Stress of a Banded Group of 
Turbine Buckets, by M. A. Prohl, 7rans, 
ASME, vol. 80, Jan. 1958, pp. 169-180, 


High-frequency Vibration of Steam- 
turbine Buckets, by F. L. Weaver and M. 
A. Probl, Trans. ASME, vol. 80, Jan. 
1958, pp. 181-194. 


F-104 Intake Ducts Save Weight and 


Complexity, by A. F. Watts, Aviation Age, 

vol. 29, May 1958, pp. 60-67. 
Testing Paper Gas Turbines: A 

the Digital Computer to Engine 


lying 
esign, 


by Paul T. Vickers and Charles A. Amann, 
General Motors Engng. J., vol. 5, April- 
June, 1958, pp. 6-16. 

Boost Systems for Helicopter Gas Tur- 


bines, by A. W. Morley, J. Helicopter 
Assoc. of Gt. Brit., vol. 12, April 1958, pp. 
66-83; Discussion, pp. 84-90. 

A Combination Powerplant for High- 
performance Aircraft, by J. Dupin, Jn- 
teravia, vol. 13, April 1958, pp. 349-351. 

Measurement of Ignition Delay at 
Simulated Altitude of Hypergolic Liquid 
Rocket Motors, by Pierre Sarrat, Lz 
Recherche Aeronautique, No. 62, Jan.-Feb. 
1958, pp. 15-25 (in French). 


Aerodynamics of Jet 
Propelled Vehicles 


Experimental Study of the Equivalenc: 
of Transonic Flow about Slender Cone- 
cylinders of Circular and Elliptic Cross 
Section, by William A. Page, NACA TA 
4233, April 1958, 45 pp. 

A Second-order Shock-expansion 
Method Applicable to Bodies of Revolu- 
tion Near Zero Lift, by Clarence A. Sy 
vertson and David H. Dennis, NACA 
Rep. 1328, 1957, 20 pp. (Supersedes Res 
Mem. L55G05.) 

Stagnation-point Heat Transfer to Blunt 
Shapes in Hypersonic Flight, Including 
Effects of Yaw, by A. J. Eggers Jr., C 
Frederick Hansen and Bernard E. Cunning- 
ham, NACA T'N 4229, April 1958, 54 pp. 

Stationary Wall Temperature of a Body 
under the Influence of Aerodynamic 
Heating Flying at Mach Numbers from 
1 to 10 and up to 30 km, by H. G. L. 
Krause and M. E. Kubler, Forschungsin- 
stitut fur Physik der Strahlantreibe E. V., 
Mitteilungen 7, Oct. 1956, 138 pp. (in 
German). 

The Structure of a Steady State One 
Dimensional Detonation Wave Supported 
by a Reversible Reaction, by B. Linder, 
C. F. Curtiss, and J. O. Hirschfelder, 
Univ. of Wisconsin, Naval Res. Lab., 
CM-91la, Dec. 1957, 11 pp. 

A Discussion of Higher-order Approxi- 
mations for the Flow Field about a Slender 
Elliptic Cone, by Roberto Vaglio-Laurin 
and M. D. Van Dyke, J. Fluid Mech., vol. 
3, March 1958, pp. 638-644. 

Hypersonic Flight and the Re-entry 
Problem (Twenty-first Wright Brothers 
Lecture), by H. Julian Allen, J. Aeron. 
Sci., vol. 25, April 1958, pp. 217, 229, 262. 

On Optimum Nose Curves for Super- 
aerodynamic Missiles, by H. S. Tan, J. 
Aeron. Sci., vol. 25, April 1958, pp. 263- 
264. 

Drag Due to Lift of a Not So Slender 
Configuration, by Hsien K. Cheng, Wright 
Air Dev. Center, Tech. Rep. 57-316, Part I. 
(ASTIA AD 118336), Sept. 1957, 38 pp. 

The Supersonic Blunt Body Problem— 
Review and Extension, by Milton D. Van 
Dyke, Inst. Aeron. Sci., Preprint 801, Jan. 
1958, 19 pp., 18 figs. 

Examples of Detached Bow Shock 
Waves in Hypersonic Flow, by P. R. 
Garabedian and H. M. Lieberstein, Jnst. 
Aeron. Sci., Preprint 817. 

The Dynamics and Certain Aspects of 
Control of a Body Re-entering the Atmos- 
phere at High Speed, by Joseph D. Welch 
and S. L. Shih, Jnst. Aeron. Sci., Preprint 
818, Jan. 1958, 12 pp., 12 figs. 


Experimental Pressure Distribution on 
an ymmetrical Nonconical Body at 
Mach Number 1.90, by De Marquis D. 
Wyatt, NACA Res. Mem. E9B03, Feb. 
1949, 58 pp. (Declassified from Confiden- 
tial by authority of NACA Res. Abstr. 
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SENIOR SCIENTISTS AND ENGINEERS 
AVCO—PIONEER IN RE-ENTRY—IS EXPLORING 
NEW APPROACHES TO SPACE AND MISSILE TECHNOLOGY | 


The Aveo Research and Advanced Development Division is 
conducting an extensive program of basic and applied research 
in the physics, chemistry and engineering associated with space 
and missile technology. Supervisory and staff positions are avail- 
able for creative senior physicists and engineers—both theo- 
retical and experimental—and for electrical engineers with a 
strong physics background. 
Unusual and challenging openings exist in the following fields: 


High-intensity arcs, gaseous discharge phenomena, properties of 
gases at high temperatures, radiation measurements 


High-temperature reflectivity and emissivity measurements 
High-temperature properties of materials 
Infrared 
Terminal ballistics of high-velocity particles 
Chemical physics, spectroscopy, surface physics, hydrodynamics, fluid 
dynamics, free molecule flow studies, upper atmosphere phenomena 
Missile detection and discrimination, advanced missile warfare 
concepts 
Ultra-high-speed electronic and optical instrumentation 
Microwaves, telemetry systems, radar, propagation through ionized 
plasmas, space communication studies 
The division’s new suburban location provides an unusually 
attractive working environment outside of metropolitan Boston. 
The large, fully equipped, modern laboratory is in pleasant 
surroundings, yet close to Boston educational institutions and 
cultural events. Publications and professional development are 
encouraged, and the division offers a liberal educational assistance 
Dr. R. W. Johnston, 


program for advanced study. _ 
Scientific and Technical Relations, 


Avco Research and Advanced Development Div., 
201 Lowell Street, Wilmington, Massachusetts 


Mesearch and Atvanced Development, 


Address all inquiries to: 
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DECONTAMINATION 
BOOTH 


Instantly ready for 
vital FIRST AID 


Miscues and accidental exposure to dan- 
gerous propellants and other chemicals can 
occur with shocking suddenness. Adequate 
water irrigation is an important key to min- 
imizing such injuries and subsequent 
claims. HAWS Decontamination Booth is 
the answer...a complete safety station for 


immediate first aid. 4 
HAWS MODEL 8600 
DECONTAMINATION BOOTH 

is made of durable, lightweight reinforced 
fiberglass plastic, and features Haws Eye- 
Face Wash Fountain, eight lateral body sprays 
and overhead spray unit. All are simultane- 
ously activated by weight on the base- 
mounted foot treadie! Contaminated victims 
are instantly “covered” with water that 
floats away foreign matter from body and 
clothing. 

At aeronautical and astronautical installa- 
tions everywhere, HAWS Safety Facilities are 
important in boosting air-age safety pro- 
grams. Find out what this equipment can 


mean to your operation. Full details sent on 
request, with no obligation. 


DRINKING FAUCET CO. 


EXPORT DEPARTMENT * 19 COLUMBUS STREET 


SAN FRANCISCO 11, CALIFORNIA 


856 


125, p. 18, 3/18/58.) 

Effect of Mach Number on Boundary- 
layer Transition in the Presence of Pres- 
sure Rise and Surface Roughness on an 
Ogive-cylinder Body with Cold Wall 
Conditions, by Robert J. Carros, NACA 
Res. Mem. A56B15, April 1956, 30 pp. 
(Declassified from Confidential by author- 
ity of NACA Res. Abstr. 125, p. 18, 3/- 
18/58.) 

Transonic Free-flight Drag Results of 
Full-scale Models of 16-inch Diameter 
Ram-jet Engines, by Wesley E. Messing 
and Loren W. Acker, NACA Res. Mem. 
52B19, April 1952, 17 pp. (Declassified 
from Confidential by authority of NACA 
Res. Abstr. 125, p. 19, 3/18/58.) 


On Optimum Nose Shapes for Missiles 
in the Superaerdyanmic Region, by I. D. 
Chang, J. Aeron. Sci., vol. 25, Jan. 1958, 
pp. 57-58. 

On Hypersonic Blunt Body Flow with a 
Magnetic Field, by Nelson H. Kemp. 
Avco Mfg. Corp., Avco Res. Lab., Res. Rep. 
19, Feb. 1958, 9 pp., 2 figs. 

Lift and Drag on Cone Cylinders, by T. 
Nark, Univ. of Calif., Inst. Engng. Res., 
Tech. Rep. HE-150-154, Dec. 1957, 32 pp., 
10 figs. 

Pressure and Wave Drag Coefficients 
for Hemispheres, Hemisphere Cones, and 
Hemisphere Ogives, by Jerome R. Katz, 
NAVORD Rep. 5849 (NOTS 1947), 
March 1958, 17 pp. 


Effect of Jet-Nozzle-expansion Ratioon 


Drag of Parabolic Afterbodies, by Gerald 
W. Englert, Donald J. Vargo and Robert 


Heat Transfer and 


W. Cubbison, NACA Res. Mem. E54B12, 


April 1954, 26 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 125, p. 21, 3/18/58.) 

Aerodynamic Characteristics of a Cru- 
ciform-wing Missile with Canard Control 
Surfaces and of Some Very Small Span 
Wing-body Missiles at a Mach Number 
of 1.41, by M. Leroy Spearman and Ross 
B. Robinson, NACA Res. Mem. L54B11, 
April 1954, 27 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 125, p. 24, 3/18/58.) 

A Linear Perturbation Method for 
Stability and Flutter Calculations on Hy- 
personic Bodies, by Maurice Holt, Jnst. 
Aeron. Sci., Preprint 793, Jan. 1958, 20 
pp. 

Similitude of Hypersonic Flows over 
Thin and Slender Bodies—an Extension to 
Real Gases, by Hsien K. Cheng, Cornell 
Aeron. Lab., Inc., Rep. AD-1052-A-6 
(AFOSR-TN. -58-87: ASTIA A D 148136), 
Feb. 1958, 13 pp., 6 figs. 

Preliminary Investigations of Spiked 
Bodies at Hypersonic Speeds, by 8S. M. 
Bogdonoff and I. E. Vas, Princeton Univ. 
Dept. of Aeron. Engng., Rep. 412, (Wright 
Air Dev. Center, TN 58-7; ASTIA AD 
142280), March 1958, 10 pp., 16 figs. 

Molecular Approach to Problems of 
High Altitude, High Speed Flight, by G. 
N. Patterson, NATO, Advis. Group for 
Aeron. Res. and Dev., Rep. 134, July 1957, 
58 pp. 

The Direct Measurement of Local Skin 
Friction on Viking No. 13 During Ascend- 
ing Flight, by Felix W. Fenter and W. C. 
Lyons Jr., Univ. of Texas, Defense Res. 
Lab., CM 917 (DRL-414), Jan. 1958, 32 
pp., 4 tabs., 16 figs. 

On the Drag of a Sphere at Extremely 
High Speeds, by V. C. Liu, J. Appl 
Phys., vol. 29, Feb. 1958, pp. 194-195. 

A Model of Supersonic Flow Past Blunt 
Axisymmetric Bodies, with Applications to 
Chester’s Solution, by M.D. Van Dyke, 
J. Fluid Mech., vol. 3, Feb. 1958, pp. 515- 
526. 

An Approximate Method of Determining 
Axisymmetric Inviscid Supersonic Flow 
over a Solid Body and Its Wake, by Sin-I 
Cheng, J. Aeron. Sci., vol. 25, March 1958, 
pp. 185-193. 

Approximate Solution for Slightly Yawed 
n-Power Bodies at Hypersonic Speeds, by 
Robert W. Truitt, J. Aeron. Sci., vol. 25, 
March 1958, pp. 206-207. 


On Optimum Nose Shapes for Missiles 
in the Superaerodynamic Region, by 
David H. Dennis, J. Aeron. Sci., vol. 25, . 
March 1958, p. 216. 


Boundary Layer Growth on a Spinning 
Body: Accelerated Motion, by Y. D. 
Wadhwa, Phil. same vol. 3, Feb. 1958, PP. 
152-158. 


Fluid Flow 


Crocco’s Vorticity Law in a Non-uni- 
form Material, by C. S. Wu and W. D. 
Hayes, Quart. Appl. Math., vol. 16, April 
1958, pp. 81-82. 

Temperature Dependence of Vibrational 
Relaxation Times in Gases, by P. G. 
Corran, J. D. Lambert, R. Salter and B. 
Warburton, Proc. Roy. Soc., vol. 244A, 
no. 1237, March 11, 1958. pp. 212-219. 


The Determination of the Lorenz Num- 
ber at High Temperatures, by M. R. 
Hopkins and R. L. Griffith, Zeitschrift fir 
Physik, vol. 150, no. 3, Feb. 14, 1958, pp. 
325-331 (in English). 

Aerodynamic Research on Accelerating 
Cascade; Cascade Test of the Nozzle 
Blades, by Isamu Wada, Proc., Japan, 
6th Natl. Congr. Appl. Mech., 1956 (pub. 
1957), pp. 333-336. 

Design of Cascade in Compressible 
Fluid and Application, by Ichiro Hama- 
moto, Proc., Japan, 6th Natl. Congr. A ppl. 
Mech., 1956 (pub. 1957), pp. 337-340. 

On the Evaporation, Ignition Lags and 
Combustion of Fuel Droplets, by Niichi 
Nishiwaki, Masaru Hirata, Sanni Hagi 
and Shin-ichiro Yamazaki, Proc., Japan, 
6th Natl. Congr. Appl. Mech., 1956 (pub. 
1957), pp. 407-410. 

Thermodynamics, by R. A. Budenhol- 
zer and Alfred Ritter, Aviation Age, vol. 
28, March 1958, pp. 44-49. 

Fluid Dynamics, by A. K. Oppenheim, 
C. L. Coldren, L. M. Grossman and C. V. 
Sternling, Ind. and Engng. Chem., vol. 50 
March 1958, Pt. 2, pp. 525-542. 

Heat Transfer, by E. R. G. Eckert, J. P. 
Hartnett and T. F. Irvine Jr., Ind. and 
Engng. Chem., vol. 50, March 1958, Pt. 2, 
pp. 543-554. 

Mass Transfer, by C. R. Wilke and J. 
M. Prausnitz, Ind. and Engng. Chem., vol. 
50, March 1958, Pt. 2, pp. 555-560. 


Thermodynamics, by J. M. Smith and 
G. M. Brown, Ind. and Engng. Chem., vol. 
50, March 1958, Pt. 2, pp. 561-568. 

Non-linear Interactions in a Viscous 
Heat-conducting Compressible Gas, by 
Boa-Teh Chu and Leslie 8. G. Kovasznay, 
J. Fluid Mech., vol. 3, Feb. 1958, pp. 494— 
504. 


Pyrolysis of Simple Hydrocarbons in 
Shock Waves, by E. F. Greene, R. L. 
Taylor and W. L. Patterson Jr., J. Phys. 
Chem., vol. 62, Feb. 1958, pp. 238-243. 


Incompressible Two-dimensional Stag- 
nation-point Flow of an Electrically Con- 
ducting Viscous Fluid in the Presence of a 
Magnetic Field, by Joseph L. Neuringer 
and William Mcllroy, J. Aeron. Sct., vol. 
25, March 1958, pp. 194-198. 


(Continued on page 859) 
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Pennsalt announces completion of its first facility This is another contribution to the missile age by 
at its Portland, Oregon plant for the production Pennsalt, a pioneer in the development of ele- 
of AMMONIUM PERCHLORATE, a leading oxidizer mental fluorine and other high energy fluorine 
used in solid propellants. Sodium chlorate, neces- chemicals, and a basic producer of quality 
sary for AMMONIUM PERCHLORATE manufacture, chemicals for more than one hundred years. We 
has been produced by Pennsalt at this location invite you to call or write Pennsalt of Washington 


for over 17 years. Division for detailed information. 


PENNSALT OF WASHINGTON DIVISION 


P NNSALT CHEMICALS CORPORATION 


TACOMA 1, WASHINGTON 
PLANTS AND OFFICES: 


Los Angeles and Menlo Park, Calif. . Philadelphia, Pa. 


DECEMBER 1958 
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Assignments now open include... 


SYSTEMS ENGINEER to develop complex devices 
in fields of servo-mechanisms, radar, or com- 
puters, and integrate these elements in weapons 
systems. Must have 3 to 5 years’ experience in 
such activity and two years’ experience in over-all 
systems analysis. Assignment involves design 
and analysis of closed-loop systems, consisting 
of inertial and radar equipment, display mate- 
rials, and digital or analog computers. 
Qualifications: B.S. or Advanced Degree in E.E. 
or A.E. 


INERTIAL GUIDANCE ENGINEER to assume 
broad project leadership in the planning and con- 
trol of development projects. Must have 3 to 5 
years’ experience in servo-mechanisms or de- 
velopment of complex devices for military appli- 
cations, including 2 years as technical leader of 
inertial guidance system development. Must have 
experience in astro-compass, with ability to 


analyze relationship of inertial equipment with 
bombing and navigation computer. 
Qualifications: B.S. or M.S. in E.E. or Physics. 


RADAR ENGINEER to analyze ultimate limits of 
present techniques and develop new concepts of 
providing topographical sensors for advanced air- 
borne and space systems; to design airborne 
radar pulse, microwave and deflection circuitry; 
to analyze doppler radar systems in order to de- 
termine theoretical accuracy and performance 
limitations. 

Qualifications: B.S. or Advanced Degree in E.E. 
and 3 to 5 years’ experience in radar system de- 
velopment, including display equipment and cir- 
cuits, control consoles, and doppler or search 
radar design. 

CONTROL SYSTEM ANALYST to perform physical 
and mathematical analyses needed to solve com- 
plex inertial control problems with real-time digi- 


tal computers. Applications in area of bombin 
navigation systems, missile systems, special-pi 
pose computer systems such as DDA. 
Qualifications: M.S. or Ph.D. in Physics or relat 
fields with strong math minor and up to 2 yea 
experience. 


COMPUTER ENGINEER to undertake logical 4 
sign of airborne digital computer equipme 
transistorizing computer circuitry for size 4 
weight reduction. Assignment entails: Mechani 
tion of computer circuitry and packaging; int 
preting problems for solution with the IBM 74 
Qualifications: B.S. or Advanced Degree 
Engineering or Physics and 3 to 5 years’ exp4 
ence in design of digital or analog compu 
equipment. Experience desirable in transis 
technology, computer logic, programming, ins' 
mentation, computer system evaluation, or se 
mechanism design. 
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ee ere is a creative engineering assignment for you wi 
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Are you making full use of your creative abilities . .. do you want 
to undertake challenging, important career assignments . . . are q 
you looking for career opportunities . . . then you should consider oq 
a position with IBM. In Owego, IBM’s vast wealth of computer 
technology is applied to completely integrated B-70 bombing- 
navigation systems; you'll work alongside such men as Engineer 
David R. Baldauf, who says: ‘“‘With IBM | have the opportunity to 
work on completely integrated B-70 bombing-radar-navigation 
systems. This means that my assignments have tremendous 
scope and versatility. For instance, I’m now evaluating several 
contractor proposals for a monopulse radar system, and at the 
same time I’m performing a study on attaching a unique radar q 
indicator to the presently produced bombing system. I’ve had 

rapid advancement at IBM, and there are still plenty of oppor- 

tunities for me to grow.” 


IBM is a recognized leader in the electronic systems field. Its 

products are used for military as well as commercial work. You 

will find opportunities for professional advancement at IBM, j 
where the ‘‘small-team’’ approach assures quick recognition of 

individual merit. Company benefits set standards for industry 

today, and salaries are commensurate with your abilities and 

experience. 


Ideally situated in rolling New York State countryside, Owego, 
in the Binghamton Triple City area, provides an excellent environ- 
ment for gracious, relaxed family life. Owego’s proximity to both 
New York City and the Finger Lakes region offers a pleasant 
variety of recreational opportunities. 


CAREER OPPORTUNITIES IN THESE AREAS... 


e Airborne digital & analog computers 7 
e Ground support equipment 

e Inertial guidance & missile systems 
e Information & network theory 
Magnetic engineering 
Maintainability engineering 

e Optics 

e Radar electronics & systems 

e Servo-mechanism design & analysis 
e Theoretical design & analysis 

e Transistor circuits 


Qualifications: B.S., M.S., or Ph.D. in Electrical or Mechanical 
Engineering, Physics, Mathematics—and proven ability to assume 
a high degree of technical responsibility in your sphere of interest. 


There are other openings in related fields. 


FOR DETAILS, write, outlining your background and interest, to: 
Mr. P. E. Strohm, Dept. 572-Z 

International Business Machines Corporation 
Owego, New York 


| | 3 M 
® 
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MILITARY PRODUCTS 


Flows in Partly Dissociated Gases, by 
Hans-Joachim Metzdorf, J. Aeron. Sci., 
vol. 25, March 1958, pp. 200-201. 

Experiments on Chemical Kinetics in a 
Supersonic Nozzle, by Peter P. Wegener, 
Jack E. Marte and Carl Thiele, J. Aeron. 
Sci., vol. 25, March 1958, p. 205. 


An Energy Principle for Hydromagnetic 
Stability Problems, by J. B. Bernstein, E. 
A. Frieman, M. D. Kruskal and R. M. 
Kulsrud, Proc. Roy. Soc., vol. 244A, Feb. 
25, 1958, pp. 17-29. 

Turbulent Heat Transfer through a 
Highly Cooled Partially Dissociated Boun- 
dary Layer, by P. H. Rose. R. F. Probstein 
and M. C. Adams, Avco Mfg. Corp., Avco 
Res. Lab., Res. Rep. 14, Jan. 1958, 52 pp. 

Tables of Radiation from High Tem- 
perature Air, by B. Kivel and K. Bailey, 
Avco Mfg. Corp., Avco Res. Lab., Res. Rep. 
21, Dec. 1957, 29 pp. 

Heat Transfer Economics as Related to 
Design of Nuclear Fuel Elements, by D. 
Kallman, Chem. Engng. Progr., vol. 52, 
1957, pp. 289-292. 

A Method of Heating Matter of Low 
Density to Temperatures in the Range 10° 
to 10° K, by F. B. Knox, Australian J. 
Phys., vol. 10, June 1957, pp. 221-225. 

Prediction of Heat Transfer Burnout, 
by Louis Bernath, Chem. Engng. Progr., 
vol. 53, 1957, pp. 1-6. 

Pressure Drop During Forced-circula- 
tion Boiling, by Max Jakob, George Lep- 
pert and J. B. Reynolds, Chem. Engng. 
Progr., vol. 53, 1957, pp. 29-36. 

Effect of Agitation on the Critical Tem- 
ag Difference for a Boiling Liquid, by 

.S. Pramuk and J. W. Westwater, Chem. 
Engng. Progr., vol. 53, 1957, pp. 79-83. 

Some Experiments Related to the Noise 
from Boundary Layers, by Harvey H. 
Hubbard, J. Acoustical Soc. of America, 
vol. 29, 1957, pp. 331-334. 

Investigation of Burnout Heat Flux in 
Rectangular Channels at 2000 Psia, by H. 
S. Jacket, J. D. Roarty and J. E. Zerbe, 
Westinghouse Electric Corp. Atomic Power 
Div. (U.S. Atomic Energy Comm.), WAP D- 
A1W (IM)-3, Dec. 1955, 1 vol. 


An Empirical Correlation for Velocity 
Distribution of Turbulent Fluid Flow, by 
B. F. Ruth and H. H. Lang, J. Am. Inst. 
Chem. Engrs., vol. 3, 1957, pp. 117-120. 


Determination of Burnout Limits of 
Polyphenyl Coolants, by K. Sato, Aerojet- 
General Corp. (U.S. Atomic Energy 
Comm.), AGC-AE-32, Feb. 1957, 67 pp. 


Measurement of Local Heat Transfer 
Coefficients with Sodium-potassium Eutec- 
tic in Turbulent Flow, by K. D. Kuczen 
and T. R. Bump, Nuclear Sci. and Engng., 
vol. 2, 1957, pp. 181-198. 


The Heat Transfer and High-tempera- 
ture Properties of Liquid Alkali Metals, 
by I. I. Novikov, J. Nuclear Energy, vol. 4, 
1957, pp. 387-408. 

Thermal Conductivity of Some Re- 
fractory Materials, by J. F. Clements and 
J. Vyse, Trans. British Ceramics Soc., vol. 
56, 1957, pp. 296-308. 

Thermal Relaxation in Carbon Dioxide 
as a Function of Temperature, by F. 
Douglas Shields, J. Acoustical Soc. of 
America, vol. 29, 1957, pp. 450-454. 


Radiative Transfer in Discontinuous 
Media, by R. G. Giovanelli, Australian J. 
Phys., vol. 12, no. 2, June 1957, pp. 227- 
239. 

Rotational Waves in Turbulent Liquid, 
by R. V. L. Hartley, J. Acoustical Soc. of 
America, vol. 29, 1957, pp. 195-196. 


Material Transport in Turbulent Gas 
Streams; Radial Diffusion in Circular 
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Conduit, by S. Lynn, W. H. Corcoran and 
B. H. Sage, J. AIChE, vol. 3, 1957, pp. 11- 
15. 


Velocity Distribution in Vortices, by A. 
Timme, /ngenieur-Archiv, vol. 25, no. 3, 
1957, pp. 205-226 (in German). 

Ground Reflection of Jet Noise, by 
Walton L. Howes, NACA TN 4260, April 
1958, 56 pp. 

Near Field Jet Noise, by M. O. W. 
Wolfe, NATO, AGARD, Rep. 112, April- 
May 1957, 41 pp. 

Cloud-droplet Ingestion in Engine In- 
lets with Inlet Velocity Ratios of 1.0 and 
0.7, by Rinaldo J. Brun, NACA Rep. 
1317, 1957, 35 pp. (Supersedes TN 3593.) 

Experiments in Liquid Atomization by 
Air Streams, by Henry F. Hrubecky, J 
Appl. Phys., vol. 29, March 1958, pp. 
572-578. 

The Prospects for Magneto-aerodynam- 
ics, by E. L. Resler Jr. and W. R. Sears, 
J. Aeron. Sci., vol. 25, April 1958, pp. 235- 
245, 258. 

Note on Hydrogen as a Real-gas Driver 
for Shock Tubes, by Paul W. Huber, J. 
Aeron. Sci., vol. 25, April 1958 p. 269. 

Hypersonic Heat Transfer to Catalytic 
Surfaces, by Sinclaire M. Scala, J. Aeron. 
Sci., vol. 25, April 1958, pp. 273-275. 

A Note on Isentropic Compressible 
Flow of Air with Variable Specific Heats, 
by K. E. Tempelmeyer and L. Self, J 
Aeron. Sci., vol. 25, April 1958, pp. 278- 
279. 

Propagation of Weak Waves in a Dis- 
sociated Gas, by F. K. Moore, J. Aeron. 
Sci., vol. 25, April 1958, pp. 279-280 

Wall Effects in Shock Tube Flow, by 
Raymond J. Emrich and Donald B. 
Wheeler Jr., Physics of Fluids, vol. 1, 
Jan.-Feb. 1958, pp. 14-23. 

Effect of Radiation on Shock Wave 
Behavior, by R. E. Marshak, Physics of 
Fluids, vol. 1, Jan.-Feb. 1958, pp. 24-29. 

Hydromagnetic Stability of a Conduct- 
ing Fluid in a Circular Magnetic Field, by 
Frank N. Edmonds Jr., Physics of Fluids, 
vol. 1, Jan.-Feb. 1958, pp. 30-41. 

Statistical Behavior of a Reacting Mix- 
ture in Isotropic Turbulence, by Stanley 
Corrsin, Physics of Fluids, vol. 1, Jan.-Feb. 
1958, pp. 42-47. 

On the Diffusion of a Chemically Reac- 
tive Species in a Laminar Boundary Layer 
Flow, by Paul L. Chambré and Jonathan 
D. Young, Physics of Fluids, vol. 1, Jan.- 
Feb. 1958, pp. 48-54. 

An Energy Principle for Hydromagnetic 
Problems, by I. B. Bernstein, E. A. Frie- 
man, M. D. Kruskal and R.M. Kulsrud, 
Proc., Roy. Soc., vol. 244A, Feb. 25, 1958, 
pp. 17-40. 

Chemical Relaxation in Air, Oxygen and 
Nitrogen, by M. Camac, J. Camm, 8. 
Feldman, J. Keck and C. Petty, Inst. 
Aeron. Sci., Preprint 802, Jan. 1958, 21 
pp., 15 figs. 

The Interaction of a Reflected Shock 
Wave with the Boundary Layer in a Shock 
Tube, by Herman Mark, NACA Tech. 
Mem. 1418, March 1958, 128 pp. 


On the Mixing Theory of Crocco and 
Lees and Its Application to the Interaction 
of Shock Wave and Laminar Boundary 
Layer, Part II: Results and Discussion, 
by Sin I. Cheng and I. D. Chang, Princeton 
Univ., Dept. Aeron. Engng., Rep. 376 

(AFOSR-TN-58- 3; ASTIA AD 148042), 
Nov. 1957, 32 pp. 


Supersonic Mixing of Jets and Turbu- 
lent Boundary Layers, by Harry E. Bailey 
and Arnold M. Kuethe, Wright Air Dev. 
Center, Tech. Rep. 57-402 (ASTIA AD 
150992), June 1957, 45 pp. 


New Approach for the Calculation of 
Laminar and Turbulent Boundary Layers 
in Compressible Flow, Case of a Pressure 
Gradient; Peculiarities—Examples, by 
Alfred Walz, France, Ministere de I’ Air, 
Publications Scientifiques et Techniques 
336, 1957, 48 pp. (in French). 

Light Diffusion through High-speed 
Turbulent Boundary Layers, by Howard 
A. Stine and Warren Winovich, N ACA 
Res. Mem. A56B21, May 1956, 46 pp. 
(Declassified from Confidential by au- 
thority of NACA Res. Abstr. 125, p. 18, 
3/18/58.) 

Approximations for the Thermodynamic 
and Transport Properties of High-tem- 
perature Air, by C. Frederick Hansen, 
NACA TN 4150, March 1958, 67 pp. 

Incompressible Two Dimensional Stag- 
nation Point Flow of an Electrically Con- 
ducting Viscous Fluid in the Presence of a 
Magnetic Field, by Joseph L. Neuringer 
and William Mcllroy, Inst. Aeron. Sci., 
Preprint 764, Jan. 1958, 19 pp., 1 tab., 4 
figs. 


Charts for Flow Parameters of Helium 
at Hypersonic Speeds; Mach Number 10 
to 20, Princeton University, James Forres- 
tal Res. Center., Gas Dynamics Lab., 
Wright Air Dev. Center, TN 57-377 (AS- 
TIA AD 142310), Nov. 1957, 62 pp., 18 
figs. 

On Reducing Aerodynamic Heat Trans- 
fer Rates by Magnetohydrodynamic Tech- 
niques, by Rudolf C. Meyer, Jnst. Aeron. 
Sci., Preprint 816, Jan. 1958, 25 pp., 7figs. 

Approximate Three Dimensional Solu- 
tions for Transient Temperature in Shells 
of Revolution, by Maurice A. Brull and 
Jack R. Vinson, Inst. Aeron. Sci., Pre- 
print 775, Jan. 1958, 29 pp., 3 figs. 

The Hecker Method of Transient Tem- 
perature Calculation, by A. W. Trimpi, 
Inst. Aeron. Sci., Preprint 774, Jan. 1958, 
24 pp. 

On the Hydrodynamic Stability of 
Curved Laminar Compressible Flows, by 
M. Lessen, Inst. Aeron. Sci., Preprint 812, 
Jan. 1958, 16 pp., 5 figs. 

The Shock Wave in a Nonmonatomic 
Fluid, by Phrixos Theodorides, Maryland 
Univ., Inst. Fluid Dynam. and Appl. 
Math., TN BN-123 (AFOSR-T N-58-156; 
ASTIA AD 152182),Feb. 1958, 20 pp. 

Numerical Analysis of Inviscid Super- 
sonic Flows, Utilizing the Streamsurface 
Characteristics Method, by C. Cobb, H. 
G. Loos and S. Manus, Propulsion Res. 
Corp., Rep. R-279 (AFOSR-TN-58-105; 
ASTIA AD 152014), Jan. 1958, 31 pp. 

Critical Study of Certain Extensions in 
Rational Mechanics; Thermodynamics 
and Viscosity, by Henry Du Boseq De 
Beaumont, France, Ministere de I’ Air, 
Publications Scientifiques et Techniques, 
Notes Tech. 72, 1957, 46 pp. (in French). 


Experimental Study of Turbulent Flow 
in a Two Dimensional Variable Diffuser, 
by Jean Pierre Milliat, France, Minis- 
tere de l’ Air, Publications Scientifiques et 
Techniques 335, 1957, 134 pp. (in French), 

A Numerical Method for Calculating 
the Starting and Perturbation of a Two 
Dimensional Jet at Low Number, 
by R. B. Payne, Gt. Brit., Aeron. Res. 
Council, Rep. & Mem. 3407, (formerly 
ARC Tech. Rep. 18846), 1958, 50 pp. 

Boundaries of Supersonic Axisymmetric 
Free Jets, by Eugene S. Love, Mildred J. 
pag and Louise P. Lee, NACA Res. 


Mem. L56G18, Oct. 1956, 98 pp. (De-. 


classified from’ Confidential by authority 
of NACA Res. Abstr. 125, p. 24, 3/18/58.) 


An Analogue Computer for Convective 
Heating Problems, by H. G. R. Robinson, 
Gt. Brit., Aeron. Res. Council, Curr. Paper 


374 (formerly ARC Tech. Rep. 19098; 
Roy. Aircr. Estab., TN G. W. 434), 1957, 
15 pp., 10 figs. 

The Conservation Property of the Heat 
Equation on Riemannian Manifolds, b) 
Matthew P. Gaffney, Northwestern Univ. 
(AFOSR-TN-58-145; ASTIA AD 152- 
172), Feb. 1958, 12 pp. 

Heat Transfer in Turbulent Pipe Flow, 
by Raul R. Hunziker, J. Franklin Inst.. 
vol. 265, March 1958, pp. 205-225. 

On Time Smoothing, a Model for Mo- 
mentum Transport, by Marshall Fixman, 
J. Chem. Phys., vol. 28, March 1958, pp. 
397—400. 

Thermal Conductivities of Gases at Low 
Pressures, I: Monatomic Gases, Helium 
and Argon, by F. G. Waelbroeck and P. 
Zuckerbrodt, J. Chem. Phys., vol. 28, 
March 1958, pp. 523-524. 

Thermal Conductivities of Gases at Low 
Pressures, II: Rotational Relaxation 
Times in Hydrogen and Oxygen, by F. G 
Waelbroeck and P. Zuckerbrodt, J. Chem. 
Phys., vol. 28, March 1958, pp. 524-526. 

Drop-size Distributions for Impinging- 
jet Breakup in Airstreams Simulating the 
Velocity Conditions in Rocket Combustors, 
by Robert D. Ingebo, NACA TN 4222, 
March 1958, 23 pp. 

Thermodynamics of Electrically Con- 
ducting Fluids and Its Application to 
Magneto-Hydromechanics, by Boa-Teh 
Chu, Wright Air Dev. Center, TN 57-350 
(ASTIA AD 142039), Dec. 1957, 32 pp. 

The Shock Tunnel and Its Applications 
to Hypersonic Flight, by A. Hertzberg, 
NATO, Adv. Group for Aeron. Res. and 
Dev., Rep. 144, July 1957, 32 pp. 

Correlation of Turbulent Heat Transfer 
in a Tube for the Dissociating System 
N.0, 2NQ,, by R. S. Brokaw, NACA 
Res. Mem. E57K19a, March 1958, 17 pp. 

Effect of Jet Temperature on Jet-noise 
Generation, by Vern G. Rollin, NACA 
TN 4217, March 1958, 13 pp. 

A Theory about the Secondary Flow in 
Cascades, by Shintaro Otuka, Proc., 
Japan, 6th Natl. Congr. Appl. Mech., 1956 
(pub. 1957), pp. 327-332. 


Combustion, Fuels and 
Propellants 


Relation of the Turbojet and Ramjet 
Combustion Efficiency to Second-order 
Reaction Kinetics and Fundamental Flame 
Speed, by J. Howard Childs, Thaine W. 
Reynolds and Charles C. Graves, NACA 
Rep. 1334, 1957, 13 pp. 

Development of an Automatic System 
for Measuring Time in the Determination 
of Burning Rates of Propellants, by H. M. 
Burns and W. A. Hendricks, Hercules 
Powder Co., Allegany Ballistics Lab., Rep. 
ABL/B-16, Jan. 1958, 28 pp. 

The Equations of Motion in a Multicom- 
ponent Chemically Reacting Gas, by 
Sinclaire M. Scala, General Electric Co., 
Tech. Info. Series, Rep. R58SD205, Dec. 
1957, 38 pp. 

Irreversible Stochastic Thermodynam- 
ics and the Transport Phenomena in a 
Reacting Plasma, by H. J. Kaeppeler and 
G. Baumann, Forschungsinstitut fur Phy- 
sik der Strahlantriebe, E.V. Mitteilungen 8, 
1956, 104 pp. 

re and Thermodynamic Func- 
hae of the Chemically Reactive Combus- 
tion Gas: Hydrocarbon Air Mixture, by 
H. J. Kaeppeler and G. Baumann, For- 
schungsinstitut fur Physik der Strahlan- 
triebe, E. V. Mitteilungen 10, April 1957, 
221 pp. (in German). 
(Continued on page 864) 
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Servomechanisms’ proven ability to design and produce consi 
smaller and more reliable Central Data Computers, has establis 
SMI as the pace-setter in this exacting science. This experience, 
coupled with our major scientific advances in new materials research 
and deposited film circuitry techniques, will meet the requirements 
of even smaller and more reliable subsystems for the 

spacecraft of tomorrow. 


SUBSYSTEMS DIVISION, Hawthorne, California 
MECHATROL DIVISION, Westbury, L.!., New York 
SPECIAL PRODUCTS DIVISION, Hawthorne, California 
RESEARCH LABORATORY, Goleta, California 


GENERAL OFFICES « 12500 Aviation Boulevard, Hawthorne, California 
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By the application of properly chosen alternating and 
static electric fields, electrically charged particles can 
be maintained in dynamic equilibrium in a vacuum 
against interparticle and gravitational forces. This is 
illustrated in the above photograph of the orbit of a 
charged dust particle. During the time of exposure the 
particle traversed the closed orbit several times, yet it 
retraced its complicated path so accurately that its 
various passages can barely be distinguished. 

The range of particles of different charge-to-mass 
ratios which can be contained in this manner is deter- 
mined by the gradients of the static and alternating 
electric field intensities and by the frequencies of the 
latter. In the absence of static fields and for a given 
electric field strength, the minimum frequency required 
for stable containment of the particles is proportional 
to the square root of their charge-to-mass ratios. Thus, 
charged colloidal particles require the use of audio fre- 
quencies, atomic ions need HF frequencies, while elec- 
trons require the use of VHF and higher frequencies. 

Under the confining influence of the external fields, 


The Ramo-Wooldridge Corporation 


LOS ANGELES 45, CALIFORNIA 
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the particles are forced to vibrate with a lower fre- 
quency of motion which is determined by the external 
field intensities, space charge, and the driving fre- 
quencies. If the initial thermal energy is removed, a 
number of particles may be suspended in space in the 
form of a crystalline array which reflects the symmetry 
properties of the external electrodes. These “space 
crystals” can be repeatedly “melted” and re-formed by 
increasing and decreasing the effective electrical bind- 
ing force. These techniques offer a new approach in 
the study of plasma problems and mass spectroscopy in 
what may be properly termed “Electrohydrodynamics:’ 

At The Ramo-Wooldridge Corporation, work is in 
progress in this and other new and interesting fields. 
Scientists and engineers are invited to explore curren 
openings in Electronic Reconnaissance and Counter- 
measures; Microwave Techniques; Infrared; Analog 
and Digital Computers; Air Navigation and Traffic 
Control; Antisubmarine Warfare; Electronic Language 
Translation; Radio and Wireline Communication, and 
Basic Electronic Research. 
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REVOLUTION IN DESIG 


’ As revolutionary in design and 
engineering as the technological transition from flying 
machine to space vehicle... Achieved a quarter of a 
century ago, the only regulators which could be used 
in the early space labs of Annapolis...still today the 
standard for the critical control of high pressure fluids 
... Now on the threshold of space, Grove is uniquely 


* 
prepared—through advance research and development 


to meet the challenge of even higher pressures ahead. 


GROVE POWREACTOR DOME REGULATOR MODEL GH-408—50-6000 PSI INLET...5-3000 PSI REDUCED PRESSURE 


GROVE VALVE and REGULATOR COMPANY 


6529 Hollis St., Oakland 8, California + 2559 W. Olympic Bivd., Los Angeles 6, California 
Offices in other principal cities 
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SPECIALISTS 


Bell Aircraft Corporation offers you 
outstanding high level positions in 
the Propulsion Department of its 
Space Flight and Missiles Division. 
Experienced specialists are needed 
to assume the responsibility for con- 
ducting analytical and/or experimen- 
tal investigations in specified fields as 
related to advanced propulsion tech- 
nology. These challenging and unique 
positions offer you opportunities for 
technical advancement as well as pro- 
fessional and financial growth. 


AEROTHERMODYNAMIC SPECIALIST 
to work on supersonic aerodynamics, 
boundary layer flow, heat transfer and 
combustion problems. 5-7 years of ex- 
perience related to supersonic inlet and 
exhaust nozzle design including disso- 
ciation effects required. Ph D desirable. 


RAMJET SPECIALIST for preliminary 
design, evaluation and selection of ram- 
jet engines and components for ad- 
vanced missile programs. MS degree 
or equivalent plus 10 years progressive 
and intensive experience in all phases 
of ramjet engine research, design and 
development desired. 


TURBOJET SPECIALIST to assume 
responsibility for preliminary design, 
evaluation and selection of turbojet 
engines and components for advanced 
missile programs. MS degree or equiv- 
alent plus 10 years progressive and 
intensive experience in all phases of 
turbojet engine research design and 
development desired. 


THERMOCHEMICAL SPECIALIST to deal 
with transport properties of gases at 
high temperatures, dissociation and re- 
combination phenomena and rates, high 
temperature gas and surface interaction 
and boundary layer flows involving sub- 
limation or combustion. PhD plus 5-7 
years directly related experience desired. 


SOLID PROPELLANT ROCKET ENGINE 
SPECIALIST to engage in preliminary 
design, evaluation and selection of solid 
propellant rocket engines for advanced 
missile programs. MS degree or equiv- 
alent plus 10 years progressive and in- 
tensive experience in solid propellants, 

rain configuration, combustion and 

dware desired. 


Salaries are commensurate with 
your background. Living and 
working conditions are excellent 
and employee benefits outstanding. 
Write: 
SUPERVISOR, ENGINEERING 
EMPLOYMENT, DEPT. K-62 
BELL AIRCRAFT CORPORATION 


BUFFALO 5, NEW YORK 


Reaction Kinetics, Thermodynamics, 
and Transport in the Hydrogen-bromine 
System; A Survey of Properties for Flame 
Studies, by Edwin S. Campbell and Robert 
M. Fristrom, Chem. Rev., vol. 58, April 
1958, pp. 173-234. 

The Ignition of Combustible Gases by 
Flames, by H. G. Wolfhard and D. S. 
Burgess, Combustion and Flame, vol. 2, 
March 1958, pp. 3-12. 

Application of Well-stirred Reactor 
Theory to the Prediction of Combustor 
Performance, by H. C. Hottel, G. C. 
Williams and A. H. Bonnell, Combustion 
and Flame, vol. 2, March 1958, pp. 13-34, 

A Review of Some Unusual Stationary 
Flame Reactions, by W. G. Parker, Com- 
bustion and Flame, vol. 2, March 1958, pp. 
69-82. 

The Methyl Nitrite Decomposition 
Flame, by E. A. Arden and J. Powling, 
Combustion and Flame, vol. 2, March 
1958, pp. 55-68. 

Use of an Electro-optical Method to 
Determine Detonation Temperatures in 
High Explosive, by F. C. Gibson, M. L. 
Bowser, C. R. Sommers, F. H. Scott and 
C. M. Mason, J. Appl. Phys., vol. 29, 
April 1958, p. 628-632. 

Measurement of Rate Constants of 
Fast Reactions in a Supersonic Nozzle, by 
Peter P. Wegener, J. Chem. Phys., vol. 28, 
April 1958, pp. 724-725. 

Heats of Combustion of Some Organic 
Compounds Containing Chlorine, by G. C. 
Sinke and D. R. Stull, J. Phys. Chem., vol. 
62, April 1958, pp. 397-401. 

Cyanogen Derivatives as Rocket Pro- 
pellants, by Eckart W. Schmidt, Raketen- 
technik und Raumfahrtforschung, vol. 2, 
no. 1, 1958, pp. 20-21 (in German). 

Flow Visualisation Techniques Applied 
to Combustion Problems, by FE. F. Win- 
ter, J. Roy. Aeron. Soc., vol. 62, April 
1958, pp. 268-276. 

The Point of Oxygen Attack in the Com- 
bustion of Hydrocarbons, I: The Origin 
of Carbon Monoxide, by C. F. Cullis, F. 
R. F. Hardy and D. W. Turner, Proc., 
Roy. Soc., vol. 244A, April 22, 1958, pp. 
573-580. 

Reactions of Nitrogen Atoms, II: Hb, 
CO, NH;, NO, and NO,, by B. G. Kistia- 
kowsky and G. G. Volpi, J. Chem. Phys., 
vol. 28, April 1958, pp. 665-667. 


Materials of Construction 


A Vibration Manual for Engineers, by 
R. T. McGoldrick, 2d ed., David Taylor 
Model Basin, Rep. R-189, Dec. 1957, 27pp. 

New Graphites Beat Missile Hot Spots, 
by Irwin Stambler, Aviation Age, vol. 29, 
April 1958, pp. 150-154. 

Thin Pressurized Shells Look Best for 
Space Structures, by Joseph S. Lewin, 
Aviation Age, voi. 29, April 1958, pp. 86- 
92. 


Orbital Re-entry Will Intensify De- 
mands on Structures, by J. 8. Butz Jr., 
Aviation Week, vol. 68, April 21, 1958, pp. 
50-51, 53, 55, 57, 59. 

The Minimum Weight of a Structure 
Protected Against Short Duration Aero- 
dynamic Heating by Means of Thermal 
Diffusion, by Richard A. Dobbins, Jnst. 
A — Sci., Preprint 773, Jan. 1958, 20 pp., 
9 figs. 

Ignition of Kel-F and Teflon, by Lewis 


Greenspan, Rev. Sci. Instr., vol. 29, Feb. ° 


1958, pp. 172-173. 

Cermets as Potential Materials for High 
Temperature Service, by O. A. Sandven, 
NATO, Adv. Group. for Aeron. Res. and 


Dev., Rep. 101, April 1957, 25 pp 

Structures and Materials, by Paul E 
Sandorff, Aviation Age, vol. 28, March 
1958, pp. 50-63. 

High Temperatures Spur Use of Nickel- 
base Alloys, II, by T. E. Kihlgren, Avia- 
tion Age, vol. 28, March 1958, pp. 130-137. 

1958 Missile Materials Review, by 
Alfred J. Zaehringer and Raymond M. 
Nolan, Missiles and Rockets, vol. 3, March: 
1958, pp. 69-75. 

Aluminum for Missiles in Production, 
by Don Fabun, Missiles and Rockets, vol. 
3, March 1958, pp. 85-87. 

Materials Build a New Technology, by 
W. C. Rous Jr., Missiles and Rockets, vol. 
3, March 1958, pp. 91-92, 95-96, 98, 100. 

Welded Stainless Steel Hollow Core, by 
Dr. Michael Watter, Missiles and Rockets, 
vol. 3, March 1958, pp. 104-105, 107-110. 

High Temperature Brazing Looks Good 
for Missile Parts, by John V. Long and 
George D. Dremer, Aviation Age, vol. 29, 
May 1958, pp. 30-31, 33. 

Lithia Agents Boost Heat Resistance of 
Ceramic Coatings, by Paul A. Huppert, 
Aviation Age, vol. 29, May 1958, pp. 56- 
58. 

Ceramic Nozzles for Uncooled Rocket 
Motors, by W. L. Wroten, Raketentechnik 
und Raumfahrtforschung, vol. 2, no. 1, 
1958, pp. 21-22 (in German). 

Axial - temperature - gradient Bending 
Stresses in Tubes, by F. G. Hammitt, J. 
Appl. Mech., vol. 25, no. 1, March 1958, 
pp. 109-114. 

Thermal Stress Analysis of a Cylinder 
of Semi-plastic Material, by Donald 
Hunter and Glenn Murphy, Ames Lab. 
(U.S. AEC) ISC-839, Dec. 1956, 31 pp. 


Instrumention, Telemeter- 
ing, Data Recording 


Extending Transducer Transient Re- 
se by Electronic Compensation for 
High-Speed Physical Measurements, by 
F, F. Liu and T. W. 
Instrum., vol. 29, no. 1, 
pp. 14-22. 

Description of a Sensitive Micromanom- 
eter, by R. Eichhorn and T. F. Irvine 
Jr., Rev. Sci. Instr., vol. 29, no. 1, Jan. 1958, 
pp. 23-27. 

Industry Strives for Missile Telemetry 
Economy, by John Kinkle, Missiles and 
Rockets, vol. 3, no. 2, Feb. 1958, pp. 87-88. 

Industry’s Answer to the Data-process- 
ing Challenge, by H. H. Rosen, Missiles 
and Rockets, vol. 3, no. 2, Feb. 1958, 
pp. 95-102. 

Challenge to Industry: Spaceship Te- 
lemetry, by Henry P. Steier, Missiles and 
Rockets, vol. 3, no. 2, Feb. 1958, 
pp. 105-110. 

Vortex Tube Performance Data Sheets, 
by R. Westley, Cranfield, Coll. Aeron., 
Note 67, July 1957, 7 pp., 65 figs. 

Digital Computer Calculation of Trans- 
ducer Frequency Response from Its 
Response to a Step Function, by Ralph B. 
Bowersox and Joseph Carlson, Calif. Inst. 
Techn., Jet Propulsion Lab., Progr. Rep. 
20-331, July 1957, 12 pp. 

Information Criteria for Estimation of 
Telemetering Systems, by M. M. Bach- 
metiev and R. R. Vasiliev, Automatika 7 
Telemekanika, vol. 18, no. 4, 1957, pp. 
371-375 (in Russian). 

Static Transfer Device of a Pulse-fre- 
quency Telemetering System, by A. M. 
Pshenichnikov, Automatika i Telemekanika, 
vol. 18, no. 5, 1957, pp. 444-448 (in Rus- 
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Kodak reports on: 


why 5,000 $5 bills were stuffed into our till...a large optical device... 
the boys who do their arithmetic in advance 


The stickum of last resort 


For Eastman 910 Adhesive, of which 
the active ingredient is methyl-2- 
cyanoacrylate, we now have a slogan 
—‘‘the adhesive to try if no other 
will do.” 
: We are just being sensible. Its 
price by the ounce is $10; special 
pound price, $75. It does make pos- 
sible some distinctly advantageous 
new assembly techniques in a large 
number of industries. So we gather 
| from the correspondence incidental 
to the 5,000 orders filled during the 
past year. It bonds virtually every- 
thing (except silicones and polyole- 
fins, else how could we deliver it?). 
The biggest plus is the speed at 
which enormous bond strength is 
developed within minutes after 


application of this thin, clear liquid 
to one of the adherents. Within 
hours the tensile strength of its bond 
to steel is in the thousands of pounds 
per square inch. In the case of glass, 

_ rubber, or wood, bonds are stronger 
than the material itself. There is 
virtually no shrinkage on setting. 
No heating, no great squeezing, no 
evaporation is required. The bond, 
~ however, should not be depended 
‘ on for too many weeks at tempera- 
_ tures above 175°F, particularly in 
-- the presence of much moisture. _ 
That’s the minus. 
In case Eastman 910 Adhesive sounds — 
_ more interesting now than it did when we 
_ practically swamped our boat by offering 
samples at $5 an ounce, write to Eastman 
Chemical Products, Inc., Department 
E 910A, Kingsport, Tenn. (Subsidiary of 
Eastman Kodak Company). 


For blades (also vanes and 
buckets) 


This is a blade from a jet engine. 
Many mathematical minds, mighty 
mathematical machines, and much 
aerodynamical experimentation 
have created its shape. Violation of 
the plan to the extent of a few 
thousandths of an inch in a single 
cross-section of a single blade sucks 
at efficiency like a little leech. And 


there are so many blades in a single 
compressor or turbine that the total 
number of them made in the brief 
span of air-breathing non-recipro- 
cating history must compare with 
all the wooden spokes in all the 
wagon wheels of all the supply 
trains in all armies since Alexander 
the Great. Tolerances on wooden 
spokes have always been broad. 


Therefore we have been busy 
lately building this large optical 
device. 

It works like this. 


CONTOUR 


\LRRESNEL 
BI-PRISM 


+ 
o> 
CS) mercury varor Lamp 


Not long before this periodical 
reached its subscribers, the two 
mercury lamps were turned on and 
the first cross-section of the first 
blade was seen thus in magnifica- 
tion against its tolerance envelope 
scribed on the screen. Inspection 
from now on should go rapidly 
and well. 


The device has been named Kodak 
Section-Profile Projector, and it is 
enough to restore a man’s faith in the 
future of geometrical optics. Inquiries 
go to Eastman Kodak Company, Mili- 
tary and Special Products Sales, Roch- 
ester 4, N. Y. 


A fast film with a quiet name 


Dearly as a hunting-and-fishing edi- 
tor loves to describe some high 
Andean stream which only a major 
general has ever had occasion to 


This is another advertisement where Eastman Kodak Company 
probes at random for mutual interests and occasionally a little 
revenue from those whose work has something to do with science 
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whip, so the photographic press 
dotes on miracles in making pictures 
with a parsimony of radiant energy 
that few occasions demand. We 
would no more discourage this than 
we would the game of bridge. 

The fastest film that can be ob- 
tained by walking up to a film 
counter without previous notice and 
plunking down the money is labeled 
Kodak Royal-X Pan Film. It comes 
either as roll film or sheet film. It 
is grainy. The average photo-hobby- 
ist needs it about as badly as he 

~ needs a fishing trip to the Andes, 
except that a roll of the 120 size 
lists for only 75¢. 

On the other hand, take a fellow - 
whose hobby is bridge but whose > 
job is to record the path of a piece - 
of metal hurtling through the twi-— 
light a thousand miles above him. 
Or take his classmate who needs | 
to analyze the first microsecond of 
a thermonuclear explosion. These a) 
boys are not apt to mail in their 
negatives to us with a letter asking 
how come they turned out so dark. — 
They have been in close touch with — 
us before going out to their stations. — 

In their behalf we have very care- — 
fully considered all the data on the 7 
Royal-X Pan emulsion and con- 
cluded that for the present putting — 

it on 16mm, 35mm, or 70mm film 
would not be a good idea. Further- — 
more, the arithmetic indicates that — 
the most sensitive regular-produc- 
tion films in these widths, Kodak — 
Tri-X Film and Kodak Linagraph — 
Pan Film, might not be quite sen- _ 
sitive enough. 


No, this is not leading up to some new 
“Quadruple Grand Ducal Whoosh !-X 
Film.” In these cases of speed-above-all 
what we often recommend bears the 
quieter name, Kodak Spectroscopic 
Film, Type I-D(2). Jt is one of 111 
special emulsions we make for astro- 
physicists and the like. By some tests it 
is faster than Royal-X Pan, but we are 
not prepared to prove this. By advance 
arrangement through Eastman Kodak 
Company, Special Sensitized Products 
Division, Rochester 4, N.Y., as litile as 
100 feet of it in 35mm perforated form 
can be furnished for $10, or two 100-foot 
16mm rolls for $6.60 apiece (list). 


Prices stated are subject 
to change without notice. 
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Issue Pages 
January 1-72 

February 73-148 
March 149-212 
April 213-292 
May 293-364 
June 365-436 
July 437-504 
August 505-572 
September 573-640 
October 641-712 
November _ 713-788 
December 789-872 


Apams, Mac C., and ProssTEein, RONALD 


F. 

On the Validity of Continuum Theory 
for Satellite and Hypersonic Flight 
Problems at High Altitudes........ 

See also Prospstein, RONALD F. 

ADLER, ALFRED A. 
Calculation of Re-Entry Velocity Profile 
Acoston, Georce A., Woop, BERNARD J. 
and Wise, Henry 

Influence of Pressure on the Combustion 

ANDERSEN, W. H., Britis, K. W., DEKKER, 
A. O., MisHucx, E., Mog, G. and 
Scuuttz, R. D. 

The Gasification of Solid Ammonium 

ARENs, M. 
A Comparison of Turbojets and Ram- 
jets for High Speed Flight........ 
a Generalized Optimization Proce- 
dure for N-Staged Missiles........ 
AROESTE, and Rupin, Morton 
Relaxation Lags in Gas Flow.... 
Bajek, J. J., see Sears, G. A. 
Baker, Rosert M. L., Jr. 
Note on Interplanetary Navigation... . 
See clso SAMUEL 
Bartz, D. R. 

Factors Which Influence the Suitability 
of Liquid Propellants as Rocket Motor 
Regenerative Coolants............ 

See also Hastrup, R. C. 

BAUMANN, Rosert C. 
Design, Fabrication and Testing of the 


Benney, D. J. 
Escape From a Circular Orbit Using 
BENTON, MILDRED 
Artificial Satellites—A Bibliography of 
Recent Literature 


Bercer, J. 
Celestial Iconospherics, the 


Ultimate 


_ Bercer, W. J. and Ricuprro, J. R. 


Visibility of Orbital Points.......... 
Berwin, T. W., see Liv, F. F. 
Bevutter, F. J. and Ravucn, L. L. 
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Precision Measurement of Supersonic 

Rocket Sled Velocity—Part II...... 

Bitts, K. W., see ANDERSEN, W. H. 
Buioxsom, D. E., Jr 

Supersonic Aerodynamic Experiments 
Using Very High Temperature Air 

Use of Capacitor Discharges to Produce 
High Temperature, High Pressure Air 

Bocponorr, S. M., sce Vas, I. 
BoLLINGER, LOREN E. 
Evaluation of Flame Stability at High 
Reynolds 
Bonp, J. W., 
Plasma Physics and Hypersonic Flight 
Bracc, S. L. and Bywortn, S. P. Q 
Combustion Chamber Pressure Loss... . 
S. L. and Rartcirre, H. 
Characteristic Velocity as a Measure of 
Bren, R. L., see Dunvap, R. 
Broseck, W. M., CLEMENSEN, R. E. and 
VoRECK, w. 

A Recording Sodium-Line Reversal Py- 

rometer 
BroMBeErG, R. and Lipxis, R. P. 

Heat Transfer in Boundary Layers With 
Chemical Reactions Due to Mass 

See also Conen, C. B. 

Bryson, A. E., Jr. and Ross, STANLEY E. 

Optimum Rocket Trajectories With 
Aerodynamic Drag.............. 

BurcEss, D. S., see RicHMonp, J. K. 
BussarD, W. 

Concepts for Future Nuclear Rocket 

Bywortn, S. P. Q., see Brace, S. L. 
CaAMBEL, Att BULENT, see ZIEMER, 

RicHarp W. 

Cuanc, Curen-CHren and Hsu, CHENG- 

Tinc 

Solutions to Stoolman’s External Dif- 
fusion Equation for Instability of a 
Normal Shock Inlet Diffuser...... 

Curist, O. J. W. and SMa tt, B. B. 

Direct Digital Read-Out of Missile Role 

Foom: Film Records.............- 
Cau, 5.. T. 

An Iterative Method of Determining 
Equilibrium Compositions of React- 

CLEMENSEN, R. E., see Bropecx, W. M. 
Conen, C. B., BromBERG, R. and Lipxis, 

Boundary Layers With Chemical Reac- 
tions Due to Mass Addition........ 

Cotsourn, JosEPpH L., see ROTHBERG, 

SIDNEY 
CourTNEy, WELBy G. 

Ignition by Flow Over Hot Surfaces. 

Cox, 
Comment on “Stresses in Rocket Grains” 
Crocco, 

Comments on the Zucrow-Osborn Paper 

on Combustion Oscillations........ 
Crocco, Lu1c1, Grey, JERRY AND 

T. 

On the Importance of the Sensitive Time 
Lag in Longitudinal High-Frequency 
Rocket Combustion Instability... . . . 

Dawtey, RicHarp E., see GEDEON, 

Geza S 
Dean, Leo E. and SHurtey, Lucian A. 

Analysis 


of Regenerative Cooling 
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Rocket Thrust Chambers......... > 

Dexker, A. O., see ANDERSEN, W. H. 

Desster, A. J., Hanson, W. B. HERTz- 
BERG, M., McKrpsin, D. D. and 
Warictey, R. C. 

A New Instrument for Measuring At- 
mospheric Density and Temperature 
at Satellite Altitudes.......... 

Dietz, RoBert Jr. and HINNERS, 
ArTHuR H., 
Ramjet Test Facility Planning........ 
Dosrowo A. 

Satellite Orbit Perturbations Under a 
Continuous Radial Thrust of Small 

Dona.pson, W. F., see RicHMonpD, J. K. 
Duncan, D. B. 
Analysis of an Inertial Guidance System 
Dunvap, R., BrEHM, R. L. and NICHOLLS, 

A Preliminary Study of the Application 
of Steady-State Detonative Combus- 
tion to a Reaction Engine........ 

Dussourp, Jutes L. and SHAPIRO, 
AscHER H. 

A Deceleration Probe for Measuring 
Stagnation Pressure and Velocity of 
a Particle-Laden Gas Stream...... 

Ecxert, E. R. G., ScHNnewer, P. J., 
Haypay, A. A. and Larson, R. M. 

Mass- Transfer Cooling of a Laminar 
Boundary Layer by Injection of a 
Light-Weight Foreign Gas 

Evans, W., GIVEN, 
and Mutter, G. M. 
Ignition of Electrolytic Monopropellants 
by Submerged Electrical Discharge. . 
FARBER, MILTON 
Thermodynamics of Al2Os.......... 
Farre., E. C., see KiuGer, P. 
FINKELSTEIN, A. B., see Yuan, S. W. 
Friep, Burton D. 

Corrections to “Comments on the Pow- 

ered Flight Trajectory of a Satellite” 
GamBiLt, W. R. and Greene, N. D. 

A Preliminary Study of Vortex-Boiling 

Burnout Heat Fluxes.......... 
Garwin, L. 

Solar Sailing—A Practical Method of 

Propulsion Within the Solar System 
Gepeon, Geza S. and Dawtey, RIcHARD 
E 


Frank I, 


The Influence of the Launching Condi- 
tions on the Orbital Characteristics 
GERARD, GEORGE 
An Evaluation of Structural Sheet Ma- 
terials in Missile Applications. ..... 
Giept, W. H. and Ratt, D. L. 


Bore-Surface Temperature Variation 
During Rapid Firing of a 40-mm 

Gin, W 


Calculated Viscosity of a Solid Propel- 

lant Rocket Exhaust Gas Mixture... . 

Given, FRANK I., see Evans, MARJORIE 
W. 


GoETTELMAN, R. 
Measurement os Satellite Erosion Rates 
by the Backscattering of Beta-Rays 

M. 
The Optimization of Nozzle Area Ratio 
for Rockets Operating in a Vacuum 
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Visicorder Record—ectual size 


Wyle Laboratories in El Segundo, California, have 
used a battery of four Visicorder consoles like the 
one shown (right) to run a series of tests on a vital 
missile component. In the Wyle test project, the 
unique Visicorder consoles are easy to operate. Most 
parameters are low frequency, requiring response on 
the order of 5 to 60 cycles. 

The two calibrator control panels in each Visi- 


corder console accommodate 10 plug-in balance and Tom Jackson, Wyle engineer, examines Visicorder record 
matching units—designed to match tachometer gener- The HONEYWELL VisIcoRDER is the first high- 
ators, pressure transducers, thermocouples, expanded- frequency, high-sensitivity direct recording oscil- 
_ scale voltmeters, etc., to the Heiland galvanometers. lograph. In laboratories and in the field every- 
4 Dick Johnson, Instrumentation Branch Head at where, instantly-readable Visicorder records are 
Wyle, says, “This system, I feel, is one of the most pointing the way to new advances in product de- 
_ efficient instrumentation consoles in operation. Set-up sign, rocketry, computing, control, nucleonics 
- and calibration time has been reduced by the use of ...in apy field where high speed variables are 
_ _Visicorders by approximately 50%. This is due to the under study. 
simplicity of operation and trouble-free performance. To record high frequency variables—and 
There are no inking pens to clean, high-gain ampli- monitor them as they are recorded—use the Visi- 
fier maintenance, and so on, and we can also use corder Oscillograph. Call your nearest Minne- 
these consoles together to form systems of more than apolis- Honeywell Industrial Sales Office for a 
six channels.” demonstration. 


Reference Data: Write for Visicorder Bulletin 
Minneapolis-Honeywell Regulator Co., Industrial Products Group, Heiland Division, 5200 E. Evans Ave., Denver 22, Colorado 
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ound Support 


‘2 Continental's TC-106 turbine air compressor, developed in conjunc- 
tion with the United States Air Force, is now available for commer- 
cial jet liner ground support. The unit, supplying low pressure 
air, is especially suited to engine starting, cabin air conditioning 

and actuation of electrical generat- 

ing equipment for ground opera- 

tions of the aircraft . . . Continental 

Turbo-Compressors have been in 

operation with the Air Force for 

SNOW AND ICE more than four years, compiling an 
REMOVAL NOZZLE enviable service record for depend- 
ability and long life . . . Continen- 

tal’s new Snow and Ice Removal Nozzle unit, designed to be used 

oo with the TC-106 air compressor, reduces man hours and material 

a costs for ice removal from aircraft surfaces. Warm bleed air from 

__ the compressor is converted to a high-energy air shaft which erodes 

through the ice, lifting it from the surface and blowing it away. 

The Snow and Ice Removal Nozzle is equipped with a glycol spray 

attachment for quick application of anti-ice protection after the 
surface has been cleared. 


CONTINENTAL AVIATION & ENGINEERING CORPORATION 


12700 KERCHEVAL AVENUE, DETROIT 15, MICHIGAN 
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